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ABSTRACT 


Coodetiie pos 1 L 1 oning using range, integrated Doppler, and inter'* 
foromotrie observations from a eons to] 1 a cion of twenty-four Global 
Positioning System satellites is analysed, A summary of the proposals 
for geode fie positioning and baseline determination is given which 
Includes a description of measurement techniques and comments on rsiua 
deficiency and error sources. An analysis of variance comparison of 
range, Doppler, and interferometric time delay to determine their rela- 
tive geometric strength for baseline determination is included. An 
analytic examination of the effect of a priori constraints on posi- 
tioning using simultaneous observations from two stations is presented. 

Dynamic point positioning and baseline, determination using range 
and Doppler is examined in detail. Models for the error sources influ- 
encing dynamic positioning are developed. Included is a discussion of 
atomic clock stability, and range and bonnier observation error statis- 
tics based on random correlated atomic clock error are derived. 

Criteria for establishing observation schedules for optimum geometric 
strength for positioning solutions are. examined. Results of geodetic, 
positioning simulation studies are presented. 

Satellite, interferometry results based on the double differ- 
encing of simultaneous Interferometric phase measurements from two 
satellites are given. The effects of ephemeris and refraction errors 
and the. nonsimultaneity of observation are considered. 
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1. INTRODUCTION 


1 . 1 Background on a Global Positioning, System 

A Global Positioning System (GPS) is a passive all-weather 
navigation satellite system proposed for operation after 1985. The 
system uses the concept of passive satellite navigation based on highly 
accurate atomic frequency standards to enable the navigator to determine 
his three-dimensional position, velocity, and time Instantaneously on a 
continuous worldwide basis. Range and range-rate measurements taken 
simultaneously from four satellites will be reduced to determine these 
parameters [Milllken, 1978). A total of twenty-four satellites in three 
orbit planes will be available for navigation giving accuracies and 
availability far exceeding the current Navy Navigation Satellite System 
or Transit System [Stansell, 1978a] which GPS is designed to replace for 
navigation. With the number of satellites in view always exceeding the 
required number for navigation, the user may select a subset of four 
based on some criterion which optimizes the geometric strength of the 
navigation solution. 

The GPS system consists of three major segments: Space System 

Segment, Control System Segment, and User System Segment. Each segment 
is developed over three separate phases, each being a logical extension 
of the previous phase in an integrated and cohesive manner. 
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l?hasc X encompasses the Initial design and evaluation of system 
components Including the development of user equipment satisfying the 
various navigation applications [Borel, 1978], testing of user equipment 
at a ground based simulation facility [Denaro, 1978], and the space 
based system as satellites become available. These satellites are pro- 
totypes of operational satellites which will validate a new ranging 
technique and the stability of atomic frequency standards In a space 
environment [Bartholomew, 1978] . This Initial constellation will pro- 
vide four-ln~vlew geometry similar to the complete system for up to 
three hours each day over selected geographic areas. An Initial ground 
tracking network will be developed and tested during Phase X as a proto- 
type of the operational ground system [Russell, 1978]. Certain limited 
demonstrations of operational scenarios are to be conducted. 

Phase II consists of the Initial production of low cost user 
equipment and development of operational satellites. During this phase 
additional satellites will augment the Phase I constellation. This will 
result In a constellation of four satellites In each of three orbit 
planes providing eight hours of continuous four-ln-vlew geometry each 
day. These satellites will later be maneuvered to provide continuous 
worldwide two-dimensional navigation. 

Phase III builds upon this two-dimensional capability augmenting 
the constellation until a total system of twenty-four satellites in 
three orbit planes exists. Orbital periods are twelve hours. The 
ground tracking stations will become operational and modified as neces- 
sary to accommodate full system operation. 
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Summarizing, Phase 1 Is the concept validation period, Phase II 
is the system validation period, and Phase III consists of production 
and operation. Initial worldwide operational capability should become a 
reality after 1985. Phase I has been completed. 

The final. Space System Segment will consist of twenty-four 
satellites deployed In three orbit planes separated in right ascension 
by sixty degrees. Eight satellites are equally spaced within each 
plane. Integrated into each satellite will be at least two atomic fre- 
quency standards to maintain stable time and frequency required for pre- 
cise ranging. 

The Control System Segment is composed of a master control sta- 
tion, an upload station, and three monitor stations [Russell, 1978). 

The master control station and the upload station are currently located 
at Vandenberg Air Force Base in California and three monitor stations 
are located on Guam and in Alaska and Hawaii. These monitor stations 
measure the range and range-rate of the satellites, collect meteoro- 
logical data and forward this information to the master control station. 
Every monitor station is equipped with a cesium frequency standard. The 
master control station processes the data collected at the monitor sta- 
tions and its own tracking data to obtain best estimates of satellite 
ephemerldes and time synchronization offsets for the system. Predicted 
ephemerides and clock corrections are forwarded to the upload station 
for transmission to the satellite. 

The User Control Segment consists of the development and testing 
of electronic receivers and associated equipment required to perform 
navigation. The function of this equipment is to detect and to acquire 
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the GPS satellite navigation signals, to extract range and range-rate 
information, to perform corrections for ionospheric refraction, and to 
compute three-dimensional position and velocity and time, The expected 
positional accuracies of the system are nine meters in each horizontal 
component of position and ten meters in the vertical component ninety 
percent of the time. These estimates of accuracy are based on a single 
determination of position using four satellites based on the expected 
error budget and optimum satellite geometry [Mllllken, 1978], 

An eventual replacement of the Transit System by GPS would pos- 
sibly curtail geodetic positioning currently available with the former 
system using integrated Doppler observations and precise satellite 
ephemerides [Sims, 1972], At the present time Doppler positioning is 
playing an Increasingly important role in many countries for network 
denslflcatlon and control as detailed in the Proceedings of the First 
and Second International Geodetic Symposiums on Satellite Doppler Posi- 
tioning [1976, 1979], The curtailment of this program could have sig- 
nificant implications within the geodetic community. 

1 . 2 Review of Previous Studies 

A Global Positioning System, although designed for navigation, 
can offer the means for continued geodetic positioning using Doppler or 
range observations, Anderle and Tanenbaum [1974] point out that a GPS 
system is orders of magnitude better in oscillator stability and sup- 
pression of ionospheric refraction and is effected less by uncertain- 
ties in the gravity field. These factors imply that the typical errors 
present in current Transit positioning would be reduced using GPS. 
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In addition Che presence of six to nine satellites in view at all times 
means that continuous data acquisition will be possible as opposed to 
intermittent data obtained from Transit. 

However the extreme altitude of these satellites, having an 
orbital semi-major axis of over 25,000 kilometers, means that the rela- 
tive velocity or Doppler shift between a satellite and an electronic 
receiver on the earth would be smaller limiting the amount of posi- 
tioning information available from each integrated Doppler observation. 

A comparison of simulated range difference data from GPS and Transit 
demonstrates this geometric dilution of information. GPS range dif- 
ference data which have a maximum value of around 17 kilometers for a 
thirty-second integration period are approximately an order of magnitude 
smaller than typical Transit observations which can have a maximum range 
difference of 150 kilometers over the same integration Interval. 

Figures 1.2.1 and 1.2.2 Illustrate thirty-second integrated Doppler 
range differences for a typical Transit pass and for a high elevation 
GPS satellite pass respectively. The elevation angle of the satellite 
is given at the endpoints of the curve and at the time of closest 
approach (TCA) . In addition the maximum length of a GPS satellite pass 
is about six hours whereas a Transit satellite pass lasts about twenty 
minutes. Thus GPS range differences are smaller in magnitude than cur- 
rently obtainable Transit observations and, due to the length of a pass, 
range differences from consecutive Integration periods will vary less. 
This implies that continuous tracking of GPS satellites over a complete 
pass may not represent an optimum data acquisition procedure. A sequen- 
tial tracking approach in which a number of satellites are tracked over 
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Figure 1.2,1. NAVSAT Range Differences 
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Figure 1.2.2. GPS Range Differences 



segments of a pass may give a more geometrically significant collection 
of observations. 

Thun it la evident that GPS integrated Doppler observations 
offer certain real advantages over Transit observations but lack In 
geometric strength of observation. However the GPS system offers addic- 
tions! observational approaches) namely » ranging and, as will be dis- 
cussed below, the potential for interferometric observation. 

The majority of the investigations made to date have centered on 
the navigational capabilities of the OPS system. These studies consist 
of both simulations and analysis of actual observations to determine the 
accuracies achievable In numerous navigational applications . Donaro 
[1978] describes the initial testing of aircraft and land-based 
navigation receivers using the Inverted Test Range at *uma, Arizona, 
These tests Involved the use of ground-based transmitters simulating the 
satellite system. Stansell [1978b] considers the civil marine applica- 
tions of GPS and Cox [1978] describes the augmentation of an Inertial 
navigation system with GPS observations. Miller [1977] gives results of 
an analysis of ocean navigation using GPS range observations, and 
Kruczynskl [1978] considers aircraft navigation using a limited opera- 
tional phase of the GPS system. 

Numerous additional studies have centered on the theme of navi- 
gation using the GPS system. However only restively few studies have 
examined the possible geodetic or geophysical potential of this system. 
One of the earliest papers, given by Anderle [1978a], discusses the 
major error sources effecting GPS range and Doppler observations and 
arrives at anticipated accuracies for geodetic positioning and baseline 
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components by extrapolating results of a limited analysis based on 
single pass solutions for two components of position, Andcrle 11978b] 
again gives estimates of precision of relative station positioning based 
on GPS range observations. The results were again based on the projec- 
tion of limited results. Fell [1979] gives an Indication of the effect 
of atomic clock stability errors on positioning based on the use of 
range and Doppler observations obtained from one or two GPS satellite 
passes. These limited studies comprise the present results Indicating 
the potential of GPS range and Doppler observations for geodetic posi- 
tioning derived using a dynamic point positioning /ipprooch. 

In addition to dynamic positioning, Interferometric approaches 
have been proposed which utilize radio signals broadcast by GPS satel- 
lites to determine baseline components by measuring the time difference 
of arrival or phase of these signals at two stations. Counselman [1978] 
proposes to utilize interferometric observations derived from a series of 
continuous wave signals transmitted by equipment which would augment the 
GPS system satellites. Using this approach baselines ranging up to a few 
hundred kilometers would be measured. Counselman presents baseline 
uncertainty estimates for this syjtare based on the geometry of the 
satellite passes. These results are then adjusted to reflect the effect 
of unmodeled tropospheric refraction. Applications of the system are 
discussed. 

MacDoran [1979] proposes to derive interferometric observations 
from broadcast GPS satellite radio signals In a manner similar to that 
used in very long baseline interferometry [Dermanls, 1977] or In the 
portable ARIES system [MacDoran et al., 1978] both using 
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quasar sources » MocDornn gives a summary of the proposed SERIES system 
and estimates of the effects of random and systematic error sources, A 
graph of estimated baseline accuracy derived from SERIES is given. 

Finally , Bender [Letter to I. I. Mueller, 1979] proposes an 
interferometric approach In which the phases of the reconstructed OPS 
carrier frequencies with respect to a local oscillator are measured at 
two stations in order to monitor crustal movements. As with the pre- 
vious two interferometric proposals this approach remains in an early 
stage of development and the exact magnitudes of the error sources can 
only bo conjectured at present, A more detailed examination of all pro- 
posed systems of usage is presented in Chapter 2, 

1 * 3 Description of Present Study 

The major objective of this study is to present an analysis of 
geodetic positioning obtained from both dynamic point positioning using 
GPS range and integrated Doppler observations and from interferometric 
satellite observations . One of the basic aims of geodesy is the pre- 
cise and consistent determination of the coordinates of points of 
interest in an adopted earth-fixed frame of reference. How well this 
can be accomplished using GPS satellite observations will depend on 
many factors which must be examined in detail. 

The first; sfcep in this study, described in Chapter 2, is to 
examine the proposed methods for the geodetic implementation of Global 
Positioning System observations, These proposals are divided into two 
basic classes, dynamic positioning with range and Doppler observations 
based on the use of satellite ephemerldes and satellite interferometry, 
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A discussion of these techniques is presented giving the mathematical 
description of the observing technique. A brief discussion of rank 
deficiency Is presented for each system along with a discussion of the 
error sources effecting each. 

The second phase of the study Is a comparison of range , Doppler 
and Interferometric observations to determine their relative geometric 
strength for baseline component and chord length determinations. 

Hanging observations are treated In three distinct modes, as range, cor- 
related range difference and as interferometric observations. A 
description of the adjustment procedure Is given and an examination of 
the effect of a priori constraints on positioning using simultaneous 
observations from two stations is given for each approach. This analy- 
sis is presented in Chapter 3. 

Dynamic positioning using range and Doppler observations is 
addressed In Chapter 4. A detailed description of the error sources 
influencing dynamic positioning is presented and error models for these 
sources are developed. Included are a discussion of atomic clock error 
modeling and the development of the statistics for range and Doppler 
observation errors due to random atomic clock error. Ephemerls, atmos- 
pheric refraction and instrumental error sources are considered. Simu- 
lation of GPS range and Doppler observations is discussed along with 
criteria for the selection of satellites to be tracked which yield opti- 
mum geometric strength of solution. A sequential algorithm is derived 
for the estimation of geodetic station coordinates from range and 
Doppler observations with fully correlated weighting. Results of geode- 
tic positioning simulation studies are presented. 
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Satellite 1 interferometry re mil is are presented in Chapter 
based on the double d i f forone ing of inte.rforomet i to plume measurements 
from two satellites observed simultaneously at two locations. This 
observation procedure, is designed to eliminate, the effect of timing 
errors on the determination of baseline components. The effects of 
ophemoris and tropospheric refraction errors and the nonslmul tnnolty of 
observation are considered. 

A final summary and recommendations for additional analysis 
are presented in Chapter 6. 


SUMMARY AND CONSOLIDATION OF PROPOSED 


2 

SYSTEMS OF USAGE 


2,1 Introductory Remarks 

All currently proponed methods for the geodetic implementation 
of a Global Positioning System of navigation satellites have centered on 
the use of three basic types of measurement. These observations are 
range, integrated Doppler or range difference, and the Interferometric 
delay in time of reception or difference in phase of electromagnetic 
signals at two sites. Ranging and Doppler techniques discussed by 
Anderle [1978a] are suitable, for dynamic point positioning applications 
In which the coordinates of the tracking receiver are determined in an 
adopted earth-fixed frame of reference. Coordinate differences, or 
baseline components, may also be, obtained from such observations 
acquired at two or more stations. The interferometric approaches 
advanced by MacDoran [1979], Gounselman and Shapiro [1979], and Bender 
[Letter to I. I. Mueller, 1979] , although differing greatly in methodology , 
are. proposals for using the measured time delay or phase difference at two 
stations to determine baseline components in order to density existing 
geodetic control and to monitor crustal movements. 

In this chapter a discussion of these techniques is presented 
which summarizes each observational procedure and gives a mathematical 
description of the observation equations, A brief discussion of rank 


deficiency is presented for each system and the error sources effecting 
each are addressed. 

2.2 Dynamic Positioning Using Range 
and Doppler Observations 

The concept of dynamic point positioning using satellite obser- 
vations is nearly two decades old. At present geodetic point posi- 
tioning using integrated Doppler observations from Navy Navigation 
Satellites forming the Transit System is performed on a worldwide basis 
primarily for network denslfication and control. Stansell [1978a] and 
Laurlla [1975] give overviews of this system and its applications in 
geodesy and navigation. Although differing philosophies exist for the 
exact implementation of Doppler observations for geodetic positioning, 
as seen in the discussions of Brown [1976], Anderle [1974, 1976], and 
Colquitt [1979], where differences in methodology exist in such areas 
as parameter definition and procedures for treating Doppler observations 
either as uncorrelated range differences or as biased range, this system 
has made a great impact on geodesy. 

With a Global 'Positioning System of navigation satellites both 
range and Doppler observations are available for point positioning, 
although the electronic technology required to acquire these observa- 
tions differs greatly from current Doppler measurement methods . 

2.2.1 Measurement of Range and Doppler 

2. 2. 1.1 Range Measurement Procedure 

Each GPS satellite broadcasts on two L band frequencies, 

1575.4 MHz and 1227.6 MHz, called and L ^ respectively, to allow for 
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precise first: order ionospheric compensation. Modulated on the car- 
rier are two pseudo random noise (PRN) code sequences known us the pre- 
cision (P) code and the course acquisition (C/A) code. The P code is a 
binary random sequence generated at a rate of 10.23 me nab its each second 
and may be considered as a square wave whose frequency is 10.23 MHz and 
whose amplitude is randomly taken as plus or minus one every cycle 
depending on the code sequence. The C/A code is generated at a rate of 
1.023 megabits each second and may be considered as a square, wave simi- 
lar to the P code but having lower frequency. The C/A code repeats 
itself approximately every millisecond; whereas the P code lias a repeti- 
tion rate of approximately 38 weeks, although in practice the. sequence 
will be reset every week. 

Lindsay [1973] discusses the general properties of digital 
sequences known as pseudo random noise sequences for use in ranging 
applications. The desired properties of these sequences are: 

(i) the complete code cycle length must be long enough to 
avoid ambiguities in range, measurements; 

(11) the code symbol repetition rate must be high enough to 
obtain the required resolution of the range measurement; 

(ill) the autocorrelation function of the code should be simi- 
lar to that of band limited white noise having two distinct levels; 

(iv) to improve efficiency in radio frequency (IIP) transmis- 
sion the code should have a balanced number of ones and zeros over a 
complete period of the sequence so that the power of the modulated sig- 
nal is more evenly distributed about the carrier frequency. 
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The. significance of these properties will be apparent shortly. 

The signal transmitted by satellite 1 has the following form 
as given by Spllker [1978] due to the biphase modulation of the PPN 
codes and phase quadrature 


s /(t) 


a ApP^CtjD^tJcosCu^t+cj)) 


+ A^t^COsinCw^-hj)) 


( 2 . 2 . 1 ) 


where the carrier has the form 

L^(t) » cos(w^t+cj>) . (2.2.2) 

In equation (2.2.1) P'(t) Is a +1 pseudo random noise sequence. Thus 

whenever the P code changes sign the, phase of the cosine component is 

reversed by 180 degrees or biphase modulated. These phase shifts occur 

at the positive zero crossings of the L, carrier. The factor C("(t) has 

i 1 

an amplitude of plus or minus one and has the property that when the C/A 
code is minus one, the phase of the second term in equation (2.2.1) is 
reversed by 180 degrees. Thus the first and second term in that equa- 
tion will remain out of phase by 90 degrees or retain phase quadrature 

regardless of the code values. The factors A and A represent the 

p c 

amplitude of each signal when transmitted. The factor D^(t) is an addi- 
tional data code of amplitude +1, modulated on the carrier at a rate of 
50 bits per second, which gives the navigation message along with the 
information required to determine the time shift between the epoch of 
the received C/A code and the epoch of the received P code. Figure 
2.2.1 taken from [Butler, 1978] displays the biphase modulation of a 
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Figure 2.2.1. Time Domain Waveforms: (a) Unmodulated 

Carrier, (b) PRN Code Sequence, 

(c) Biphase Code-Modulated Carrier 
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Figure 2.2.2. Effect of PRN Code Modulation 
on the KF Spectrum 



carrier with a PRN code. The resulting RF signal and Its power are 
spread into a frequency Interval centered on the carrier whose distri- 
bution depends on both the bit rate of the code and on the code itself. 
Figure 2.2.2 also taken from [Butler, 1978] demonstrates this spread 
spectrum effect where f is the frequency of the, code. 

The I lr , signal is biphase modulated by either the P code or the 
C/A code. Assuming the former the L n signal has the form 

S T i (t) a B P." ( t) D ( t) cos (w 0 t+(j>) (2.2.3) 

P *1* 4- *+ 

where 

L 0 (t) a cos(w 0 t+ijO . (2.2.4) 

Both the L. and L 0 signals and all codes are in synchronization with 

X 

one another when generated. 

To measure range a ground receiver must generate the same PRN 
codes that are broadcast by the tracked satellite. This requires 
a priori knowledge of the codes selected for broadcast by each satel- 
lite during the current week. With the receiver generating the appro- 
priate P and C/A codes the range measurement is obtained by first 
shifting the C/A code in time, compensating electronically for the 
Doppler shift, until a maximum correlation with the received signal is 
obtained. Thus the C/A code is shifted In time by t'* and biphase modu- 
lated with the received signal giving 

(^(t-OSj. 1 (t) c A 5 C^(t-t , ')P^(t)D i (t) cos(w^t+<j>) 

1 (2.2.5) 

+ A c C^(t-t^)C^(t)D i (t)sin(w^t-hj)) . 
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Wien a maximum correlation of the C/A codes is reached the second term 
on the right side of equation (2.2.5) will have its power compressed 
into a much narrower band about the carrier frequency since the product 
C^(t-t'*)G^(t) is one, demodulating the signal. Since C^t-t") and 
P£(t) do not correlate the power of the first term is spread into an 
even wider band. Since the code correlation functions t(P,C/A) 
and <J>(G/A,C/A) are essentially two valued with a distinct maximum as 
discussed above, a value t' can be determined where maximum C/A code 
correlation occurs. 

Since the C/A code has a short period, t' may be multivalued, 
but a maximum correlation can be obtained readily. The data code D^(t) 
then provides the receiver with information relating the epoch of the 
broadcast C/A code to the epoch of the P code. Thus the approximate 
time required to shift the receiver generated P code to correlate with 
the broadcast P code modulated signal can be determined based on t'' and 
the data message. The P code correlation processes is performed until 
a maximum correlation occurs as in the C/A code correlation process. 

The unique time x for which P(t-x) correlates with the signal is the 

i . . 

measured quantity. By performing a second correlation on S_ (t) an 

u 2 

estimate of the first order ionospheric refraction may be obtained and 
applied to correct X as described in Section 4.1.3 and in [Spilker, 
1978] . The corrected value of x multiplied by the speed of light c is 
known as the pseudo range measurement. It represents the geometric 
range between the receiver and the transmitter plus the effect of the 
synchronization error between the receiver and satellite clocks. In 
addition the measurement is subject to other error sources 
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discussed below. For the moment, ignoring these error sources, the 
observation equation for pseudo range is 

R « cT ■=> Ip’ -pi + cAt 

s (2.2.6) 

« [(u - u) 2 + (v - v) 2 + (w - w) 2 ]^ 2 + cAt 

8 S 8 

where u , v , w are the coordinates of the satellite in an adopted 
s s s 

earth~fixed reference frame. The quantities u, v, w represent the 
receiver coordinates In the same frame and At represents the synchroni- 
zation error between the satellite and receiver clocks. 

2. 2.1.2 Doppler Measurement Procedure 

In the range measurement process both carriers are reconstructed 
since the C/A and P codes are correlated and biphase modulated with the 
received signal. In addition the data code is deciphered by the 
receiver and removed from the carrier. The result is a continuous wave 
carrier subject to Doppler shift. 

Two approaches may be taken to measure the accumulated Doppler 
shift over an Interval of time. First, in forming the range measurement 
the P code must be correlated with the received signal. Because of the 
relative motion of the satellite with respect to the receiver the sig- 
nal is subject to a varying Doppler shift and the electronic correlation 
process must time shift the receiver code at rates proportional to the 
range rate to maintain correlation. Thus a Doppler measurement can be 
obtained by monitoring the code sequence shift rates over an interval. 
The second procedure is to difference the reconstructed carrier f with 
a frequency generated by the receiver £q and count the zero crossings 
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of the resulting signal over a time interval. This second procedure In 
the standard Doppler technique In use currently, GPS Doppler receivers 
however could theoretically use either approach. 

In either case the observation equation for Integrated Doppler 
can be expressed as the range difference over the integration interval 
[t^tj]. The equation is 


AR 


|p B <V - p| - fito) - p| 


(2.2.7) 


“ * 


(£ 0~ V (t j " ^ 


where is the accumulated Doppler count over the interval. The 
measurement is subject to errors due to oscillator frequency variations 
and atmospheric refraction. As with range this measurement la made on 
two frequencies to allow for ionospheric refraction correction. 


2*2,2 Comments on Rank Deficiency of 
Range and Doppler Approaches 

Dynamic point positioning solutions are obtained from range and 
Doppler observations by linearizing equations (2.2,6) and (2,2,7) about 
an initial estimate of station and satellite position 

V *» AX + L (2.2.8) 


T 

and minimizing V PV with respect to the unknown parameters X. This 
minimization leads [Uotlla, 1967) to the least squares normal equations 

NX + U « 0 (2.2.9) 

where 

N *« A T FA (2.2.10) 
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and 

U « A T 1>L , (2. 2, IX) 

Because of a lack of coordinate system definition a unique solution to 
equation (2,2,9) is not possible since N is a singular matrix with rank 
less than the number of parameters. Despite the dynamical constraints 
imposed on satellite motion a unique solution to equation (2,2.9) can 
only be achieved if origin and orientation constraints are Imposed on 
the solution. In dynamic point positioning solutions these necessary 
constraints are usually imposed through the use of previously estimated 
satellite ephemerides. The satellite positions appearing in equation 
(2,2,6) and (2,2,7) are Included in the normal equations (2.2.9) with 
weighted constraints based on the accuracy estimates of the satellite 
ephemerldes utilized. If range or Doppler observations are made at two 
sites the station position solutions may be transformed into estimates 
of the parameters of the baseline connecting the sites. 

Arur [1977] performed a rank analysis of Doppler observations 
and found that the vector of coordinate differences between the 
observing station and the mid-arc state vector of the satellite pass, 
the velocity components of this vector, and the frequency offset (f^-f^ 
are estimable. The components of station position only become estimable 
if constraints are imposed on the ephemerls. For coobserving stations 
the Interstation coordinate differences are estimable quantities. A 
theoretical rank analysis carried out for ranging [Van Gelder, 1978] 
showed that the rank deficiency for the short arc mode is two. Thus 
without the use of sufficient constraints unique solutions to equation 
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(2.2.9) are not possible, PavllB [1979] discusses the general problem 
of rank deficiency and procedures for obtaining solutions. 

for positioning applications of the GPS system, satellite 
aphemerldes will be estimated based on ranging observations from four 
stations. The projected accuracy of these ophcmerldcs Is discussed in 
Section A. 1,2. Range and Doppler positioning studies described in 
Chapter A will Incorporate weight constraints based on assumed ephemerls 
accuracy . 


2.2.3 Range and Doppler Error Sources 
The accuracy of satellite ephemerldes and tropospheric 
refraction modeling and the stability of satellite and receiver atomic 
clocks will have Important consequences In the application of range 
and Doppler observations to geodetic positioning. An additional factor 
will be the precision of the electronic receiver. These sources of 
error are discussed in detail In Section 4.1. Their effect on geodetic 
positioning are discussed In Section A. 5. 

2 . 3 Satellite Interferometry 

Radio signals transmitted by GPS satellites have been proposed 
as a new resource for the application of Interferometric techniques to 
baseline determination. The interferometry technique Is based on 
observing the time (phase) difference of arrival of radio signals from 
a single source at two or more coobserving sites. Three different 
satellite Interferometry proposals have been advanced. MacDoran [1979] 
proposes to utilize the broadcast GPS spread spectrum signals by 
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cross-correlating the recorded signals at two sites as In very long 
baseline Interferometry, The observed quantity Is the time difference 
of arrival of the signal at the two sites subject to a time synchroni- 
zation error, A second proposal [Counselman, 1979] would derive Inter- 
ferometric phase observations from a series of continuous wave signals 
transmitted by equipment which would augment the GPS satellites, Obser- 
vations wouJd be made from at least four satellites simultaneously at 
each site to recover the components of the baseline in near real time. 
This technique relies on measuring the phase of up to ten continuous 
coherent signals broadcast from each satellite to eliminate the 2tt phase 
ambiguity which occurs when a continuous wave Is used. The phase 
measurements are differenced at both observing sites to form the Inter- 
ferometric phase difference. Bender [Letter to I, X. Mueller, 1979] 
proposes an alternative approach based on measurement of the phase of 
the reconstructed GPS carrier frequencies at two sites. The phase of 
the reconstructed carrier Is measured with respect to a signal based on 
the receivers local frequency standard. Bender proposes making such 
measurements from three or more satellites simultaneously or within a 
relatively short time Interval so that the local frequency standard 
stability is not a serious limitation. The use of a water vapor radio- 
meter Is proposed as In the MacDoran approach to virtually eliminate 
tropospheric refraction effects. This approach Is also subject to the 
2ir phase ambiguity which must be resolved. 

Thun three separate proposals have been advanced for an inter- 
ferometric determination of baselines. The first Is based on observing 
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the time difference of arrival of random signals at two sites. The 
ocher two proposals are based on the measurement of phase of continuous 
wave signals, 

2.3.1 Measurement of Interf e rometric 
Time Delay and Phase 

The Interferometric time delay Is the difference in the time of 
arrival of radio signals from a common source at two sites. In very 
long baseline Interferometry the sources are the extremely distant 
quasars. For the proposals described above the sources are radio 
signals emitted from GPS satellites. Using the notation of equation 
(2.2.6) the time difference of arrival at sites 1 and i is given by 

6t = ^ R i “\)/ C " (iP s ~ pj ~ Ip s “’PjJ)/® + (At^-At^,) . (2.2,12) 

In equation (2.2.12) the earth-fixed coordinates of the satellite 
appear since the radio signals received at each site are not incoming 
along parallel paths as with quasar sources. The last term 
in equation (2.2.12) Is the clock synchronization error of the two 
observing sites. 

If the observation is interferometric phase based on continuous 
wave radio signals, either broadcast or reconstructed, the observation 
equation has the form 

0 1 = T [ |P S "Pil " V + GAT i ] (2.2.13) 

where X is the wavelength of the signal, m^ is the integer number of 
wavelengths comprising the geometric range and At^ is the 


24 


synebronl nation error of equation (2*2.0). The wavelonfitH X Is a func- 
tion of time due to the Doppler ohlft. The difference in phase at two 
silt ea Is? Riven hy 

AO a 0 t - 0 ? [ !p‘ f; - p'J ~ |p 8 - “pj - (m 1 - ny)A + e(ATj - At^) ] . (2.2.1/*) 

.Since the Doppler shift in frequency will not be identical at both 
ohsevvinp. sites equation (2.2.1/*) is an approximation to the order of 
accuracy that the Doppler shift is known a priori. The third term in 
equation (2.2.1/*) is the 2ir ambiguity mentioned previously. Its 
a priori uncertainty will be a function of the initial accuracy of the 
observing station’s coordinates. 

Finally an examination of equations (2.2.12) and (2.2.1./i) shows 
that the. time delay and the difference in interferometric phase are 
related by 


fir 


AIL , 

2ir 


t m, 



(2.2.15) 


2.1, 2 Commen t' s on Rank Deficiency of 
Satellite Interferomet ry 

Equation (2.2.12) and (2.2.14) reveal that satellite interfero- 
metry observations are a function of satellite position unlike quasar 
observations and are related to the. difference in range between the 
satellite and the two observing sites. If normal equations for station 
position are formed from such observations the normal matrix N will 
not have full rank. Unless sufficient information is available on 
satellite position, a unique solution for earth-fixed station 


coordinates la not possible. Even with ouch constraints the normal 
oquationu can still tend to become singular no t*ho base! Ine distance 
decreases. Thin lo demonstrated In Chapter II, 

The following approximation [Counselmau, 1978] can he uoed In 
equations (2,2,12) and (2. 2. 14) to recast these equations In terms of 
b a a e 1 in e e o mp o nen t o 
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where n is the unit vector In the direction of the satellite, For 
h o 

short baselines defined as bavin?; 
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the second term in equation (2.2.16) may be deleted. Then equations 
(2.2.12) mul (2.2.14) become 

dT m (p^ “?(,) * P R /o (At ( - At^) (2.2.1b) 

and 

AO a -y- [ (p'j - ?(,) * - (»V { - m £ ) A + e (A t ;L - At ? ) 1 . (2,2.19) 

An examination of the derivatives of equation (2,2.18) with respect to 
baseline components and the time synchronisation error reveals that 
these parameters are estimable if four satellites are observed which do 
not lie on the same circle in the sky (Gounselman, 1978], 
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liquation (2.2.19) has the same form as equation (2. 2. IB) except for the 
.hi ambiguity term which cannot be separated from the synchronization 
error unless special procedures are implemented [Counselman, 1979]. The 
double differencing approach examined in Chapter 5 is another technique 
for handling this problem, 

2 . 3 . 3 Interferome tric Error S ources 

MacDornn [1979] and Counselman [1979] outline the systematic 
and random error sources effecting their proposals. Included are the 
frequency stability of the receiver clocks, transmission media errors 
consisting of tropospheric mul ionospheric refraction, GPS satellite 
positional accuracy and the precision of the instrumentation. 
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3. PRECISION COMPARISON OF RANGE, DOPPLER, AND 
INTERFEROMETRIC APPROACHES FOR 
BASELINE DETERMINATION 

A complete comparison of the positioning accuracies obtainable 
from range, Doppler, and interferometric satellite observations would be 
difficult to perform since proposals based on the latter approach remain 
in an early stage of development. For instance, the exact nature of the 
instrumental error sources associated with satellite interferometry can 
only be, conjectured at present. However range and Doppler geodetic 
receivers are currently being tested and estimates of measurement error 
are available. Therefore, for the range and Doppler proposals a 
detailed error analysis is presented in Chapter 4 . Then in Chapter 5 an 
Interferometric observation technique for baseline determination is con- 
sidered which has a distinct advantage over the range and Doppler 
approaches . 

In this chapter a comparison of the geometric strength of the 
three approaches is given based on the processing of range observations 
as range, correlated range difference and as interferometry or dif- 
ferenced ranges from two stations. This analysis will give an indica- 
tion of the relative geometric strength of each approach for the 
determination of coordinate differences and baseline distances using 
observations from a constellation of high altitude satellites. 
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1 . 1 Mathematical Model and Adjustment ^J’rojmluve 


3.1.1 Mathematical Mode l 

Let 11^ be the topocentrlc range from any ground station 
R l^ u i ,v l ,w l^ t0 an ^ satG -'--'‘- :I - te position ( ^ u j» v j ,w j^ 38 shorn in 
Figure 3.1.1, where the earth-fixed coordinate system (u, v,w) , Is 
oriented towards the Greenwich mean astronomical meridian (u-axis) and 
the Conventional International Origin (w~axls) with the v-axls forming a 
right-handed coordinate system with u and w, this coorulnate system 
being defined by the Bureau International de 1 Mleure (Bill). From 
Figure 3.1.1, the following equation can be written for the topocentrlc 
range 

R ij = [(u j“ U i )2 + (v j" V l )2 * (w j" w i )2ll/2 ' (3.1.1) 

From two consecutive topocentrlc ranges, R^ and R^, the range dif - 
ference is defined as 


AR , . . =11,. - R . . 
ijlc lie ij 


(3.1.2) 


and from simultaneous range observations, R, , and R.., taken at two 
stations 1> ; j_( u ^ ,v i ,w i^ and P L^ u L ,v Jl ,w Jl^ ’ Interferometric observation 
is defined as 


6R Uj ~ R ij " R Lj * 


(3.1.3) 
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3.1.2 Adjustment Procedure 

The mathematical models (3.1.1) through (3.1.3) of the general 

form 


L c F(X ) (3.1.4) 
a a 

may be linearized by a Taylor series expansion about preliminary values 
for station and satellite coordinates X q to obtain the observation cqua“ 
tions [Uotlla, 1967) 


V = AX + L (3.1.5) 

where V is the vector of observation residuals defined by 

V ° L - L, • (3.1.6) 

a d 
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The elements of L are the differences between the function F evaluated 
at the preliminary values for the coordinates and the observed quanti- 
ties 1,^ and A is a matrix of partial derivatives of P with respect to 
the coordinates. The vector X, representing corrections to the prelimi- 
nary coordinate values, will be estimated from the observations 
giving 


X = X + X 
a o 


(3.1.7) 


using least squares minimum variance estimation. 

For range observations R. . and R , made simultaneously at two 
stations the rows of the design matrix have the form 


and 


OR. 


A, . “ 




lj 0u^,0u^,3uj 

‘“n 1 0 Kj 1 


(3.1.8) 


where 


and 


OR 


X M. 


ij 0u jL ,0u j ,,0u^ 

[ ° I a £ j I “ a £j - 1 


±i 


T 

X o [du i ,dv i ,dw i ,du Ji ,dv jl ,dw Jl ,dUj ,dVj ,dw^3 . (3.1.11) 


u i ~ U 1 V i " V 1 w i ~ W 1 
R ij ’ R ij ’ R ij 


(3.1.9) 


(3.1.10) 


The index j ranges over the number of satellite positions where range 
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observations are acquired. For single station tracking the parameters 


du^, dv^, and dw^ are naturally omitted. 

For range difference observations AR,^ 
design matrix takes the form 


the 


contribution to the 


A ..... 8AI V _ 
ijlc 3u i ,3u k ,3Uj 

(3.1.12) 

ta lk ~ a ij I " a lk I a lj ^ 

where 

T 

X e [du i ,dv i ,dw 1 ,du k ,dv k ,dw k ,du^ ,dVj ,dw^] . (3.1,13) 


And finally for interferometric observations iSR,^ from two stations 
the contribution to the design matrix for each observation Is 


with 


i a ^ 1A.J 

i&j 9u i ,3l V 3u j 

[a ij I " I " a lj + a Jlj ] 

T 

X » [du i ,dv i ,dw i ,du^,dv^,dw Ji ,du^ ,dv_j ,dw^] . 


(3.1.14) 


(3.1.15) 


In the analysis presented In this chapter, which is Intended to 
compare the geometric strength of these three observational approaches, 
the satellite ephemeris will be assumed known and excluded from the 
normal equations. 

The least squares minimum variance estimate of X based on a set 
of observations is obtained by minimizing the function 

<J> = V T PV - 2K T (AX+L-V) (3.1.16) 
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with respect to the unknowns V, K, and X. P Is the weight matrix, for 
the observation sot. After minimizing <}> and eliminating the unknowns K 
and V, the least squares estimate for X Is given by the solution of the 
normal equations 

NX + U Q 0 (3.1.17) 

where 

N ° A T PA (3.1.18) 

and 

U «=* A T PL . (3.1.19) 

Solving equation (3.1.17) gives 

X *=> N*”^TJ . (3.1.20) 

The covariance of the parameter estimates Is given by the inverse of the 
normal matrix provided P Is the inverse of the observation ccvaiiance matrix, 

£ « N _1 . (3.1.21) 

X 

For observations from two stations the uncertainties in the base- 
line components are obtained by the linear transformation 

E Ax = GE x G T (3.1.22) 

where 

G = [-1 I] (3.1.23) 

and 

AX T = [du^ - du^, dv^ - dv_^, dw^ - dw^] = [Au,Av,Aw] . (3.1,24) 
The uncertainty in the chord length d is given by the transformation 
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(3.1,25) 


where 


2 T '1' T 

a, » hx. ir n iigx gV 

u Ax X 


ii » 






(3.1.26) 


3.1.3 Weight Matrix 

In minimum variance least squares estimation the weight matrix 
is taken to be the inverse of the covariance matrix of the observational 
errors 

P » X7 1 , (3.1.27) 

h 

For statistically Independent range observations with constant variance 
the covariance matrix is given by 

X R = a 2 I . (3.1,28) 

The dimension of this matrix is equal to the number of observations 
acquired. 

For N independent range observations taken from a single sta- 
tion, whose statistics are given by equation (3.1.28), the least squares 
normal equations for station coordinate Improvement are 

<3 - 1 - 29 > 

where 

X T = [du,dv,dw] . (3.1.30) 


For (N-l) correlated range difference observations, defined as 
the difference between successive ranges, the least squares normal 


34 


equatlono can be directly obtained from Che matrix components In cqua 
tlon (3.1.29) by the transformations 


a ar " BA ll 

(3.1.31) 

l ar d bl r 

(3.1,32) 

s AW a bk„b t 
AR R 

(3.1.33) 


where the matrix B la defined by 


■1 1 0 ... 0 
0 *-1 1 . . . 0 


(3.1.34) 


L o 


o o 


-l 


(N-l x N) 


with the range observation covariance matrix given by equation (3.1.28). 
The weight matrix for correlated range difference observations is given 
by 


..**1 


,Ts“l 


1'.; = (BE n B x ) a -i. (BB a ) 

a 


T v -I 


AR AR, s ~~R~ ' “ Jl 


For unit variance equation (3.1.35) becomes 


(3.1.35) 




2 

•1 

0 

0 

.0 


-1 

2 

-1 


0 

-1 

2 


12-1 

0-12 


(3.1.36) 


•J(N-lXN-l) . 
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The normal equations lor (N-l) correlated range difference observations, 
according to equation (3.1.17), may be written as 


(A AR ? AR A AR )X * A AR r AU L AR ” ° (3.1.37) 

or, using equations (3.1.31) through (3,1.36), equation (3.1.37) 
becomes 


m m n’ m m m a 

(Ajirojin a ba r )x + a^t(bb a ) 


bl r » o 


(3.1.38) 


Consider now N independent range observations taken simul- 
taneously at two stations at times t^t^...,^ 

R T “ . (3.1.39) 


The least squares normal equations for the parameter sat 

X T *=> [ du^ , dv^ , dw_^ , du^ , dv^ , dw^ ] (3.1.40) 

are given by equation (3.1.29) with modifications to allow for the addi- 
tional set of parameters, Defining N Independent satellite interfero- 
metry observations as the difference between simultaneous ranges, the 
least squares normal equations for Interferometry can similarly be 
developed from the matrix components of equation (3.1.29) by the trans- 
formations 


SR “ “Si 

(3.1.41) 

'SR " ,1L R 

(3.1.42) 

T 

' C3 Mg M 

l 5R W V 

(3.1.43) 
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where 


M - [“II] (NX2N) • (3.1.44) 

The weight matrix for the statistically Independent Interferometric 
observations is given by 


° Eli a » ( 20 2 !)** 1 


1 §R oR 
“ ( 21 )* 


(3.1.45) 


for unit variance. The normal equations for Interferometry arc 

(A 6R P SR A SR )X + A 6R P SR Tj 6r ° 0 (3.1.46) 

or, using equations (3.1.41) through (3.1.45), are equivalent to 

(aJm T (>IM T )“SiA r )S + ^M T (MM T )" 1 ML r « 0 . (3.1.47) 

Thus the weight matrices for range difference and Interfero- 
metric observations are obtained using the same linear transformation 
matrices which convert the range observations to the alternative data 
form. The range difference observations are correlated since each suc- 
cessive range difference observation Is formed using a common range. 
This is reflected by the off-diagonal elements in matrix equation 
(3,1.36). Finally, It was shown that the range difference and Inter- 
ferometric normal equations are directly obtained from the range normal 
equations if the weight matrix is also modified accordingly. In equa- 
tion (3.1.38) ne modified weight matrix becomes B T (Bb' P ) ^B and in 
equation (3.1.47) it is M T (MM T )~ 1 M. 
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3.1.4 Effect of A Priori Constraints on Positioning 
Based on Simultaneous Observation 

In general, for simultaneous observations from two stations, the 
normal equations are developed for two sets of parameters, station coor- 
dinates and ephemeris variables. Although the analysis In this chapter 
assumes the latter to be known, it Is of Interest to examine the more 
general form of the normal equations to arrive at an understanding of 
what effect a priori Information on either the ephemeris variables or 
the coordinates of one station has on the variance of the coordinates of 
the second station and baseline components. This situation would natu- 
rally arise In network denslficatlon using any of the observation types 
considered herein. 


3. 1.4.1 Range and Doppler 

For simultaneous range and Integrated Doppler observations from 
two stations, the least squares normal equations for station coordinates 
and ephemeris parameters have the following form for measurements that 
are either uncorrelated or are correlated by errors at Individual 
tracking stations 


" N u 0 N is" 


I 

X 

I - 1 

i 


J 

i — { 

i 

0 N 22 N 2S 


X 2 

+ 

U 2 

-^Sl N S2 N SS- 


1 

CO 

X 

i 


1 

CO 

I 


(3.1.48) 


The covariance matrix for station coordinates based on the observations 
and on a priori knowledge of the first station's coordinates and the 
ephemeris parameters is 
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■1 


)i a 

'x 


N ll + 1, l" N lS (N SS + V \l ” N 1 S (N SS + V N s?, 


~^2S^SS * P S' J N S1 


N 22~ N 2S (N SS + V N S? 



,! n 

B 12 

-1 


v-i 

•r- i 

or 

t 

'hz 

A3 




» 




> 

f) r) 

i* ** 



“a 

•* 

{ hz_ 


where the parameter set X is defined by 
T 

X a [du l ,dv 1 ,dw 1 ,du 2 ,dv 2 ,dv^ 2 ] . 

The covariance matrix elements are given by 


^11 ” ^ B ll *" B 12 B 22 B 21^ 

" t»ll + P l “ N JS (N SS + V"Vl " N 1S< N SS + P S ) " ll, S2 
^ N 22" N 2S fN SS + V N S2^ ^S^SS^V N S1^ 

(J 22” B 22 + B 22 B 21 Q ll B 12 B 22 

= [N 22~ N 2S^ N SS + 1 ’s ) N S2~ N 2S (N SS + 3 V N S1 
(N ll + P l” N lS < ‘ N SS + P S ) N S1^ N 1S N SS N S2^ 


} l2 ra "% B 12 B 22 


aQ U N lS (N SS + P S ) N S2 [N 22" N 2S (N SS + P S ) N S2^ 


-1„ i“l 


Q 21 aQ I2 * 


The matrix Q ^ is the covariance matrix for the first station 1 


(3.1.49) 


(3.1.50) 


(3.1.51) 


(3.1.52) 

(3.1.53) 
s 
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coordinates T, v , and 
X 1 

Cion's coordinates X, 


Q 

2 


22 


Is 


Che covariance matrix for the second sta- 


The covariance matrix for coordinate differences, or baseline 
components, is obtained from equation (3.1, 49) by the linear transfor- 
mation 


'AX 


[- 11 ] 


Q 11 

^12 

Q 21 

Q 22 

'22 “ 

^12 


-x 


'21 


(3.1.54) 


In terms of equations (3.1.50) through (3.1.53), equation (3.1.54) 
becomes 



Q u + « 22 5* B 


+ B 2?2Al 

B 22 + ^ B 22 B 21 


22 B 21 Q 11 B 12 B 22 + Q 11 B 12 B 22 
+ I]Q 11 CB 22 B 21 + I]T • 


(3.1.55) 


This equation may also be written in the form 


2 AX “ ^ + X X 2 + £ X 1 B 12 B 22 + B 22 B 21 ) ‘X 1 


(3.1.56) 


Consider now the effect of a priori information on the covari- 
ance matrices given by equations (3.1.50), (3.1,51), and (3.1.56). The 
following cases are considered: 

Case (1). Ephemeris parameters constrained and no knowledge of 
station 1 coordinates (P g = P^ “ 0) . Under these assumptions equa- 
tions (3.1.49) through (3.1.56) reduce to the results 
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^r 


}] « N 

**v i 

XX 

>%, “ N 22 
-1 -1 

V n M 4* N . 

‘ J AX 11 w 22 ‘ 


(3.1.57) 

(3.1.58) 

(3.1.59) 


Case (11). Epliemerls parameters partially known and station 1 


coordinates constrained (0 < Pg < •*>, P^ “ <*>) . 


Z a 0 

X 1 

“X,, " ^ N 22 “ N 2S^ N SS + V N S2 ] 
£. 

S A X “ *- N 22 " N 2S (N SS + V 


(3.1.60) 

(3.1.61) 

(3.1.62) 


Case (111) . Epliemerls parameters constrained and station 1 


coordinates constrained (P n - <», P,. ~ <*>) . 

& x 


2 V = 0 
X 1 


T a N 
6 X 2 N 22 

-1 

y a n 
^ AX iN 22 


(3.1.63) 

(3.1.64) 

(3.1.65) 


Case (Iv) . Epliemerls parameters and station 1 coordinates are 
partially known (0 < P g < 03 , 0 < P^ < «) . 


^ " [N n + p i - V N ss + VSi - W V \ 2 

^ N 22 " 8 2S^ fI SS + V N S2^ H 2S^ N SS + 'W 


(3.1.66) 
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V [H 22- ,, 2 S < H SS + P S> 1Ns2 ’ 1 

+ [ N 22 -li 2S^ , SS + V N S2 ! N 2S^ N SS + V "sl^Xj 

hs^V P S > ' lN S2' N 22 ' W V V’" 


(3.1.67) 


X AX = E X 1 +E X 2 ’ , ’ E X i n i2 B 22 + B 22 B 21 S X 1 * 


(3.1.68) 


A comparison of these results indicates that the uncertainties 
in the coordinates u 2> v 2 , w 2 are equivalent in cases (1) and (ill) 
where the ephemerls was assumed known and that this uncertainty can be 
expected to Increase as the orbit uncertainty increases as in case (li) 
and Increase further as the uncertainty in station 1 also increases. In 
terms of eigenvalues of the covariance matrix, or parameter uncertain- 
ties, the following relationship can be established among the cases 


^(iil) _ ^ 
X 2 



< 

X 2 


< X 


(lv) 

X n 


(3.1.69) 


The uncertainties in the coordinate differences u 2 - u^, ~ v^, 

W 2 ~ W .l are ^^ cew ^ se a function of the assumed a priori information. 
Comparison of the results indicates first that if the coordinates of one 
of the observing stations are known, increasing the ephemeris uncer- 
tainty increases the uncertainty of the baseline components; and second, 
that if the ephemerls is known, increasing the uncertainty of the first 
station's coordinates also increases the baseline component uncertainty. 
The relationships among the baseline component covariance matrices are 
established to be 


, (Hi) 
A AX 


< X 


(li) 

AX 


(3.1.70) 
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and 



>(i) 

AX 


(3.1, 71) 


Tlu» most important result however is obtained by noting that for rela- 
tively close stations the submatrices N^g and Ngg in equation (3.1.48) 
are approximately equal. Thus the last term in aquation (3.1.68) of 
ease (iv) will be negative definite insuring that the covariance for 
coordinate differences will be smaller than the sum of the coordinate 
covariance matrices, as opposed to cose (i) when the ephemeris is con- 
strained. Thus in general 


> < 4 lv) + 4 iv) 


(3.1.72) 


This demonstrates how baseline, component determinations may bo obtained 
successfully in the presence of ephemeris errors which cause larger 
uncertainty In the coordinates themselves. 


3. 1.4. 2 Interferometry 

For satellite interferometry observations from two stations the 
least squares normal equations for station coordinates and ephemeris 
parameters have the form 


— - 




** *" 

N N N 

1N 11 12 a lS 


X 1 


U 1 

N 21 N 22 N 2S 


X 2 

+ 

U 2 

N N N 

_ SI S2 lN SS_ 






(3.1.73) 
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After algebraic elimination of the ephemcria parameters the covariance 
matrix for station coordinates including a priori information is 


N ll + P l - V N SS + V\l "is - N 1S <N SS + V'Siz 

2 a 

V ~1 —1 

A N -N (N +r ) 1 N -N (M +F ) l 

lN 21 "2S u, SS "SI 22 n 2S v "SS S' S2 


B, „ B_ „ 
11 12 


B 21 B 22 


°*11 Q 12 


^21 Q 22 


(3.1.74) 


where the covariance matrix elements are given by 
**11“ tl) ll _I! 12 li 22 1! 21 1 

" ' N ll + W N SS + V '\l ■ (N 12 - N 1S < N SS + V\ 2> (3 • 1 • 75) 

< N 22 - N 2S (N S S + - "iS^SS + V'S/ 1 * 1 


^22 “ B 22 + B 22 B 2X^11 B 12 B 22 

“ t N 22- ,) 2S (,! SS + V' 1[, S2- V'V (3 - 1 ' 76) 

<' f U + P 1 - H lS (H SS + ,, s rlN S i rl ( t, 21- ,, 2 S (,i S S + P s rlN Sl )T l' 1 


‘ 7 2l"^U )! 12“22 

“ _<) ll (N 12 -t! lS^ >l SS + r S ) I W (N 22"' , W N SS + 1 V ”s2 ) 


(3.1.77) 


Q 21 = Q 12 # 


(3.1.78) 


Again find are the covariance matrices for the station coordi- 


nates . 
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The ovarlance matrix for the baseline components Is obtained 
using the transformation equation (3,1.54) and has the same form as 
equation (3.1.56): 


S AX “ Q 11 + Q 22 ~ Q 12 “ Q 21 


E X X + h z + 2 X 1 B 12 B 22 + B 22 B 21 J: X 1 * 


(3.1.79) 


For the cases considered previously, equations (3.1.75) through (3.1.79) 
reduce to the following 

Case (1). Ephemerls parameters constrained and no knowledge of 

station 1 coordinates (P„ c P 1 <=> 0) . 

o X 


S X 1 = ^11 " N 12 N 22 N 21’ ! 

S X 2 “ ^22 ~ N 21 N 11 N 12' 1 

S AX “ + S X 2 + S X 1 N 12 N 22 + N 22 N 21 Z X 1 ' 


(3.1.80) 

(3.1.81) 

(3.1.82) 


Case (11) . Ephemerls parameters partially known and station 1 

coordinates constrained (0 < P < P = “) . 

S X 


E = 0 
X 1 

S X 2 = [N 22 " N 2S (N SS f V N S2'* 

£ M ■ [N 22 - Vss + P S ) ' ll, S2J' 1 • 


(3.1.83) 

(3.1.84) 

(3.1.85) 


Case (111) . Ephemerls parameters constrained and station 1 
coordinates constrained (P g =* °°, P^ = °°) . 
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y 


(3.1.86) 


(3.1.87) 


(3,1.88) 


Case (iv). Ephemorls parameters and station 1 coordinates are 
partially know (0 0 < P. < <») . 


2 V ■ (see equation (3.1.75)) 

X 1 

E Y » (see equation (3.1.76)) 

X 2 

S AX - \ + S X 2 + Yu"” + A ' 


(3.1.89) 


An examination of these cases reveals the following relation- 
ships In terms of the eigenvalues of the station covariance matrices 


.(111) .(11) .(iv) 

A v ^ A y A y 

X 2 X 2 2 


(3.1.90) 


.( 111 ) .( 1 ) 

Av s Av . 


(3.1.91) 


For baseline compor 2 nt determinations 


.(ill) ,(ii) 

a ax a ax 


(3.1.92) 


In cases (1) and (iv) , however, the equations reveal how interferometry 
is suited for the determination of baseline components for close sta- 


tions. Under those circumstances 
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12 K “ N 22 

(3.1.93) 

12 K “ B 22 * 

(3.1.94) 


Thus, the last two terms In equations (3.1.82) and (3.1.8Q) are nearly 
the additive Inverse of the sum of the first two terms yielding an 
acceptable baseline component covariance even In the presence of large 
ephemeris error. 

Interesting also, under conditions where approximations 
(3.1.93) and (3.1.94) are valid, are the results obtained In cases 
(1) and (Iv) for the coordinate covariances, equations (3.1.80), 
(3.1.81), (3.1.75), and (3.1.76). In these cases, even when the satel- 
lite ephemeris is known, the covariance matrix tends to be singular as 
the baseline distance decreases. For Interferometry, station coordinates 
are not estimable under these conditions; however, baseline components 
are. 

3.2 Comparison of Range, Doppler, and 
Satellite Interferometry 

In this section, a comparison of range, Integrated Doppler, and 
satellite Interferometry techniques for the determination of baseline 
components and chord distances Is described. The basic Intent of this 
analysis Is to compare the relative geometric strength of each techni- 
que and obtain a measure of how the results themselves vary under dif- 
fering circumstances of usage. The analysis Is based on statistically 
independent range observations of unit variance taken simultaneously 
from two sites using a single-channel receiver as shown in Figure 3.2.1. 
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Figure 3.2.1, Simultaneous Range Measurements 
From Two Ground Stations 


These observations are treated as range, correlated range difference, 
and as Interferometric observations, Observations are included In the 
analysis If the satellite elevation exceeds 10 degrees. An analysis of 
variance Is performed using least squares minimum variance estimation 
Incorporating the weight matrices of Section 3.1.3. The parameters are 
the corrections to the baseline components 

Au = du^ - du^ 

Av e dVy ~ dv^ (3.2.1) 

Aw « dw^ - dw^ 


and the chord distance d, defined as 

d - [(u £ - Ui ) 2 + (v^-v^ 2 + (w^-w^ 2 ] 172 . (3.2.2) 

No time synchronization parameters are included. 

The orbital elements used In this study are given in Table 
3.2.1. With a 24-satelllte constellation five to nine satellites are 
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In view of a station at all times. For simultaneous observations from 
two stations, the number In view decreases with an Increase In the 
separation distance. Simultaneous three*»statlon tracking was not con- 
sidered since each baseline would not be determined as well, especi- 
ally for stations of great separation where the number of satellites In 
common view Is less. 


TABLE 3.2.1. GLOBAL POSITIONING SYSTEM ORBITAL ELEMENTS 


EPOCH: 197S DAY 116.0 


SATELLITE 


1 

2 

3 

4 

5 

6 
7 
0 
9 

10 

11 

12 

13 

14 

15 

16 
17 
10 

19 

20 
21 
22 

23 

24 


a 26560km 

1 B 63* 

B 0.0 

UJ« 0* 


M 

O 



0* 

0* 

45 

O 

90 

0 

135 

O 

100 

0 

225 

0 

270 

O 

315 

0 

345 

120’ 

30 

120 

75 

120 

120 

120 

165 

120 

210 

1 20 

255 

120 

300 

120 

15 

240 • 

60 

240 

105 

240 

150 

240 

195 

240 

240 

240 

205 

240 

330 

240 


Two station groups are considered. The first is a mid-latitude 


group of three stations whose geodetic coordinates are given In Table 
3.2.2, The chord distances separating the station pairs 1001-1002 and 
1001-1003 are approximately 100 kilometers. The second group of sta- 
tions is the so-called "Iron Triangle" very long baseline 
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interferometry (VLB!) stations whose geodetic coordinates are given in 
Table 3.2.3 and whose chord distances are found in Table 3.2,4. The 
maximum baseline distance for this group is nearly 4000 kilometers, and 
the minimum is 1500 kilometers. 

Since variations in the tracking scenario are possible with 
multiple satellites in view, a criterion for satellite selection is 
adopted. Given the normal matrix N based on all prior observational 
data, the next satellite to be selected for observation will be the one 
whose observations, when Included with prior data, minimize the trace 
of the parameter covariance matrix. For baseline components this trace 
is the sum of the baseline parameter variances 

ir<! w ■ °L + * °L • <3 - 2 ' 3) 

For the chord the trace is the variance of the estimated chord length. 
For each type of observation these criteria are virtually independent 
since minimizing the trace of the baseline component covariance matrix 
does not guarantee that the chord length variance is a minimum. That 
will depend on the correlations between the baseline components. 

3.2.1 Short Baseline Comparison 

An analysis of variance study was made for the mid-latitude sta- 
tions with parameter sets consisting of the baseline components and 
chord length. The observation schedules for the two baselines con- 
sidered, the north-south baseline 1001-1002 and the east-west baseline 
1001-1003, were based on ranging measurements taken every five minutes. 
Range observations were processed as range, correlated range difference 
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TABI,i: 3.2,3. IRON TRIANGLE STATION COORDXN/TES 


STATION 

GEODETIC COORDINATES 

LATITUDE 

LONGITUDE 

HEIGHT (m) 

WESTFORD IWS) 

42° 

36' 

46.518" 

7.33° 30' 22.720" 

■B 

OWENS VALLEY IOV> 

37 

13 

53.287 

241 43 2.441 


FORT DAVIS (FD) 

30 

38 

44.924 




TABLE 3.2.4. BASELINE DISTANCES (kw) 


ws 

























and as Interferometry. The results are based on an equivalent number of 
observations of each type, The satellites selected for tracking were 
chosen using the criteria defined above which are a function of the 
observation type and the parameter set definition. Several observa- 
tion schedules were considered where the time allotted for simultane-' 
ously tracking each satellite was fined at one, two, or three hours. 

One tynical observation schedule Is given in Figure 3,2.2 where the 
satellite tracking Interval is three hours. The analysis of variance 
results for the mid-latitude station group are given In Tables 3,2.5 
through 3.2.7. The results are based on 24 hours of continuous obser- 
vation with Intermediate results given at either 8 or 9 hours . No 
a priori knowledge of the station coordinates was assumed in these 
results. The range observations were taken as statistically independent 
having unit variance or a one meter standard error. To obtain an esti- 
mate of the ultimate precision obtainable for a particular observation 
type the results found In the tables must be scaled by the ratio of the 
assumed standard error in centimeters of that observation type, con- 
verted to the uncertainty of an equivalent range observation, to the 100 
centimeter standard error used to obtain the results. For instance, if 
it were assumed that correlated range differences may be measured with 
a standard error of 10 centimeters then, noting the defining equation 
(3.1,3) for range differences and equation (3.1.35), the standard error 
of an equivalent range measurement would be 10 centimeters divided by 
the square root of two. The standard error of an equivalent range 
measurement is defined as that /alue which when utilized in equation 


(3.1.35) or (3.1.45) would yield the assumed standard error for range 
difference or Interferometric observations respectively. 
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Figure 3.2.2. Typical Observation Schedule for a Three-Hour 
Satellite Tracking Interval (Stations 1001 
and 1002, Range Observations) 


Using 24 hours of range observations the baseline components 
are determined with an uncertainty of approximately 15 centimeters with 
slight variation as a function of the time interval each satellite is 
tracked. The chord distance has a standard error of approximately 11.5 
centimeters and increases, but not more than 10 percent, as the tracking 
interval increases to three hours. This increase is due to an increase 
in the correlations between baseline parameters. There are no discern- 
ible trends due to the orientation of the baseline. 
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For correlated Doppler reservations based on one-hour tracking 
Intervals there la a six-fold Increase In the baseline component uncer- 
tainty compared to the range results as seen from Table 3,2,5. How- 
ever as the satellite t,. eking Interval Is increased to three hours this 
standard error decreases dramatically. The same Is true of the uncer- 
tainty in chord length. No variation in the Doppler results is seen on 
the basis of orientation except with chord length where the uncertainty 
in the length of the east-west chord remains significantly larger in 
all cases, ranging from a difference of 2 parts to 0.9 parts per million 
(ppm). A comparison of the best Doppler results from Table 3.2.7 with 
the corresponding range results from Tables 3.2.5 through 3.2.7 indi- 
cates that results obtained from range observations with a one meter 
standard error can be equivalently obtained using correlated range dif- 
ferences if the standard error of the latter observation type is 
approximately 49 centimeters. This result is obtained by determining 
the uncertainty of an equivalent range observation (35 centimeters) 
which when used in the range difference weighting equation (3.1.35) will 
scale the Doppler results of Table 3.2.7 to be equivalent to those of 
Table 3.2.5 obtained using range observations with a one meter standard 
error. To obtain an equivalent uncertainty in estimated chord length, 
correlated range difference would require a standard error of 54 centi- 
meters . 

Also of importance is the ratio of the uncertainty of the esti- 
mated parameters to the observation uncertainty. This ratio is obtained 
by dividing the parameter uncertainties found in Tables 3.2.5 through 
3.2.7 by the standard error of the appropriate measurement type. 
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The standard errors for range, range difference, and interferometry obser- 
vations are 100, 141, and 141 centimeters respectively which may be veri- 
fied using equations (3.1.33) and (3.1.43) assuming the one meter 
standard error for range. Based on 24 hours of observation this ratio 
is approximately 0.15 for range and 0.28 for correlated range difference 
considering the best results for the latter. For the chord length the 
ratios for range and range difference are 0.12 and 0.14 to 0.21, respec- 
tively. The last two ratios for range difference reflect the variation 
in the results in Table 3.2.7 for the two orientations. These ratios 
are of importance as scale factors which can be applied to assumed 
observational uncertainties to obtain estimates for parameter uncer- 
tainties. For Instance, if correlated range differences had a measure- 
ment uncertainty of X centimeters, the uncertainty in the derived 
baseline components would be approximately 0.28X centimeters instead of 
the approximately 40 centimeters as given in Table 3.2.7. 

With Interferometric observations the resulting uncertainties of the 
baseline components are approximately twice as large as the range obser- 
vation results after 24 hours. The uncertainty increases about 25 per- 
cent as the tracking interval Increases from one to three hours . The 
uncertainty in the chord length is about 2.5 times greater than the 
range-derived chord. The trace of the covariance matriv from inter- 
ferometry shows little variation with orientation but variation in the 
distribution of the uncertainty among the parameters exists:, The chord 
length uncertainty is nearly equivalent for the two orientations. To 
produce baseline component uncertainties equivalent to the range 
results, the standard error of interferometric observations would be 
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requited to be 71 centimeters. Again this result is obtained by deter- 
mining the uncertainty of an equivalent range measurement (51 centi- 
meters) which when used in the interferometry weighting equacion 
(3.1.45) will scale the Interferometry results of Table 3.2.5 to be 
equivalent with the range results of Table 3.2.5 based on a one meter 
standard error of observation. For comparable chord results with iiiter- 
ferometry a standard error of 54 centimeters or an equivalent range 
uncertainty of 38 centimeters would be necessary. The ratio of para- 
meter uncertainty to measurement uncertainty is approximately 0.21 for 
baseline components and 0.38 for the chord length. 

The covariance computations for the one-hour tracking interval 
were repeated to obtain a measure of how knowledge of one station's 
coordinates could improve the results. The expected change in the base- 
line component covariance is given by a comparison of equations (3.1.59) 
and (3.1.65), which predict a square root of two decrease in the coordi- 
nate difference uncertainty for range and Doppler observations and by a 
comparison of equations (3.1.82) and (3.1.88) for interferometric obser- 
vations. In the latter comparison the exact decrease in the uncer- 
tainty to be expected is not as obvious. An examination of Tables 
3.2.5 and 3.2.8 show in fact that the uncertainty of the coordinate dif- 
ferences and also of the chord length decrease by the square root of two 
for range and Doppler. For interferometry the baseline component uncer- 
tainties decrease by approximately the square root of three and the 
chord uncertainty by approximately the square root of seven. Notice that 
the uncertainty in the chord based on Interferometry with one station held 
fixed is equal, to the number of digits given, to the chord uncertainty 
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based on range with no a priori constraints. For shorter baselines 
these results will be nearly equal, and this can be shown mathematically 
using equations (3.1.43), (3.1.54), (3.1.59), and (3.1.88) keeping in 
mind the partial derivative equations (3.1.10) and (3.1.14). It can be 
concluded for network densiflcatlon that a priori knowledge of the 
coordinates of the existing control point has a greater Impact on 
interferometry than on range and Doppler. 

3.2.2 Long Baseline Comparison 

For the long baselines of the Iron Triangle (Table 3.2.3) a simi- 
lar analysis was performed to determine the relative geometric strength 
of each observation type for determining baseline components and chord 
length. The results are based on an observation schedule of simultane- 
ous observations taken from two stations every five minutes for a 
full day. Here, one- and two-hour satellite tracking intervals are 
examined. Parameters corresponding to each side of the triangle are 
determined using only observations from the two stations forming that 
side. This allows the greatest flexibility In satellite geometry. The 
results are given In Tables 3.2.9 and 3.2.10. 

Based on 24 hours of range observations with one-hour tracking 
Intervals, the uncertainties of the baseline components range from 13.5 
to 17.0 centimeters showing minor variation with triangle side despite 
the different orientations and lengths. Increasing the satellite 
tracking Interval to two hours produces a marginal Increase in these 
uncertainties. The chord length uncertainty also Increases slightly 
with side length but in terms of ppm decreases significantly. 
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Increasing the tracking interval to two hours produces a small increase 
in the standard error of the chord* 

For correlated range differences the uncertainties are again 
much larger than the range results after one day of observation. The 
results Improve relative to the range results as the tracking interval 
increases. The trace of the baseline covariance matrix shows more varia- 
tion with baseline length than range, and the chord uncertainty is about 
three times larger than the range case, comparing best results . 

Interferometry observations give baseline component results 
which are better than correlated range difference results; however, the 
uncertainty in the chord length can be determined better from range dif- 
ferences from longer satellite tracking Intervals. The uncertainties of 
the parameters tend to increase as the tracking interval is Increased 
and a pronounced Increase in parameter uncertainty is noticed as the 
baseline length Increases. 

For range difference observations to yield equivalent base- 
line component uncertainties to range observations after 24 hours the 
standard error of range difference observations would need to be 
reduced to approximately 41 centimeters. In that case the range dif- 
ference results given in Table 3.2.10 would be reduced to approximately 
the level of uncertainty given in Table 3.2.9 for range observations 
with a one meter standard error. For chord length a range dif- 
ference uncertainty of 46 centimeters would be necessary to achieve 

equivalent results with range. 
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For Interferometry, a statement of the required observational 
uncertainty necessary to produce range equivalent reou3 ts Is more com- 
plicated since the results based on Interferometry are more variable no 
a function of station separation. All 24-hour Interferometry results In 
Table 3.2.9 would be at least as good as the range results in that table 
if the Westford-Owens Valley results were equivalent. The parameter 
uncertainties are greatest for this baseline. For this to occur a 
measurement uncertainty of 41 centimeters would be required for equi- 
valent baseline component results and 30 centimeters for chord length. 

Ratio? of parameter uncertainty to observational uncertainty 
may likewise be developed from Tables 3.2.9 and 3.2,10, 

3.2.3 Summary 

Some general conclusions can be drawn from an examination of 
the results. For the observation types considered it Is evident that 
ranging measurements provide the best geometric strength of solution. 

The two other derived observation types, correlated range difference 
and interferometry, are. geometrically weaker although the results 
obtained from these latter procedures can be Improved upon by Increased 
observational precision. Correlated range difference observations give 
the best geometric strength of solution If observed satellites are 
tracked over longer time Intervals. With this type of tracking proce- 
dure both the baseline component and chord length uncertainties are 
minimized. For range and interferometric observations shorter satellite 
tracking Intervals produce the least uncertainty In the baseline 
parameters. Lengthening the tracking Interval for these observation 

types Increases the resulting parameter uncertainties. However the 
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rate of Increase Is smaller than the variation in Doppler results pro- 
duced by decreasing the tracking Interval. And finally the interfero- 
metry approach becomes geometrically weaker as the baseline length 
increases r.o become a more significant percentage of the distance to the 
satellite, although the relative error in parts per million decreases 
for the baselines considered. 

The analysis presented above considered the relative geometric 
strength of three observation types, two derived from basic ranging. 

The results were based on the assumptions that satellite positions in 
space were known and that the basic ranging measurements were subject 
to uncorrelated stationary random noise. Some additional comments con- 
cerning these assumptions are appropriate. If the ranging measurements 
are in addition subject to a receiver timing bias, then timing parameters 
would be required to augment the current parameter set, at least one for 
each observing station. Under these conditions the range and correlated 
range difference results would be approximately equivalent depending on 
the satellite geometries sampled, tracking interval adopted, and the 
a priori uncertainty of the timing parameters. The interferometry nor- 
mal equations will also have to include these parameters and the base- 
line parameter uncertainty will be increased. For close stations the 
effective error introduced into interferometry observations would be 
the difference in the timing error at each station. The effect of 
timing error on interferometry can be greatly reduced by tracking addi- 
tional satellites simultaneously as considered in Chapter 2. If the 
first assumption concerning the accuracy of satellite positions is 
violated, the resulting baseline parameter uncertainty will increase. 


This was shown analytically. For short baselines the effect of the 
ephemcris error will be minor for all three approaches. 

Finally a comment concerning Doppler observations is necessary. 
In the precision comparison study the derived range difference observe- 
tions were correlated since successive range dd ff erences were formed 
using a common range. If Doppler results were obtained from independ- 
ent Doppler counts over the same time intervals the correlations in the 
weight matrix would vanish and the resulting parameter uncertainties can 
be expected to increase. 
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4. DYNAMIC POSITIONING USING RANGE 


AND DOPPLER OBSERVATIONS 

In this chapter a study Is presented which attempts to determine 
the accuracy of dynamic point positioning using range and Doppler obser- 
vations from a constellation of twenty-four Global Positioning System 
satellites. Two positioning problems have been addressed. These are 
the determination of the geodetic coordinates of a station and the 
determination of baseline components for stations which lie 100 to 2000 
kilometers apart. An error analysis is performed to determine what 
effect various systematic and random modeling errors have on tracking 
station positions determined by a least squares adjustment using simu- 
lated observations. All results are based on the use of a single chan- 
nel, dual frequency, sequential receiver whereby only one satellite is 
tracked at a time on two frequencies to virtually eliminate ionospheric 
refraction. 

The observations analyzed consisted of range and integrated 
Doppler measurements. For both data types the assumption is made that 
the observations are subject to two random noise processes, namely 
uncorrelated white noise with a normal distribution and correlated error 
due to integrated fractional frequency errors in both satellite and 
.elver atomic oscillators. 




In order to perform rigorous simulations of dynamic point posi- 
tioning a complete adjustment model must be adopted. Model parameters 
Included In this study are receiver coordinates, polynomial clock error 
models for station clocks, satellite state vector components, and a 
polynomial satellite clock model for each pass. A priori weighting con- 
sistent with the model error levels introduced Is included for satellite 
ephomerls and clock parameters allowing station coordinates to be esti- 
mated. Since two sources of random error are present in the observn- 
t.ion-5 weight matrices used in the adjustment account for each random 
process, the complex one being correlated atomic clock error. An analy- 
tical method is developed to give the statistics of this random process. 
The procedure starts with either actual or models of the Allan variance 
for a particular oscillator or class of oscillators and develops the 
statistics of range and Integrated Doppler observations based on the 
two oscillators used in deriving the measurements. Statistics for 
residuals to polynomial clock models are then obtained by a transfor- 
mation. These residual statistics are incorporated into the adjustment 
weighting. 

To further define the adjustment procedure, several studies were 
performed and are described in this chapter. A study was made to deter- 
mine if it is possible to perform a sequential adjustment of the con- 
tinuously observed measurements. Since all observations based on the 
receiver clock are correlated through random atomic clock error, must 
all data be processed simultaneously using a fully correlated weight 
matrix or i an the measurements be divided into fully correlated blocks 
each with independent clock models requiring adjustment? Secondly, a 
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study was made to determine the time sp>* er which a given polynomial 
clock model might be adopted. Further, tests were conducted for range 
and Doppler tracking to determine the optimal selection of satellites 
which produce the least uncertainty in derived station coordinates, a 
selection which produces the best geometric strength. And finally, the 
use of two-body analytic partial derivatives for orbit improvement 
rather than rigorous numerically integrated partlals based on a spheri- 
cal harmonic gravity field complete through degree and order eight was 
examined . 

Results from numerous computer simulations of station posi- 
tioning are included to demonstrate the effect of the error sources and 
evaluate the full weight matrix concept. In general the results were 
computed for cases where observations are six second ranges smoothed 
over 300 second intervals and 60 second integrated Doppler observations 
aggregated over 300 seconds. 

4 . 1 Error Sources Influencing Dynamic Positioning 

In this section the dominant systematic and random error sources 
influencing dynamic point positioning using CPS range and Dopplar obser- 
vations are described in detail and error models for these sources are 
developed. Also included is background information on atomic clock fre- 
quency error characterization required for an understanding of the dis- 
cussion An Section 4.4. 

'Ahe error sources considered here are believed to be the domi- 
nate ones effecting dynamic positioning. They include: 
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(,i) systematic errors in tin* computed ephemet’Idos ol GPd 
satellites represented by periodic radial, a long- track, and 
out-of-plane errors and in addition by .1 quadratic along-- track error, 
Models f'or these errors! were chosen to produce aatellite position 
errors having t '10 same signature an errors! produced in simu.1 tat ions oi 
ephemerisi accuracy by invert i gators currently involved in coinput iny iTt 
navi national and post-fit ephemeridesi . The niroulatiomf performed Iv 
these invent I y,a tors were done to establish probable error levels; in 
future er-t mated orbits;; 

(Li) residual systematic bias; and drift in ('PS satellite, 
clocks; after ephemeris and clock error estimates are obtained in the 
orbit determination problem. The levels; oi these errors were extracted 
from the same references; as in (i) ; 

(iii) correlated random satellite clock errors due to the 
inadequacy of polynomial satellite clock models. These errors are based 
on Allan variance frequency stability models for the rubidium oscilla- 
tors on GPS satellites;; 

(iv) uncorrelated random noise from the tracking receiver. 
Nominal range measurement uncertainty is one meter for a six-second 
measurement and three centimeters for 60-second integrated Doppler; 

(v) systematic bias and drift of the receiver's cesium 
atomic clock; 

(vi) correlated random receiver clock error based on an Allan 
variance model for the re ver cesium frequency standard; 

(vii) systematic tropospheric refraction error equivalent to 
five percent of the tropospheric model predictions. 
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Residual second-order Ionospheric refraction errors were not Included 
since this error Is shown to be only a few millimeters at GPS fre- 
quencies. 

4.1.1 Atomic Clock Errors and Frequency Stabili ty 

A clock Is a device which counts the cycles of a periodic 
phenomenon and among the most stable clocks In use are the atomic 
clocks which form the basis for atomic time scales such as Interna- 
tional Atomic Time (TAI) . Atomic time is used primarily as a measure of 
time Interval and is based on the electromagnetic oscillations produced 
by quantum transitions within the atom. An excellent reference on time 
and frequency Is edited by Blair [1974]. 

Global Positioning System satellites will incorporate rubidium 
frequeitcy standards to provide short-term frequency stability for the 
navigator, and ground tracking receivers for geodetic utilization will 
be assumed to incorporate cesium frequency standards to insure good 
long-term stability. The precise definition of stability is found in 
Blair [1974]. Basically it is a measure, usually given statistically, 
of the random fluctuations in frequency which can occur in a clock’s 
oscillator over specified periods of time. For a given time interval a 
particular oscillator is considered best if the expected level of fre- 
quency fluctuation is a minimum in terms of the Allan variance defined 
below. 

This paragraph deals with the characterization of typical errors 
associated with atomic clock time scales and statistical measures of 
frequency stability. This information provides the general background 
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required for the discussion related to the development of observation 
statistics and their use in geodetic positioning studies. 

4. 1.1.1 Characterisation of Atomic Clock Errors 
Let Cj represent an Ideal clock whose oscillator frequency f ^ is 
constant. The period of this oscillator Is by definition 

Tj « 1/f j . (4.1.1) 

In (t-t ) seconds of Ideal time cycles are counted and the time 
registered by the clock Is 

NjT “ Nj/fj (4.1.2) 

where is given by the integral 

t 

Nx - / fxdT = f (t- t Q ) . (4.1.3) 

S 

Thus the time elapsed from t is 

N_T = f_(t-t )T = t - t . (4.1.4) 

I I 3. o I o 

Consider now a typical atomic clock C^ whose frequency is sub- 
ject to error. From t this clock has a frequency represented by the 
model 

f (t) = fx + Af + f(t- t Q ) + £(t) (4.1.5) 

where Af is a frequency bias, f is a drift in frequency, and f(t) are 
random fluctuations in frequency. The clock C^ records cycles in the 
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time interval [t ,t] where 
o 


N » / £,(l)dX 
fc o 


(4.1.6) 


f(t-t)^ t „ 


£ (t - t ) + Af (t - t ) + — + / f(t)dx . 


In addition the clocks and may not be synchronized at t Q intro~ 
duclng a time or phase error at t Q represented as Al^T^.. Each count 
is incorrectly assigned the period T^. giving at t Ideal the time t^ as 


Af 

t, =N,T t + AN T = (t- t ) +~ (t - t ) 
ill o I o' fj 

. 0 

f (t - t )“ t 

+ — 2f f- + -t / «T)dT + AVl • 

I x t 


(4.1.7) 


From a comparison of equations (4.1.4) and (4.1.7) the time error at t is 

2 


Af 


£(t - t ) “ , t ^ 

T,(t) = t, - t » (t- t ) + :• >- -■ - + AN T + f- f f(T)dT (4.1.8) 

i i it 

o 


or, rearranging terms and Introducing new notation, 


T i (t) “ \ D i (t " t 0 )2 + R l (t “ t o ) * T i ( - t a ) + * (t) 


The quantity x(t) Is the random time error at t defined by 


(4.1.9) 


x(t) = — f f ( t) dx = / y (x)dx 
It t 


(4.1.10) 


where y(t) is the random fractional frequency error of oscillator i. 
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Consider the quantity T^(t) written as 


D, 9 e 

T.(t) = T,(t ) + R,(t~ t) + (t- t V + / y(T)dT . (4.1.11) 

X 10X04 0 . 

•=0 

Suppose an estimate of T^Ct) was made at £ as shown in Figure 4.1.1 

based on available data taken prior to t . If a clock correction based 

s 

on this estimate was applied to the time scale, then at l the error 

T. (t ) is due to the error in the prior estimate T, (t ) which was 
1 o is 

applied to the time scale, the effect of fractional frequency errors 

over the interval [t ,t ] 

s o 


~ o 

x(t Q ) “ / y ( t) dr 

s 


(4 J .12) 


and, systematic contributions to the time error in the form of a time 
drift and ageing , the quadratic term in equation (4.1,11). The error 
Ti(t Q ) with no clock correction at t is approximately given by 


- Ti ( t g ) + R i (t o -t g ) + ^-t/* /° y(T)dT (4.1,13) 

t s 


since T (t ) is an estimated quantity, If the time scale is corrected 

u S 

/s 

at t^ then T^(t Q ) would be an estimated offset independent of the cur- 
rent oscillator random error y(t) for t greater than t Q . 

The error equation (4.1.9) is the model used to describe the 
types of error present in atomic time scales. The deterministic errors 
consist of bias, drift, and ageing terms modeled as a quadratic poly- 
nominal in time. The ageing term is usually not observable for clocks 
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whose long-term stability Is good such as cesium. The additional term 
In equation (4.1.9) represents the random time error due to the integra- 
tion of random fluctuations In frequency. The magnitude of this term 
depends on the interval of tim which has passed since the scale was reset 
or calibrated and on the stability of the clock. Table 4.1.1 lists the 
error terms associated with atomic clock time scales. 


TABLE 4,1.1. ATOMIC TIME SCALE ERROR TERMS 

DETERMINISTIC NOTATION 


TIME BIAS 
TIME DRIFT 
AGEING TERM 


T l ( l o> 
R ; 


RANDOM 


X (t) INTEGRATED FRACTIONAL FREQUENCY 


4, 1.1. 2 Frequency Stability Measurement 
and Characterization 

Hellwlg [1977] points out that "the characterization of the 
stability of a frequency standard Is usually the most Important informa- 
tion to the user especially to those interested in scientific measure- 
ments and In the evaluation and intercomparison of the most advanced 
devices (clocks)." Since the frequency stability of a standard depends 
on a variety of physical and electronic Influences both Internal and 
external to the standard, measurement and characterization of frequency 
stability are always given subject to constraints or environmental and 
operating conditions. In addition frequency stability depends on the 

exact measurement procedure used to determine stability, 
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Frequency stability characterization Is done In both the fre- 
quency and time domain. In the time domain a frequently used measure of 
stability Is the Allan variance or Its square root, In the frequency 
domain It Is the power spectral density. 


4, 1.1, 2,1 The Allan Varianc e. The Allan variance as a time 
domain measure of frequency stability is found especially useful in 
practice since It Is obtainable directly from experimental measurements 
of fractional frequency error y(t) and because It contains all Informa- 
tion on the second moments of the statistical distribution of fractional 
frequency error. The Allan variance Is defined as follows; let 
y o ,y l ,y 2’ ,y k >y k+ 1 ,y k+2* * ’ ' obse5:vecl fractional frequency errors 
separated by a repetition interval of T seconds. For each Integer N 
greater than or equal to two calculate y , from 


y m C N 


(nrt-l)N- 1 
2 

k a niN 


?k 


m “ 0,1,2, .. ,M . 


(4.1.14) 


This is an average over N consecutive values of y. . The Allan varl 





P3 


2M 


M - 1 
2 

m B 0 




(4.1.15) 


An examination of this equation reveals that the Allan variance for a 

particular sampling Interval NT Is the average two-sample variance of 

the y (N) , 

'm 
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For frequency standards the square root of the Allan variance In 
usually given in graphical form on a log-log oeale, For Individual 
clannoR of frequency atandarde model u for the Allan variance are used 
which portray general frequency stability eharacterlntico, Ilcllwlg 
[1975] given examples of such models for many oscillator types. Figure 
4.1.2 shows the typical form. In this form, cr^Cx) is the square root of 
the Allan variance for the sample interval x. The quantity Qg is called 
the flicker floor and x 1 , x 2 > are the break points of the plot. The 
constants associated with this figure are usually specified for each 
type of frequency standard. A comparison of 3uch information can facil- 
itate the selection of a frequency standard for a specific application. 

The stability characteristics shown in the three regions of 
Figure 4.1.2 are typically present in many Allan variance plots of 
specified oscillator performance. The first part, region X, reflects 
the fundamental noise properties of the standard. This behavior con- 
tinues with increased sampling time until a floor is reached corres- 
ponding to region IX. After X 2 the performance deteriorates with 
increased sampling time. Hellwig [1977] outlines the error sources 
corresponding to each portion of the graph. The magnitude and slope of 
each segment will depend on the particular category of standard. 

Figure 4.1.3 details the performance specifications fur the 
Allan variance for the GPS satellite rubidium oscillator and for the 
cesium oscillator used in tracking receivers supporting orbit determi- 
nation. Tills latter oscillator Is an example of the type which will be 
used in range and Doppler geodet'-' receivers. 
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Figure 4.1.3. Allan Variance for Satellite 
ami Station Oscillators 


4,1.1 .2.2 Power Spectral Penalty . An alternative procedure for 
specifying the stability of a frequency standard, In the. frequency 
domain, is the use of the power spectral density (PSD) of instantaneous 
fractional frequency fluctuations y(t). Allan et nl. [1974] have Riven 
a useful model to represent the PSD for various categories of frequency 
standards. Tills model is in the form of a power low spectral density 
having the form 



(to) 


a 


a 


V2ir / 




0 < w < w, 


ill > O), 


h 


(4. 1.1(0 


where a takes on the integer powers between -2 and 2 inclusive depending 
on how the interval [0,in^] is to he divided into subintervals, one for 


83 , 


each a to be used. The quantity h^ Is a scaling constant and the PSD is 
assumed to be negligible beyond the frequency range [ 0 , 0 ^]. 

Barnes et al. [1971] and Meditch [1975] give the transformations 
between the time, domain measures of frequency stability in the form of 
the Allan variance and the power law spectral densities. Table 4.1.2 
taken from Mcdltr'.i gives these conversions for three types of frac- 
tional frequency error sources. 

4. 1.1. 3 Range and Doppler Observation Errors 
Due to Random Atomic Clock Error 

As previously discussed an atomic clock's time scale can be 
expected to differ from ideal time due to both deterministic and random 
errors. The. random component is due to integration of fractional fre- 
quency errors. A range observation determined by correlating the. PUN 
signal broadcast by the satellite with a similar signal generated in the 
receiver is subject to the random errors of both atomic frequency 
standards. The effective range error at time U due to the timing error 
in one of the time scales t^ is 

dl^Ct) = cl^Ct) (4.1.17) 

with the random component being the random walk 

t 

D (t) = c / y(T)dt (4.1.18) 

t 

s 

where c the velocity of light. The random component is due to the 

accumulaced effect of fractional frequency error since the clock's 

start or reset at t . 

s 
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The 4 random error (t) is correlated In time. Consider two 
measurements of ranpo UU j) and R(tj.) based on the use of the one 11 la- 
tor Ln the satellite, and assume momentarily that the receiver's oscil- 
Inter is free from random error. The covariance between there meamired 
ranp.es due to correlated fractional frequency error in the satellite 
oscillator is 


K[R(t. )R(t k >] - KCnCt j)n(t k )l 


t, 


“k 


- c **)•:[ f yCOdt f y(Odr"] 

t t. 


B 

t. t 


S 


lc 

n c ** / / E[y(T)y(T'') ^ 

t t 

8 S 


(4.1.1«) 


t. t 


k 


» C ~ / / $ (T-T‘')dTdT" 

t t n 
a s 


whei'e i is the autocorrelation function for fractional frequency 

yy 

error y(t) defined by 


^ yy (T-' r 0 a E[y(T)y(T"*) ] 

00 CO 

° f f yy“i:(y,y'V r jOdy dy‘" . 


(4.1.20) 


The function f (y,y'*,T,T'‘) is the joint probability density function for 
fractional frequency error. Here it is assumed that y(t) is a mean aero 
stationary random process. The function <i^(T-TO can be obtained by 


the Inverse Fourier transform of the given power spectral density 



Cm) * 


<I> (t) » ~ / S (w)e iWt: dw 

yy\w 2ir yy v ' 


(4.1.21) 


where 

t n t -* T'* . (4,1.22) 

A procedure for obtaining the autocorrelation function <I> (t) from the 

J J 

Allan variance Is given in Section 4.4. 

The variance of a range observation Is obtained from equation 
(4.1.19) by taking t^ equal to t^ 

t, t. 

al » c 2 / / 4» (T-TOdTdT-* . (i. .1.23) 

R j t t yy 

Allowing random frequency error In the receiver oscillator introduces 
additional, but similar, terms into equations (4.1.19) and (4.1.23) 
which must be considered when assessing the range uncertainty due to all 
random clock errors effecting the measurement. 

For Integrated Doppler or range difference observations the ran- 
dom measurement error associated with system clocks is the Integral of 
fractional frequency error over the Doppler integration interval. The 
random error In range difference due to one oscillator is 


(ft. 1 .24) 


n tj “ "'‘.I 5 " 

« e / y (i) tlx . 

4 

Notice in equation (A . 1 . 24) that the random error r)^,j la a function of 

t . , t., and v(t). The error does not depend on t . Range difference 
i ft 

measurements have the following correlation from euch oscillator 


<1 1 1 f r 


(4 . 1 . 25) 


« c“ f r T (T-T')dTdf 
t* t ■ V ‘' r 


with the variance 


O ‘1 J 1 

o?n a c.“ / / $ (t - T'')dTdx'' . (4.1.26) 

A « h yy 

Olnterve that the random range difference errors, whose statistics are 
Riven by equations (4.1.25) and (4.1.26), are, stationary; however, ran- 
dom range errors, whose statistics are given by equations (4.1.19) and 
(4.1.23), are not. A stationary random process is one whose statistics 
are invariant in time. 

For the oscillator performance specifications shown in Figure 
4.1.3 examples of the contribution to the. range error are given for both 
oscillators in Figures 4.1.4 and 4.1.5 over a five-day span. The clocks 
are assumed to be. perfect initially. Also included is the standard 
error for the random walk Tl(t) obtained using equation (4.1.23). The 
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procedure used in simulating the random range error is discussed in 
[Mcditch, 1975]. 


TABLE 4.1.2. ALLAN VARIANCE AND POWER SPECTRAL DENSITY 
FOR COMMON ERROR SOURCES 


ERROR SOURCE 

y (t) 

ALLAN VARIANCE 
a 2 y (?) 

TWO SIDED SPECTRAL DENSITY 

Syy (W) 

WHITE NOISE 

3l 

r 

No 

FLICKER NOISE 

2Ni In 2 

N, 


it 

iwi 

INTEGRAL OF 

N 2 1 

N 2 

WHITE NOISE 
(RANDOM WALK) 

3 

UI 2 


4.1.2 Ephemerls Error 

The ultimate accuracy of Global positioning System satellite 
ephemerldes and satellite clock solutions is difficult to predict since 
many factors influencing the final error budget have to be resolved. 
Among these are the number and location of tracking sites to support 
orbit determination, the exact estimation algorithm to be used including 
force modeling, and the final geometry of the satellite constellation. 

At present, errors in computed ephemerldes significantly exceed accu- 
racy design goals, especially in the prediction region used in naviga- 
tion as reported by Schalbly [1979]. 

In order to establish bounds on expected ephemeris and satellite 
clock errors simulations of orbit determination were performed assuming 
expected levels of model error for gravity, solar radiation pressure, 
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pole position, tracking site coordinates, and system clock errors. 

The results of these studies reported by Sehalbly [1976] indicate the 
following : 

(i) the radial component of position error is a twelve-hour periodic 
function whose amplitude ranges form one to three meters; 

(li) the along- track component of position error has two components, 

The first is a twelve-hour periodic function whose amplitude 
ranges from two to five meters. The second component is a quadra- 
tic function in time introducing maximum errors of up to twenty 
meters. In most coses this error appears to average five meters; 
(Hi) the cross-track component of orbit error is a twelve-hour periodic 
function whose amplitude ranges from seven to twenty meters. This 
error and the periodic radial orbit error appear to have zero 
mean; 

(iv) satellite clock solution errors have systematic components which 
may be modeled as a bias and drift; 

(v) ephemerls and satellite clock errors will be correlated in the 
sense that the net effect of all error sources on an observation 
residual will be smaller than the sum of the individual error 
sources . 

This analysis of expected orbital accuracy will be used as a basis for 
developing error models and a priori statistics for ephemerls state vec- 
tor components and satellite clock parameters in simulation studies 
designed to predict accuracies for dynamic point positioning. 
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4. 1.2,1 Epheincrio Error Model 

Using the results of the simulations described above, models 
for oplsomerls error can be developed for use in positioning studies. 
These models will consist of variations In the osculating orbital ele- 
ments for each satellite which will produce radial, a long- track, and 
cross-track orbit errors comparable with the simulation results. The 
magnitudes of these variations can be approximated by noting the errors 
introduced by changes In Keplerian orbital elements. For Instance, 
radial orbit error Is primarily a function of errors In the scmiroajor 
axis a and eccentricity e of the orbit. The model for radial orbit 
error will be developed as follows: 

Taking 

Aa « 6a c,ob(M + d) + 3 ) (4.1.27) 

and 

Ac « <$e , (4.1.28) 

where da and <Se are errors In a and e respectively and 3 fl is the phase 
of the error signal Aa, and differentiating the equation for the radius 
of a Keplerian orbit 


r = a[l-e cos (E) ] (4.1.29) 

with respect to a 

“ « 1 - e cos(E) k l(e« 0) (4.1.30) 

del 

one arrives at the error Introduced Into r by Aa 
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Ar « Aa . 


(4,1.31) 


Choosing An as In equation (4,1.27) gives periodic radial orbit error of 
amplitude 6a and phase 3 Q « The period of the error Is orbit period. 
Differentiating equation (4.1.29) with respect to e gives 

•j" u -a cos(E) (4,1,32) 

or 


Ar w -oAe cos(E) 

B -aAa cosCM+C') 


(4,1.33) 


for e and w approximately Kero, Introducing a phase error and using 
equation (4.1, 2-3) yields 


Ar £! -afie cos(H'>w + 3 ) * 

( 3 * 


(4.1.34) 


Thus an error In orbit eccentricity <Se Introduces a periodic radial 
orbit error whose amplitude is a<Se and whose phase Is 3 e - 

Along-track orbit error can be produced by variations In mean 
anomaly M and argument of perigee w as well as eccentricity Con- 
sidering M + w as a single element the tallowing model will be adopted 

A(M+w) «=f 6(M + w)cos(M + W + 3j,j^) • (4.1.35) 

For nearly circular orbits 

r a a[l-e cos(E)]aa (4.1.36) 

and the along-track error due to equation (4.1.35) Is approximately 
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As n aA(H + w) 


(4,1,37) 


« a5(M + w)eon(M + o) ) * 

This orbit period error lino nil amplitude of o<S(M+uj) and phaoe 

An error In eccentricity will aioo Introduce a periodic error 
lu the along-trnck satellite position, Thlo error will have the form 

Aa a 2nSe aln(M + w + 6 e ) , (4,1.38) 

Notice that the along-track error due to So has twice the amplitude of 
Induced radial error due to Sc. 

An error in the ascending node P of the orbit plane produces an 
alonp,- track bias 


As ° a6$2eos(i) (4,1.39) 

independent of the in-plane satellite position. 

In addition a quadratic along-traek error polynomial will be 
Introduced to produce the quadratic error to fit found in the orbit- 
accuracy simulation studies. This error is developed through an error 
in mean anomaly M of the form 

Afl o (A + BM+CM 2 )/a (4.1.40) 

where A is an epoch error given in meters and B and G are determined to 
allow the quadratic error to be symmetric over a seven-day span; 


B ° -5040 C 
C ° 2A/ 6350400 . 


91 


(4.1.41) 


The cross-track or out-of-plane orbit error Is due to error In 
Inclination 1 and ascending node ft, The cross-track error due to an 
error In orbital Inclination has the form 

A(rxv) » a A 1 sin (M + w + 0^) (4,1,42) 

for a circular orbit. The error model chosen for Inclination Is 

A1 ® 51 . (4,1,43) 

The cross-track error due to a change in ascending node Is 

A(rxv) -aAnsin(l)cos(M + w+ (3^) (4,1,44) 

where Aft Is given by 

Aft * 6ft . (4.1,45) 

This error has orbit period and amplitude a6ftsin(l) , 

Notice that the two radial error signals are 180 degrees out of 
phase when 0 and 0 are zero, Therefore let 

0 » 0 + 180° (4.1.46) 

& G 

In order that the total radial error have the functional form found in 
th>; orbit simulation results, 

Since no along-track bias was present in the simulation results 
no error will be Introduced Into the ascending node in positioning simu- 
lations. As a result all eross-track error will be attributed to an 
error in inclination. 
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In addition the phase will be adjusted by 90 degrees to 
allow the along-track error due to Inclination and mean anomaly plus 
argument of perigee to be In phase. 

4.1.3 Refraction Errors 

4. 1.3.1 Tropospheric Refraction Error 

The nonlonlzed portion of the atmosphere slows the passage of 
electromagnetic signals introducing an error Into electromagnetic 
measurements such as range, Doppler and Interferometric phase. Compu- 
tations of receiver positions for geodetic control utilizing such 
measurements must Incorporate either refraction error modeling of suf- 
ficient accuracy to eliminate these atmospheric effects or Incorporate 
corrections based on radlometry measurements as suggested by MacDoran 
[1979]. 

Currently, the Hopfield model [Hopfield, 1969] is used exten- 
sively to correct for tropospheric refraction present in satellite 
observations of range and Doppler. This model requires a toiowledge of 
surface weather conditions at receiver sites to ensure proper scale. 
These weather observations are of surface pressure, temperature, and 
humidity. The error in this model is generally assumed to be less than 
five percent of the total refraction. 

Alternatively, a water vapor radiometer may be used to measure 
the tropospheric refraction (wet component) in the slant range direction 
to the satellite at the time of observation. MacDoran [1979] has 
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I ml ten ted that: Hitch mi approach may have an accuracy of 2 
mill lino tors, virtually eliminating thin error aourcc. 


cm to a low 


4. 1.2. 1.1 Refraction Mode VI up, Approach. The tropospheric 
refraction model adopted in this analysis is a modified vernion of the 
Uopfield model [Anderle, 1474] involving a change in the form of the 
quartic polynomial to allow more rapid calculation. 

The theoretical form of the first-order tropospheric refraction 
correction for range measurements t« Riven by 


where 


a at 


/ (n~ I)r(r 2 - kV 


■ 1/2 


dr 


(4,1.47) 


and 


k « r sin a 

H 


cot; 


w.i.m 




The vectors r and r represent the position veetors for the 

U Sul 

observing station and satellite, respectively. The aenith nnp,l e a Is 

measured from the ellipsoidal normal ]i through the station (see 

s 

Figure 4.1.6) 


The index of refraction n is computed using surface weather 
measurements. Given the centigrade temperature T, the surface pressure 
V in millibars, and the, relative himidity U, the index n is computed 


using the equations 


rSATUllll 


s^v signal path / 
'Vs. / 


A/ 
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Figure 4 .1.6* Geometry of Tropospheric Refraction 


n - l » M, + M„ 

X A< 


if r < r < r . 
s *■* — i 


1 2 *> 

1 r - r . 
s i 


0 otherwise 


where the surface refractivities and radii for the dry and wet 
components are given respectively by 


Ni » C,776 xlQ“')i^T K 

N,, « C.373)K/T,? 

r. “ r, + 40,1 + .149T (km) 

4. S 

r, » r +12.0 (km) 


The quantity T K represents the Kelvin equivalent of the Centigrade tem- 
perature T. The quantity E Is the water vapor pressure given by 

E= (H/100)exp{-37.2yif>5+ .213166T k - .000256908T i; 2 } . (4.1.56) 

For radio frequencies up t>» 1 . ,JHz tropospheric refraction Is not a 
function of frequency. 

For Integrated Doppler the refraction correction is taken as the 
difference In equation (4 .1.47) applied at the end times of the Integra- 
tion Interval. The magnitude of tropospheric refraction on range observa- 
tions Is given in Figure 4 .1 . 7 as a function of zenith angle. This error 
grows rapidly when the elevation angle of the satellite falls below ten 
degrees . 



ZENITH ANGLE (DEGREE) 

Figure 4.1.7. Tropospheric Refraction Profile for Range Observations 
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The tropospheric refraction model described above attempts to 
predict the refraction error as a function of surface weather condi- 
tions, How adequately this is done in practice is difficult to deter- 
mine especially for non-vertical measurements. Ilopfield [1972] has 
compared this model to values of tropospheric refraction computed using 
meterologlcal balloon data and found good agreement for zenith measure- 
ments with the contribution due to the wet component suffering the 
largest error. However the dry component is predicted accurately from 
surface pressure alone. 

An adopted technique in utilizing this model is to include in 
the mathematical model for the observation equation a scaling parameter 
as an unknown to be determined in the adjustment procedure with an 
a priori uncertainty. For range observations the mathematical model 
aquation (2.2.6) becomes 

R “ [(u - u) 2 + (v - v) 2 + (w -w) 2 ] 1/2 + cAT + (1 + CJ6R . (4.1.57) 

s s s xv 

Fell [1975] used such an approach for orbit determination using Doppler 
observations from Transit satellites. Although this procedure tends to 
weaken the normal equations it reduces the level of unmodeled error and 
improves the accuracy of the estimated quantities. 

4. 1.3. 2 Ionospheric Refraction Error 

The lonosophere, the charged portion of the atmosphere which 
extends above 3.00 kilometers, has a variable index of refraction which 
is a function of both the frequency of the passing electromagnetic sig- 
nal and the altitude along the signal path since the electron density 
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varies with location. The electron density distribution which deter- 
mines the refractive index at a particular frequency is quite variable 
with a dominate diurnal variation due to earth rotation and with long 
term variations with the solar cycle. 

Clynch [1979] reviews the second order expansion of the refrac- 
tion index n as a function of frequency 


n « 1 


AjNOi) A 2 N (h) 


(4.1.58) 


where the A. are constants and U is the electron density, a function of 
i 

height along the signal path. 

The total ionospheric refraction for range observations is the 
difference between the integral of the refractive index along the opti- 
cal path and the geometric range: 


(SR = /lids - /ds . 
o g 


(4.1.59) 


Using equation (4.1.58) in equation (4.1.59) gives to second order the 
ionospheric refraction as 


(SR = 



(4.1.60) 


where 


B 1 = A 1 • / ’ Nl ( h ) ds • 
o 


(4.1.61) 


If two known frequencies are transmitted from a satellite, f^ and f^, 
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then each signal Is refracted according to equation (4,1.60). These 
two equations raav be combined to eliminate the first order effect: 



/i' \ 2 


i- _ 

/f V 

R" ° 

(R+ <Sr^) “ (■jr* ) (R + <$R 2 ) 

\ *» / 

: 

, [ 1 ] 
\h) 


(4.1.62) 


where (R+iSR^) represents the range observed using frequency f^. This 
two frequem'v technique leaves a residual range error eR of 

B„ 


eR 


f 2 - 2 ’ 

^1 1 2 


(4.1.63) 


Glyneh [1979] gives examples of two frequency corrected residual range 
errors for simulated Navy Navigation Satellite passes. These iono- 
spheric errors were computed by ray tracing through an ionospheric 
model with a range of sunspot numbers considered. For elevation angles 
above ten degrees the upper bound on the residual range errors is five 
meters. Using equation (4.1.63) with CR equal to five maters and the 
150 MHz and 400 Mllz frequencies of the Transit system, an upper bound 
on B 0 Is obtained. Computing eR using the Global Positioning System 
satellite frequencies of 1227 Mllz and 1575 MHz gives an upper bound of 
4.8 millimeters for the residual range error. For observations above 
twenty degrees this residual error has an upper bound of 1.9 milli- 
meters. Since this error is small, residual ionospheric error was not 
Included in the positioning studies conducted in Section 4.5. 


4.1.4 Instrumental Error Sources 

Tracking receivers designed to measure range and integrated 
Doppler from GPS satellites introduce random measurement error in 


addition to the random clock error discussed In Section 4,1,1, The 
reasons for this error are jitter in the carrier tracking loops and In 
the code correlation process and electronic thermal noise. Specifica- 
tions of the statistical properties of this noise are usually given for 
each receiver type; and In practice estimates of these properties are 
usually obtained through examination of observation residuals. Unfor- 
tunately It may be difficult to completely separate these receiver 
noise effects from clock noise since oscillator errors are manifest In 
the residuals, 

Jorgensen [1978] attempts to predict the short-term quality of 
range observations by fitting a ninth degree polynomial to 70 minutes of 
six-second high quality range observations taken at the Hawaii and 
Vandenberg tracking stations from two satellites. In this procedure It 
Is assumed that the polynomial models all systematic trends In the 
observations. The residuals from this least squares fit appear as white 
noise. Jorgensen concludes from this investigation that two frequency 
corrected range observations of high quality are subject to 60 centi- 
meters of white noise. An extrapolation of Table 4.4.5 Indicates that 
these residuals are due almost entirely to receiver noise since the 
expected level of residual clock error Is much less than this magnitude. 

For 60 second dual frequency Doppler observations the Stanford 
Telecommunications [1978] specifications for Model 5007 NGR receiver 
indicate that the error due to jitter in the carrier tracking loops 
should not exceed 0.9 centimeters. Thermal noise may increase this 
level to at least one centimeter. 
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Baaed on them- result;? Table 4.1,3 giver, bounds on the range of 
revolver noise which can be expected from present receiver technology. 
These levels will be used in Section 4.5. 


TABLE 4.1.3. MAGNITUDE OF RECEIVER. WHITE NOISE 


DATA TYPE 

NOMINAL LEVEL 

OPTIMAL LEVEL 

RANGE (6 sec) 

1 m 

60 cm 

INTEGRATED 
DOPPLER (60 sec) 

3 cm 

1 cm 


4 , 2 Si mu l ati on of Ob servations , 

4.2,1 Range and Doppler 

Range and Doppler observations were simulated For the tracking 
stations of Table 4,2,1 over time intervals ranging from two to five 
days. The locations of the three station groups utilized are shown in 
Figures 4.2.1 through 4.2,3. Observations of topoeentrie range and 
range difference were based on satellite positions obtained from the 
numerical Integration of the satellite's equations of motion using a 
force model consisting of the WGS72 [Scppelin, 1974] geopotential coef- 
ficients to degree and order eight, solar radiation pressure, and 
luni-solar gravitational perturbations. The initial conditions for the 
orbit integrations were obtained from Table 3.2.1 for each of the 
twenty-four CPS satellites. The observation sets consisted of range and 
Doppler range differences generated every five minutes. Satellites were 
tracked sequentially and selected on the basis of criteria discussed 
below. The duration of the satellite tracking Interval varied from 
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TABLE A. 2.1. GEODETIC COORDINATES OF TRACKING STATIONS 


STATION 

NO. 

GEODETIC COORDINATES 

LATITUDE 

LONGITUDE 

HEIGHT (ml 

1001 

30' 

0' 

0.00" 

46° 

O’ 

0.00" 

0,0 

1002 

30 

64 

7.66 

46 

0 

0.00 

0.0 

1003 

30 

0 

0.00 

43 

67 

48.96 

0.0 

1004 

30 

0 

0.00 

46 

2 

11.04 

0.0 

1005 

30 

0 

0.00 

43 

26 

43,44 

0.0 

1006 

30 

0 

0.00 

46 

33 

16.56 

0.0 

1007 

30 

0 

0.00 

46 

24 

32,25 

0.0 

1008 

30 

0 

0.00 

47 

35 

27.76 

0,0 

1009 

90 

0 

0.00 

0 

0 

0.00 

0,0 

1010 

89 

6 

65.18 

0 

0 

0,00 

0.0 

1011 

0 

0 

0.00 

0 

0 

0,00 

0.0 

1012 

0 

0 

0.00 

0 

63 

53,93 

0.0 

101? 

0 

63 

63.93 

0 

0 

0.00 

0,0 

1014 

30 

38 

6.69 

46 

38 

6.69 

0,0 

1016 

23 

39 

39.70 

38 

37 

18.74 

0.0 

1016 

36 

20 

20.30 

51 

22 

41.26 

0.0 


one to three hours as a function of the adopted clock error modeling 
procedure. 

To the geometrically derived observations of range and range 
difference, equations (3,1.1) and (3.1.2), systematic and random error 
sources were added as required In accordance with Table 4.2.2. White 
noise consistent with that expected from six-second ranges smoothed over 
300 seconds or one-minute Integrated Doppler range differences aggregated 
over the same interval was added based on the adopted levels in Table 
4.2.2 as described in Section 4 .1.4. Random rubidium and cesium clock 
noise were simulated using the algorithm of Med itch [1975], based on the 
selecte.d Allan variance models for the satellite and geodetic receiver 
oscillators. Random receiver cesium clock noise was added to the obser- 
vations along with a time bias and drift as given in Table 4.2.2. 
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TABLE 4.2,2. SIMULATION ERROR, SOURCES 


ORBITAL ELEMENT 

a 

c 

i 

M + CO 
l) 

TROPOSPHERIC REFRACTION 


PERIODIC ERRORS (RMS) 


RADIAL 

ALONG TRACK 

CROSS TRACK 

2m 



2.6m 

5.3 m 

12 m 


6.4 m * 


* PLUS Sm 

QUADRATIC ERROR 



5°. o OF HOPFIELD MODEL PREDICTION 


CLOCK ERRORS 
1. STATION 

BIAS: SO TO -10 nsec 

DRIFT: .000083 TO -.00004 nsec/sec 

RANDOM: CESIUM SPECIFICATIONS 


2. SATELLITE 

RMS BIAS: 5 nsec 
RMS DRIFT: ,0002 nsec/sec 
RANDOM: RESIDUAL NOISE BASED ON RUBIDIUM 
MODEL FOR GPS SATELLITES 

WHITE NOISE 

RANGE; 100 » 60 cm (6 sec) 

DOPPLER: 3 ■ 1 cm (60 sec ) 


104 


For each satellite the total time error converted to distance 1* 


p ( lven by equation (4,1,11) ao 


cD. « 

eTjCt) « ca^C^) + ci^Ct - t 0 ) +“ (t«t o r + n 1 (t) 

(4.2.1) 

where 


t 

»b(t) n c / y. (OdT . 

M . W 

t 

0 

(4.2.2) 

Assuming the ageing to be negligible over the satellite tracking Inter** 

val, equation (4.2.1) can be written as 


cT^(t) <=> a 4* b(t~ t*') +* F^(t) 

(4.2.3) 

where 

CD . n 

a « cT i (t Q ) 4* cR i (t"- t Q ) + -g±- (t- t Q )- + n 1 (f‘) 

(4.2.4) 

b « cR^ 

(4.2.5) 

^(t) « n i (t) - n i (t") . 

(4.2.6) 


Letting cT^Ct) be the best prior linear estimate of cT^t) over the 
satellite tracking interval, the residual range error Is defined as 

c [l*^( t) - T^t)] « (a- a) + (b - b) (t - t'*) + r^t) . (4,2,7) 

A 

The expected standard errors of the residual random bias (a-a) and drift 

A 

(b-b) of the satellite clocks were taken to be consistent with the 
ephemerls simulation results described In Section 4*1.2 and are given 
In 'fable 4.2.2. The residual random range error r^(t) Is obtained from 
q^t) by linear least squares approximation. These quantities are 
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added to the geometrically derived range and range differences accord- 
ingly, assuming the random bias and drift are normally distributed 
although constant within a particular tracking interval. 

A residual tropospheric refraction error In the form of five 
percent of the Hopflold model prediction was applied to the observa- 
tions In certain cases to represent a difference between actual and 
modeled refraction. In these cases the adjustment model (Section A. 3) 
Included tropospheric scaling parameters as indicated In equation 
(A, 1.37), 

The resulting quantities arc the observed range and range dif- 
ference observations given by 

R^(t) ° R(t) + V R /vn + a r b(t- O + n(t) 2 ^ 

+ (^-ap + < b ± ” b ± ) Ct “ K"*) + r^t) 1* 06R(t) 

and 

AR^(t) °R(t) - R(t- At) + /m V Ar + bAt + n(t) 

- n(t-At) + (b i -h i )ht + r^t) - r^(t - At) (A. 2. 9) 
+ 0[6R(t) - 6R(t- At)] 

where n is the number of six-second ranges assumed smoothed, V R are 
independent zero mean Gaussian random numbers having a standard devia- 
tion equal to that of the six-second ranges, 3 is 0.05, 511(f) is the 
tropospheric refraction error, m is the number of aggregated Doppler, 
observations, V A „ are independent zero mean Gaussian random numbers 
having a standard deviation equal to that of the Doppler measurements, 
t" and t A ‘ are clock model epochs, and At is five minutes. 
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The white noise level applied to range differences according to 
equation (4.2.9) corresponds to the summation of m independent one minute 
Doppler counts. For continuous count Integrated Doppler the total noise 
o£ the aggregated count may be less and possibly Independent of the 
count inf rval; thus the coefficient of In equation (4.2.9) may be 
as small as one. The consequence of this will be examined in Section 4.5. 

In the least squares adjustment •>£ range or Doppler observations 
actual measurements are differenced from estimates of the measurementr 
The estimated observation value is obtained from the geometrically 
derived range or range difference by a linear adjustment of this quan- 
tity which Introduces the assumed level of orbit error. For range 
observations this linear adjustment has the form 


6 

R (t) « R(t) + E 
o , 

k = 


1 


5R(t) 


Ae k (t c> 


(4.2.10) 


where the Ae^(t c ) represent errors In the orbital elements of the 
tracked satellite at the midpoint of the satellite tracking Interval. 
These errors are assumed to be nrxmally distributed varying with each 
satellite of the constellation. The expected ephemeris error is given 
in Table 4.2.2. The required partial derivatives in equation (4.2.10) 
are developed in Section 4.3.2 and approximated in Section 4.3.3. for 
range difference observations an equation analogous to equation (4.2.10) 
was utilized. 

Notice that the error in the satellite clock was introduced into 

A A 

the observed quantity along with the estimates a and b. Strictly 
speaking these error estimates should be introduced into the estimated 
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observation but the net effect on the difference of estimated and 
observed quantities in either case is identical. The same argument 
holds for the tropospheric refraction error. 

4.2.2 Opt imal Desip, n for Dynamic Point Positioning 

A Global Positioning System is designed so that six to nine 
satellites are usually available for observation from any geographic 
location. Since options in the tracking geometry are available, it is 
reasonable to design a data acquisition schedule which produces Live best 
results for the adjusted station coordinates derived from dynamic point 
positioning. Some factors to be considered in such a design are the 
length of tlve tracking interval for each observed satellite, possible 
criteria for minimizing the coordinate covariance, the period of site 
occupation and the type of receiver operation anticipated, sequential 
tracking or the use of multiple channel receivers, With these factors 
defined a sequence of satellites can be selected whose observations give 
the best geometrical strength of solution according to the criterion 
adopted. 

Various procedures for selecting the satellites to be tracked 
can be defined. These include the simplest approach of random selection 
from those visible to approaches based on choices for covariance mini- 
mization. For dynamic positioning performed using a sequential single 
satellite receiver two criteria will he discussed: 

(1) the selection of satellites whose observations minimize the square 
root of the trace of the accumulated covariance matrix, and 


tiO tlu 1 select ton of satellites whoso observations mint mi no the mum 
of the vector of correlation coefficients of the accumulated 
covariance matrix. 

The first procedure' is also known as the evaluation of the geometric 
dilution of precision. This second procedure allows a more statisti- 
cally independent determination of individual parameters. In evaluatinf 
these criteria observations are assumed subject only to Gaussian white 
noise. No correlated random errors are introduced. However in addi- 
tion to the station coordinates a receiver clock error model was incor- 
porated in some cases. 

Using global ly distributed stations a geometric analysis of 
point positioning for range and Doppler tracking was made examining the 
two criteria. The use of various station locations insured that 
numerous samplings of satellite pass geometry were utilized such as 
those shown In Figures 4, 2. 4 (a) and 4.2.4(b). For one particular sta- 
tion Figures 4.2.5(b) through 4.2.5(e) give the square root of the 
trace of the covariance matrix and standard error in latitude, longi- 
tude, arid height as a function of the numbe.r of one-hour satellite 
tracking intervals of range observations having a one meter standard 
error. Range observations were assumed every five minutes and no clock 
error model was included in this case.. Figure 4.2.5(a) gives the azi- 
muth and elevation angles for the epoch of each tracking interval for 
selection based on minimizing the trace of the covariance matrix. 

Figure 4 . 2 . 6 gives analogous information for the second criterion. Obser- 
vations below ten degrees elevation were excluded from the results* 
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Using those criteria the next satellite to be tracked la that 
whose observations over the upcoming Interval, when combined with nil 
previous observations, produces the optimum coordinate covariance with 
respect to the selection criterion, Notice in these examples that tlu* 
standard error of each position component drops rapidly within the first 
day then shows only gradual Improvement with additional data. An exami- 
nation of the results using the second criterion shows some reduction In 
the parameter correlations but yields an Increase In the expected 
standard error as evidenced In Figures A. 2. 6(b) through 4.2.6(e). 

Based on a number of similar determinations the following 
general conclusions can be drawn for the optimal selection of GPS satel- 
lites for both range and Doppler. First, it is readily apparent that 
the second technique results In somewhat lower parameter correlations 
but at the cost of Increased parameter variances with respect to the 
first criterion. However, the technique adopted Is a matter of choice 
since each is Independent. For the positioning studies of Section 4.5 
the first criterion was utilized to establish the observation schedules 
In all cases. Secondly, from the results It Is noted that initially the 
variance of the estimated parameters Increases rapidly as the interval 
of tracking each satellite Is increased. For a fixed number of observa- 
tions the results obtained are quite varied when the total observation 
time Is less than six hours. With Increased observation time allowing 
more sampling of pass geometries the results become virtually equivalent 
after twelve hours. Thirdly, with range observations the introduction 
of a receiver time bias significantly weakens the variance of station 
height as evidenced by comparing Figure 4.2.5(e) with Figure 4.2.7. 
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Figure 4*2.7. Standard Error In Height when Adjustment Includes 
Time Bias (Selection Criterion Minimum Trace) 

The variance of station latitude and longitude also Increase hut not as 
significantly. The reason for this increase in height uncertainty can be 
explained by noting that a time or range bias error In a given observa- 
tion Is equivalent to a linear combination of a vertical and horizontal 
position error of the station. The horizontal error lies along the 
projection of the slant range vector onto the horizon plane. Successive 
observations taken as a function of azimuth would yield horizontal error 
components whose sum would tend to cancel. However the vertical error 
component can only be separated from an actual station height error by 
using observations of low elevation. At ten degrees elevation the ver- 
tical component of range bias Is approximately seventeen percent of the 
total bias. Therefore with the restriction of observations to elevation 
angles greater than ten degrees a weakening of the actual station 
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height uncertainty can be expected since a station height error will 
bend to be mashed by the vertical component of range bias. 

Finally, the satellite orbits have repeating ground tracks 
yielding a tracking geometry with a diurnal period. For the first 
selection criterion this tends to result In a clustering of the initial 
satellite azimuth of each tracking Interval into a series of three Imnas 
separated by 50 to 150 degrees with a sampling of different elevations 
In each band. This is most obvious in cases whore the tracking interval 
is short. This property Is not a fixed rule but a general trend as 
demonstrated somewhat In Figure 4.2.5(a). For the second criterion the 
distribution of azimuth is less consistent, although In some cases eon- 
sldered the distribution may fall almost entirely within a single band 
of 150 degrees width. 

For the determination of baseline components and chord length 
only the first criterion was examined and the results were discussed in 
Chapter 3. 

4 . 3 Adjustment Procedure 

The adjustment of range and Doppler observations using the 
method of observation equations may be developed In a mathematical form 
which accommodates the introduction of new observations and new para- 
meters. IJotlla [1967] discusses this sequential approach and it is 
emphasized that such a technique is valuable In assessing the effect of 
additional observations on current parameter estimates. This approach 
may be adopted for the analysis of GPS observations taken In a sequen- 
tial fashion as discussed in Section 4.2. The estimation equations for 
this procedure are now developed. 
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where 


4.3,1 Sequential Adjustment of Parameters 

Given N statistically Independent sets of observations L. 

°1 


h a 

$ 

* 

f 

hi a rtx 'V 


(4,3.1) 


equations are developed for the least squares minimum variance estimate 
of the primary parameters x and secondary parameters y^ using all N 
observation sets. The primary parameters of Interest are the 
earth-fixed station coordinates or coordinate differences. The second- 
ary parameters consist of orbital elements, satellite clock model para- 
meters, refraction bias parameters, and tracking receiver clock 
parameters, Formulas giving the parameter covariance matrices T, and 

X 

X and the weighted sum square of residuals after adjustment, V PV, 
y i 

are developed, The sequential forms of these equations are given as 

A 

required. In sequential form the estimated quantity based on N + 1 
sets of observations is written as 


Z N+ 1 “ Z N + AZ N + 1 


(4.3.2) 


where Az^ + ^ Is the correction to the prior estimate z^ due to the inclu- 


sion of observations L, 


N+l 


4, 3. 1.1 Estimation of Primary Parameters 

The least squares minimum variance estimate for the parameters 

of primary interest x from any one set of observations K from equa- 

1 

tlons (4.3.1) is given by 
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i - -«»- v^w ; 1 [ v v£v* 


where 


» — £ R a v 
Si^i X 


N r 


N N 
xx xyi 


N 


y l X N 


y i y l. 


“ A l P l A i 


U. 


A » 
A i 


U 


U 


- A X L i 


y iJ 

k 1 


p = z 

L 


L. = L - L, 
i o i b 


L = F,(x ,y. ) 
o. i o J ± 
i o 


(4,3.3) 


(4.3.4) 


(4.3.5) 

(4.3.6) 


(4.3.7) 


(4.3.8) 

(4.3.9) 


2 

with the a priori variance of unit weight cTq equal to one. Equation 

(4.3.3) is the solution for x based on observations L. which results 

b ± 

after algebraic elimination of the secondary parameter set y^. 

With the addition of a second observation set L, the estimate 
for x becomes 



__ 

Z“ 1 + Z~ 1 

-i 

R. + R , 


L x i x jj 


i j _ 


(4.3.10) 


Denoting the covariance matrix E as E , equation (4.3.10) may be 

X 1 X I 

written using matrix identities as 


121 


(4.3.11) 


£ ti " -is. -K K +!; . hr.+kj 

J#JL Arjr A<J« Aj Aj *»tj» 4* J 

« x T - z a: + e„ )~ l x r - (e " 1 + Er.VV, 


X x X x x j 


X x x j 


or 


X I1 “ X I + Ax n • 


(4.3.12) 


The covariance matrix 5', , based on a pair of data sets , assumed uncor* 

X II 

related set-wise, becomes 


5: x “ 

II 


-1 -l’ 

£ + E 

• x i x j. 


(4.3.33) 


In general given x M and £ based on N .sets of observations, the 
N X N 

estimate based on the Inclusion of an additional observation set is 
given by 


*N 4 3. “ X N + Ax N 4 1 


(4.3.34) 


where. 


Ax N + 1 " [1 \ +1: x rlJ N 

” n N N n+3 1 

- [s^'+ir 1 f \ +1 

X N ii+l 


(4.3.15) 


and the covariance is 


-1 -1 -1 

£ » LE +E ] 

V v* v 

*N+1 *N n+1 


-a 

* Ar 


- £ -£[£+£ ] £ 

X N *N X N x n + 1 X N 


(4.3.16) 


or 


- £ + AE 

X N + 1 X N X N + 1 


(4.3,17) 
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where 


AE « -E 
X N + 1 *N 



+ E 

X n + 1 J 



(4.3.18) 


Equations (4.3.35) mid (4.3.18), although not necessarily computa- 
tionally efficient, give a measure of the expected change to the primary 
parameter set estimates and their uncertainties if a new observation set 
is added. It Is assumed here that: 

(1) new normal equations are formed using the same Initial parameter 
estimates x ,y as used previously; and 

0 X 

o 

(11) the new data block L. is uncorrelated with all previously used 

b N+l 

observations . 


4. 3. 1.2 Estimation of Secondary Parameters 

Consider the least squares normal equations for the observation 


set 


V 


N 

N 


X 


”u 

XX 

Xy j 



+ 

X 

N „ 

N 


y . 


U 

L v 

Vi] 


3 


>> 

l 


= 0 


or 


N x + N y , + U = 0 
x:; xy^ j x 

N x + N y , + U =0. 

V y j y / J y j 


The solution for y^ Is given by 


y, » -K~ X [U +N x] 

j y j y j y j V 


where x Is based on all N observation sets: 


(4.3.19) 


(4.3.20) 


(4.3.21) 
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x «=» -£ 


x 


V (u -N N U ) 

' V W V V V ' 


it-! WV‘ 

Substituting equation (4.3.22) into (4.3.21) yields 


(4.3.22) 


y . *= -N -1 A T P.L, +N“’ ! ' N Z 
d y j y j y j d d y j y j y j X X 


£(U - N N“\ r U ) 

' V VV V V V 


--if 1 , A 1 ' F.Ii. +N -1 ,, N Z 

Vj y j J 3 Vj V * 


AW w / J 

i * ^Wi 1 

-1 T 
N A 1 ' 

L i * *>'i Vi y i 


(4.3.23) 


Equation (4.3.23) can be written as 
-1 T -1 

t j -w t» •»* t it •*» 


y, = -N A P.L, + N N E 

d y ± y ± yj J J y/j y^ *] 


A T P,L ,-N n" 1 ,, A T P.L, 

X j j Xy d y d y d y j J d Jj 


+ Q l (4.2.24) 


where Q Is a function of all data sets except L . The covariance 


j 


matrix Z is obtained from equation (4.3.23). Since all data sets are 

y j 

assumed statistically independent the covariance Z Is given by 

y i 


N 

Vl" 

p 4 * 

" dy ,i~ 

II 

dL lc 

p k 

dL k 
— — 


(4.3.25) 


Differentiating equation (4.3.23) 

L IVt - H N _1 A T F ■ 

! X L * J *y d y/j yj 3 


-1 T -1 

, -N A x P.+N N 


dyj_ 

dL, 


U l\ 


y d y J y d 3 y :i y j y d x 


k = j 


k 


(4.3.26) 


,-l 


y J y i 'j 


n z [a t : 

, y.X X X 

T <r *1 1-1 


p, -h tr 1 a t p, 


k xy k y k y k y k k Jk 


1 


ki 6 j 


Substituting equation (4.3.26) into equation (4.3.25) and summing over 
k yields the final result 


Z = N _1 + N _1 N Z N N -1 


y, y,y, y,y, y,x * x y , y.y. 


(4.3.27) 



The residuals of flu after adjustment are Riven by the 


linearized form of equation (A. 3.1) 


A A 

.* y i. 


r n i 
x 

A 


+ 1., 


(4.3.28) 


or 


V i “ Vi + L 1 • 


(4.3.?‘0 


In terms of N data sets equation (A. 3. 28) has the form 


'V 


V 2 

C3 

-v 4 , . 



x 

V 2 

x 


.N 


A 0 
y l 

0 A 


.A" 0 0 . 

x 


. . . 01 

. . . 0 

. . A I 


A 

X -1 


rh 


+ 


Lk 


(4.3. 30) 


Since the observation blocks are uncorrelated 
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(4. 3. 3.1) 


The a posteriori variance of unit weight is 


a 


2 _ v^'pv 


(4.3.32) 


where d Is the number of degrees of freedom In the total adjustment 


problem. 
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The equations developed here were utilized In the adjustment of 
simulated GPS range and Doppler observations. It is an illustracive 
approach for determining how the uncertainty and error in station 
positioning vary In time as a function of such variables as the number 
of observations, method of satellite selection, tracking interval, and 
others . 


4.3.2 Model Parameters and Partial Derivatives 

In the adjustment of range and Doppler observations the number 
of secondary parameters y^ for each data set of equation (4.3.1) Is 
subject to variation depending on the tracking schedule and the choice 
of specific clock error models and additional bias parameters. The 
secondary parameter set Includes six orbital elements for each so tel- 
llte tracking Interval, a polynomial clock model for receiver clocks 
over the time span of each observation set, a polynomial clock model for 
each satellite clock over the Interval each is tracked, and may Include 
tropospheric refraction scaling parameters for every satellite-station 
combination within an observation set. The primary parameter set x con- 
sists of the geodetic coordinates of the tracking stations in the 
adopted earth-fixed frame of reference. 

The design matrix A introduced in Section 3.1.2 is developed 
from the first partials of the data function with respect to the model 
parameters. The partial derivatives of the range observation model, 
equation (4 .1.57) , with respect to the Cartesian earth-fixed coordinates 
of the tracking station are 
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u ■=>• v,w 


(4,3,33) 


Ou 


u } (t) - u 

"' to rr~ 


where* u ,v ,w are the coordinates of the satellite at time t in the 
ass 

same frame. The partial a for geodetic coordinates are obtained using 
equation (4.3.33) and the chain rule 


MSl w M jiL+ + M A + X h (4.3.34; 
0<|> 3u W dv H 3w v ’ 

where the partials of Cartesian coordinates with respect to geodetic 
coordinates are given in Rapp [1976] as the coefficients of the dif- 
ferential equations 


du = -(M+h)simf>eosAdf|) - (N+h)cos t[> slnAdA + cosificosAdh 

dv n -(H+h)sin(j)sinAd(|)+ (N+h) costjscosAdA + cosifislnAdh •^•SS) 

dw “ (M+h) c.ostf)dtf» + sintjidh 

where M and N are the ellipsoid radii of curvature in the meridian and 
prime vertical. 

The satellite coordinates are obtained in a mean inertial system 
by numerical Integration of the equations of motion whose forces Include 
the geopotential to degree and order eight, solar radiation pressure, 
and luni-solar gravitational perturbations and are rotated into the 
earth-fixed frame. The initial conditions for the Integration are 
obtained from Table 3.2.1. 

For range difference observations the partials for station 
coordinates are given by the difference 
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u ■+ v,w 


(4.3.36) 


8AR(t) „ u s (t> - u . u 8 (t - fll:) - u 
9u R(t) R(t-At) 


where At Is the Doppler Integration interval. For geodetic coordinates 
a similar expression holds 


MgM . Mt) . 8RV - MX * + A(h 

dtp d<p dtp 


(4.3.37) 


Each interval a satellite is tracked six orbital elements are 
Introduced to model ephemerls error. The elements are represented as 
the orbital semi-major axis a, the eccentricity multiplied by the cosine 
and sine of the argument of perigee, ecostoand e slnw , the inclination 
1, the sum of mean anomaly M and argument of perigee w, and the 
ascending node 0. Letting t represent the midpoint of the satellite 
tracking Interval and t the epoch associated with the initial orbital 
elements, the partial derivatives of range with respect to the elements 
at t are given by 


3R(t) b 3R(t) 

3e* 3X*(t) 

t* s 


tp (t)ijj" ;l (t )T(t ) 
r o o c c 


where 

3R(t) l ~3R(t) 3R(t) 3R(t) n n n 

'3X*(tj" 3X (t) 3Y (t) 3Z (t) 

S Us s s 

and 


(4.3.38) 


(4.3.39) 


Mt) X s (t) ~ X(t) 
3X (t) " R(t) 


X Y,Z 


(4.3.40) 


where X ,Y ,Z are the satellite Cartesian coordinates in the mean iner- 
s s s 

tial system. The quantities X,Y,Z are the station coordinates in the 
same frame. The matrix l|/ Q (t) is the state transition whose elements are 
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obtained from the homogeneous solution of the satellite's variational 
equations [baker, 1967]. 


and T is 



3X fl (t) 
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DZ,(t) 

3a(t o ) 

WTT 

the Jacobian matrix 
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3a(t c ) 
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C4 .3.41) 


(4.3.42) 


For range difference observations equation (4.3.38) is modified using 


where 


3AR(t) K f 3AR(t) 3 AR(t) 3AR(t) Q Q Q 
3X*(t) L 3X (t) 3Y (t) 3Z (t) 

5 Boo 


(4.3.43) 


3AR(t) V t? ~ X(t ) X s (t-At) -X(t - At) 

3X (t) R(t) " R(t-At) 

s 


(4.3.44) 


For polynomial clock models of the form 


P (t) = a + a. (t - O + ... -Ha (t - t**) n (4.3.45) 

n o i n 

the partials of range with respect to the model parameters a 0 >'..> a n 
are just the parameter coefficients 
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(4,3,46) 


3R(t) 

aa l 


(t- t") 1 


i ® o, , , , ,n 


where t' Is an arbitrary epoch for each model. The constant term In 
equation (4.3,46) represents the time (range) bias or phase error at t' 
associated with the modeled clock. The second parameter a^ is propor- 
tional tu the oscillator frequency bias Af^ 

a l c T ^ (4,3,47) 


and a ^ Is proportional to the frequency drift f 


- „ j£_ f 

a 2 2f * 


(4.3.48) 


For range difference observations derived from integrated Doppler over 
the interval [t^,t^] the partials are 


9|M0 

3a 

o 


3AR(t) = 
3a, 


t l~ fc k 


(4.3.49) 

(4.3.50) 


and 


(t £ -0 2 - (t k - tj ,) 2 . 


(4.3.51) 


In terms of frequency bias and drift the partials are 


3AR(t) B c , _ v 

3Af f ^ l V 


(4.3.52) 


and 


3AR(t) = _c_ [(t, - t") 2 

. 2 f * 

3f 

according to equation (4.1.7). 


(t k - t") 2 ] 


(4.3.53) 
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For the tropospheric refraction scaling parameter the partial^ 


of range and range difference are given by 


3 R(t) 

tBawi w't i wfo n O- 

^ C H 


6R(U) 


and 


" 6R(t k ) 

using equation (4,1,57). 


(4 , '3,54) 


(4.3.55) 


4.3,3 Use of Keplerlan Partial Derivatives 
In the Adjustment Model 

Numerical Integration of the variational equations [O'Toole, 

197ft] 

ji2L. es C r 7 r ,? P,).., r*)-r ( 4 3 5 ft) 

dr dr r o o C4.3.ob) 

O 0 

gives the variation in a satellite's Inertial position and velocity at 
time t with respect to changes in the initial state at time. t^. These 
partial derivatives are used In forming the observational partial deri- 
vatives of the design matrix A when the satellite state vector at t^ is 
included in the adjustment. For satellites at extreme altitudes an 
approximation may be introduced. This approximation consists of 
replacing the numerically Integrated solutions of the variational equa- 
tions with Keplerlan two body partial derivatives . This approximation 
is both economical and valid at GPS altitudes if observation times are 
within a few hours of the time at which the satellite state vector is to 
be improved. Keplerlan partlals are analytic expressions derived from 
the basic equations of two-body motion [Mueller, 1964]. 
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The orbital elements at time t are obtained from the elements 
at time tg by the transformation 


1 0 0 0 0 0 

0 1 0 0 0 0 

0 0 1 0 0 0 

a 0 0 1 0 0 

0 0 0 0 1 0 

0 0 0 0 0 1 


(4.3,37) 


where a Is n time dependent quantity derived from Kepler's third law 


a a v/~r (t- O • 


(4.3.56) 


For Keplerlan motion the only time varying element in equation (4,3.57) 
Is the mean anomaly 


M(t) « M(t Q ) + n(t - t Q ) 


(4,3.59) 


where n is the mean motion of the satellite. Differentiating equation 
(4.3.59) with respect to the orbital semi-major axis a at t^ gives the 
rate of change of mean anomaly given a change in a at t^: 


3M(t) 

3a(t 


on \ 

T“-27To a-t 0 ) . 


(4.3.60) 


All other element variations are of the form 


9 e i (t) 


i=j 


3e d Ct o ) (0 i^j . 


(4.3.61) 
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At observation time t let (u D (t) ,v fl (t) ,v Q (t) } be the coordinate*} 
of satellite position In the earth-fixed system. Let {u,v,w} be Che 
earth-fixed Carteolan coordinates of the tracking station. For range 
observations 


3R(t) 

9ujt) 


u (t)- u 

"teT"" 


JMjsl, a 0 
9u (t) 

c> 


u v,w 


(A. 3. 62) 


For orbital element estimation at time t” the partlals of R(t) with 


respect to orbital elements are given by the matrix equation 


r 3R(ti ’ 



3R(t) 

IfcFUT 

(6x1) 

3 Vl ti 

3u*(t) 

** ” m. 


(3 x l) 


(4,3.63) 


where the matrix G Is obtained from equations (A. 3. 60) and (A. 3. 61) 


G “ *- G lj^ “ 


9e , (t) 


3e i^". 


1 1 c j 

0 l^j except 1 »1, j « A 
^ (t - t") 1 ° 1 , J “ A 


(A.3.6A) 


The last factor in equation (A. 3. 63) has row dimension three since the 
velocity partlals of range are zero. The matrix H has the form 


3u (t) 
3a(t) 


II « 



3u g (t) 
_ 3f2(t) 


3w s (t) 

3a(t) 

3w (t) 

. s 

3«(t) J . 


(A. 3. 65) 
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The development of the elements of 11 follows. The transforma- 
tion from the mean inertial (X) system of epoch t Q to the earth-fixed 
(EF) system at time t is given by 




•a. *» 

X(t) 


“ R 2 (-x p )R 1 (~y p )R 3 (GAST)NP 

y(t) 

_w 0 (U>_ 

EF 



a R 


X(t) 

Y(t) 

Z(t) 


( 4 . 3 . 66 ) 


where the coordinates in the mean inertial system are given by 


X(t) 


x(t) 

Y(t) 

“ R 3 (-«)R 1 (-i)R 3 i.-w) 

y(t) 

Z(t) 


z(t) 


and {x,y,z} are the coordinates of the satellite in an in-plane coordi- 
nate system as defined in Mueller [1964] 


x(t) = a(cos E - e) 

y(t) ° a(l-e^)"^ slnE (4.3.68) 

z(t) a o . 

The quantity E is the eccentric anomaly related to the mean anomaly by 
Kepler's equation 


H = E - e sin E (4.3.69) 

and Xj, and y p are the coordinates of the Instantaneous rotation axis of 
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the earth [Mueller, 1969]. The rotation matrices P and N account for 
precession and nutation and GAST is the Greenwich apparent sidereal 
time. In terms of direction cosines equation (4.3.67) has the 
form 


’x(t)" 


x cos(xX) + y cos(yX) 

Y (t) 

J3 

x cos(xY) + y cos(yY) 

Z(t) 

I 

x cos(xZ) + y cos(yZ) 


Equation (4.3.66) may be differentiated with respect to the orbital 
elements 


3u (tj 
s 

3e(t) 


Using 

equation (4.3.70) 

and assuming a nearly circular orbit these par- 

tials 

derivatives are 



3u (t) u (t) 

s s 

(4.3.72) 


3a a 


3v s (t) v g (t) 

(4.3.73) 


3a a 


3w (t) w (t) 

s s 

(4.3.74) 


3a ~ a 


R 


ter 


3X(t) 


9e(t) 


(4.3.71) 


3u (t) 
s 

3e 


-a[R 11 cos(xX) + R 12 cos(xY) + R 12 cos(xZ)] (4.3.75) 

- a e sin - - B ' [R cos (yX) + R^cos (yY) + R^cos (yZ) ] 
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(4. 3.7c») 


3v (L) 

" “ a coa (xX) 4- R^cosCxY) 4 R^cosCxZ) ] 


n e slnK 



[R^cosCyX) 4* R 02 cos(yY) 4- R^eosCyZ) ] 


3w (fc) 

• - * *» -a [ R .j^coh (xX) 4* R^, , 0 . 0 s (xY) 4- R^cos (xZ) 


a 0 hIhH 




tR 3] C08(yX) 4* R.,.,c.os(yY) + R.,,,eos(yZ) ] 
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Du,(L) 

0M(T)' ,r; ” a sJn J UR u coa(xX) 4-R 12 cos(xY) 4* R^ cob (xZ) ] 

4- a Vl - c." coaRtR^jCOsCyX) 4- R r) coa(yY) 4- R^eosCyX) ] 


Ov^Ct) 

m'tj' 


- a aln K[R, n 

t* J.. 

4- a Vl~e * 


cos (xX) 4- R 2t) c,as(xY) 4- R^eos (xZ) 
coaK[R 0 . eos (yX) 4- R, )0 cosCyY) 4- R 


coh ( yX) ) 


aw s (L) 


“ “« sin r.[R >n cos (xX) + R, V ,C 0 S (xY) -f R„eos (xX) ] 
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33 


4* a V 1 - 0 “ cosE[R 31 coB(yX) 4- R, 0 coa(yY) 4- R v ,eos(yX) ] 


32 


33 


3u a C t ) 


|j~ a x[RnSin t sJ.ntosinS - Rj 0 s±n. i sinwcos&i 4- R^cos i slaw] 

+ ytRjjSin i cosuainfi- R 12 ain l cosucosfi + R^eos i cosw] 


V° 


x f R 2l sln 1 sin w sinS < “ R oo sin i sin (0 cos Si 4- n cob 1 sinw] 


+ yfRpjSin 1 c.aswsinll- R ?f) slnl cosiocosSl4> R,,„co« 1 cosw] 


('1 .3.77) 


(4.3.78) 


(4.3.7") 


(4.3.80) 


(4.3.81) 


(4.3.32) 


136 


?w (t) 
0 


xtK^^ain 1 sinwslnft - R^sin 1 slnWcopfiH- R^c.oh 1 aim)] 

+ y [R^sin i coswsinQ- U^sin i coswcosQ-f U^coa i com] 




Ou (t) 

' x C R n cos (yX) +R 12 cos(yY) + R 13 coa (yZ) ] 

- y[R u cos(xX) i- R^cosCxY) + R^^cosCxZ) ] 


(4 .3. 84) 


3v (t) 


•“ x[R 0 ^cos(yX) -I- R 22 cos(yY) + R^cosCyZ) ] 
y[R 21 cos(xX) + R 20 cos(xY) + R 23 cos(xZ) ] 


3 w q (L) 

“'Tfer~ = xlR 3] cos(yX) +R 32 cos(yY) 4 R^cosCyZ) ] 

- y [Ratios (xX) + R 30 cos(xY) + R 33 cos(xZ) ] 


(L) 

4j- r - x[-Rj cob (xY) i- R|2 cos 1 

+ y [~R 1; jCOs(yY) + R^cosCyX) ] 


(4.3.8.H) 


(4.3.8l>) 


(4.3.87) 


Civ (t) 

— = x[-R, J ^cos(xY) + R 22 cos(xX) ] 


+ y t- R 21 C0S (Y Y ) + R 22 cos 1 


(4.3. 88) 


9w (t) 

x[-R 31 c°s(xY) +R 32 cos(xX)] 

+ y t-R 31 cos(yY) 4 R q9 cos(yX)] 


32 


(4.3. 89) 


137 


Partial derivatives for Integrated Doppler measurements are 
obtained by differencing range observation partlals formed at the end 
times of the Doppler Integration interval. 

The difference in positioning introduced through this approxi- 
mation was determined from simulation results of absolute and relative 
positioning using range observations from two stations 100 kilometers 
apart. Results were obtained using two days of simulated observations 
with one hour of observation each time a satellite was acquired. 

Initial positioning adjustments were made in which orbital elements 
were included as parameters with a priori weights for each hour of 
tracking. Variational equations based on the WGS72 potential model 
[Seppelin, 1974] truncated to degree and order eight were numerically 
integrated and used in forming observational partial derivatives. Then 
the adjustments were repeated using the Keplerian two body partlals. 
Except for this modification the adjustments were identical. A com- 
parison of the covariance matrices obtained in each case was made. 

For absolute positioning the standard error of the station coordinates 
obtained using the Keplerian partial derivatives averaged 2.9 percent 
more optimistic. In the determination of coordinate differences the 
solutions using Keplerian partlals had standard errors averaging 2.4 
percent more optimistic. 

As a result of this experiment it was decided that Keplerian 
two body partial derivatives could be adequately used in the adjustment 
of station coordinates when orbital elements were taken as parameters 
in the adjustment. 
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4 . 4 Development of Adjustment Weight Matrices 


The minimum variance estimate of receiver earth-fixed coordi- 
nates obtained from range or Doppler observations by least squares 
adjustment requires that the weighting matrix be developed using the 
second order statistics of the random observation errors as outlined in 
[Llebelt, 1967]. In the application of Global Positioning System satel- 
lite observations of range and Doppler to geodetic positioning the 
adjustment weighting must include the observation error statistics for 
correlated atomic clock errors in both the satellite and receiver 
clocks and for noise from the tracking receiver. In this chapter the 
observation error statistics for atomic clock fractional frequency 
error are developed from the Allan variance for each system oscillator 
by an analytic procedure which transforms the Allan variance into the 
autocorrelation function for random frequency error. The integral of 
this function provides the statistics for range or range difference 
observations based on the two oscillators used to derive the measure- 
ment. Statistics of the residuals to selected polynomial clock models 
are obtained by an additional transformation of the range or Doppler 
error statistics. These residual statistics are incorporated with the 
instrumental white noise statistics into the adjustment weighting. The 
correlations between residuals to successive polynomial clock models are 
shown to be negligible allowing the adjustment to be performed sequen- 
tially. 
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4.4.1 Range and Doppler Observation Error Statistics 


4. 4. 1.1 Fractional Frequency Autocorrelation 
from the Allan Variance 

In section 4.1.1 the equations giving the second order statis- 
tics of random range and integrated Doppler observation errors due to 
random fractional frequency errors were presented. Those equations 
require that the fractional frequency autocorrelation function he known. 
In this section discussion of a procedure for obtaining an analytic 
approximation to this function from the Allan variance is given. This 
method avoids numerical difficulties that may arise when the inverse 
Fourier transform of the power spectral density is evaluated and yields 
simple analytic autocorrelation function. 

The Allen variance models shown in Figure 4.1.3 for the satel- 
lite rubidium and receiver cesium oscillators are a function of the 
sampling time t having the form 

N 

- T - T 1 
4 r t < x < t 2 

(4.4.1) 


3 

— T < T < co 
T 3 — — 

Using the transformations in Table 4.1.2 the power spectral density for 
fractional frequency may be developed from equation (4.4.1): 



s 


yy 
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0 .< u) <. to 
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w < w < to 

O'* — , 


W X .1 u £ w 2 

w 2 ( ' J — w h 


(4 ,4.2) 


The square root of the power spectral density, or transfer function , 
corresponding to the Allan variance specifications of Figure 4 -.1.3 Is 
given in Figure 4,4.1. The constants associated with the two functions 
and formulas for computing the constants associated with the power 
spectral density function based on the Allan variance are given In 
Table 4.4.1. These formulas are developed from the transformations of 
Table 4.1.2. 

The autocorrelation function $ (t) can be obtained from the 

yy 

power spectral density using equation (4.1.21) 
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With S 
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(co) an even function equation (4.4.3) reduces 
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(4.4.4) 


Using the power spectral density model, equation (4.4.2), in 
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Figure 4*4*1= Oscillator Transfer Functions 


TABLE 4.4.1. OSCILLATOR PARAMETERS 


QUANTITY 

UNITS 

FORMULA 

SATELLITE CLOCKS 
IRUBIDIUMI 

STATION CLOCKS 
(CESIUM) 

n 

sue 


1 D0*10 3 

1 00*10° 

>2 

DOC 


1 00*10° 

1 00*10° 


SOC 


1 00*10° 

100*10 7 

“'0 

aoe ' 

VTlij 

1 73*10 6 

173*10 7 

'"1 

SOC 1 

0 In2l(ntjl 

132*10 6 

1 32*10 8 

u'2 

SOC 1 

n«Zt 1 li)2l 

227*10 3 

2 27*10 8 

»i 



6 00*10 13 

3,00*10 14 

No 

sac 

’T’l 3 

3 60*10 77 

3 00*10 73 

"i 


mtf/21112 

8 16*10 SB 

2 04*10 71 

N, 

aoc ’ 

3,< t 7 h 2 

1 08*10 77 

2,70*10 33 

No 

sue 

°(V ,r 2 

3 60*10 

9,00*10 70 

it 


Ittljfci'j)'*® 

236*10° 

1 61*10° 

W« 

SOC ’ 


2 03*10 6 

1 67*10 6 


equation (4,4,4) gives the autocorrelation function for fractional 


frequency as 
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However this form for the autocorrelation function has an oscillatory 
behavior for small t as shown in Figure 4.4.2 as a result of trans- 
forming the band limited white noise portion of the spectrum. This is 
an artificiality of the model. 
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Figure 4*4-. 2, Satellite Rubidium Standard Fractional 
Frequency Autocorrelation from Inverse 
Fourier Transform 


An alternate approach for obtaining an autocorrelation function 
is to approximate the power spectral density model with a smooth func- 
tion whose autocorrelation is expressible in simple analytic form. The 
first step in this development is to approximate the flicker noise seg- 
ment of the spectrum with a series of cascading functions whose value 
alternates between being constant and being Inversely proportional to 
the square of the frequency. This type of procedure is described by 
Meditch [1975] in constructing a linear system which simulates flicker 
noise using a white noise input. Figure 4.4.3 shows the transfer func- 
tion for flicker noise. A three stage cascading transfer function is 
superimposed consisting of the functions F^, F^, and F^ which are 
defined in Table 4.4.2. These functions are defined to have the 
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required properties and give a continuous although not smooth approxi- 
mation to the flicker noise power spectral density. 

The constants of this approximation are now derived over fre- 
quency Intervals ns given In Medltch [1975]. The general form of the 
function F^ Is 

N a 

F (w) “ Hr (4,4.6) 

U) 


between the frequencies w and aw . At frequency w » defined in 

u a ti 

Table 4.4.1, the function takes on the value 


N, 


F> ) 

An 2 W, 
w i 


(4.4.7) 


since the flicker noise power spectral density has the same function 
value at frequency w^. Solving equation (4.4.7) gives 


Sw 2 

c aw.N., 
w, 11 


(4.4.8) 


A similar analysis gives the constant N^. The function F^ has the form 


n b 

F b (u) « 2 . 

w 


(4.4.9) 


At frequency a W^, Fg has the function value 


9 N a n a 

V“ - T®I ■ JJ 

a a 


(4.4.10) 
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since at a*u) the function !?„ has the same value no function 1\ at fre~ 
a B A 

quency aw fl (ace Figure 4.4. 3), Solving equation (4. 4 JO) given 

Ny R a (4.4J1) 


uolng equation (4. 4. 8). For the function F^, 


F c Cw) 



(4.4.32) 


Its function value at frequency <A») equals the value of at fr«- 

3 

quency a w, giving 


/ *»_ N n 

F (a V) ° 

C “ a°a > 2 «V 

a a 


(4.4.13) 


resulting in the solution 

N c a a^N B c ot^w^N^ (4.4.14) 

using equation (4.4.11). Numerical values for a and W a are given In 
Table 4.4,1. The power spectral density consisting of the three eas~ 
coding functions and the remainder of the original function will be 
denoted as the second power spectral density model for each oscillator, 
The next step in the development of an analytic autocorrelation 
function is to approximate various segments of this second model with a 
first order Markov process power spectral density function, a function 
of the form 
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(4.4,15) 


S(w) 


Cl 


2(J 2 B 

w 2 + 0 2 


where B Is the Inverse of the correlation time (see [Gelb, 1974]). 

The autocorrelation function for a first order Markov process la 

t(t) » (4.4.1b) 

Graphs of these functions are given In Figure 4,4,4, Notice In equa- 
tion (4.4,15) that the power spectral density decreases as the Inverse 
of the square of the frequency which Is the type of functional behavior 
seen In the Interior of the cascading functions F^ through F^. It is 
also the behavior of the original power spectral density in the inter- 
val [u)q,Wj], In addition the power spectral density of the Markov pro- 
cess remains virtually flat until the frequency readies a point when the 
function decreases rapidly. These properties make this function an 
excellent choice for approximating the second power spectral density 
model piecewise. 

The second model Is then divided into five segments defined In 

-1 

Table 4.4.3. The high frequency cut off ey shown as 10 In Figure 
4.4,1, will be Increased so that that band limited white noise component 
of the power spectral density may be approximated better by the first 
order Markov power spectral density. 

The approximation consists then of fitting a function in the 
form of equation (4.4.15) to each subdivision of the econd model 
Sy (w) given in Table 4.4.3. There are two parameters a and B to be 
determined for each segment giving a total of ten parameters. 
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AUTOCORRELATION FUNCTION 


POWER SPECTRAL DENSITY 



Figure 4.4.4. Markov Process (First Order) 


Two procedures for this approximation were examined. The first 
was a least squares fit of the function S(w) to each segment of S'^(w) . 
The second, which was adopted for use, was an assymptotic approximation 
whereby two constraints were imposed on the Markov power spectral den- 
sity function giving a and |3 directly. The second procedure was impli- 
mented because of simplicity and because the results compared favorably 
with the least squares approach as seen by comparing Figures 4.4.6 and 
4.4.7. The assymptotic approach develops an approximation on the inter 1 


val I. , 
J 


ij ~ [vV 


(4.4.17) 


using the following constraints: 
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(i) at zero frequency the approximating Markov power spectral density 
equals the second model at frequency 

S j (0) “ 2 8) 

(11) In the limit as w increases the value of the function S. (w) con~ 

.1 

verges to the following function 

2 

lim S. (w) o - 2 £Ji (/» .4 . 1 *3) 

• to" 

w ->• « 

and at o , this limiting value is set equal to the value of S'* (w) : 

t 1 yy 

2 

. (4. A, 20) 

Equations (4,4.18) and (4.4.20) are a system of two equations in 
two unknowns . Their solution yields the parameters a and p for the 
approximating Markov power spectral density function S. (w) . The nature 
of the second constraint, aquation (4.4.20) is to force the function 
Sj (to) to assymptotically approach S'^(w) at The first constraint is 

necessary to approximate the white noise or flat component of S* (w) at 

yy 

the beginning of each subinterval. 

Finally a comment concerning the approximation in the last sub- 
division I,. is necessary. In order to obtain a good approximation to 
Syy(w) in that interval it is necessary to choose large enough to 
allow the flat portion of the Markov process spectral density to fit the 
white noise component which dominates this interval as seen in 
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SATELLITE RUBIDIUM ACF 


TABLE 4.4.3. DIVISION OF SECOND PSD MODEL 
FOR MARKOV PROCESS 
APPROXIMATION 

NOTATION INTERVAL 

^ [0, «,] 
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Figure 4.4.5. A.ssymptotic Fractional Frequency Autocorrelation 
Functions Based on Markov Process Approximations 


Figure 4.4.1. Choosing uij three or lour orders of magnitude larger than 

0.1 and S'* (u), ) two to three orders of magnitude smaller than N n enables 
yy h 0 

a good approximation to be made but adds power at these higher fro- 
quencies. The. result is an autocorrelation function whose variance will 
increase as ix)^ is chosen larger (see Figure 4.4.5) and tends to a delta 
function as goes to infinity. However, this will have negligible 
effect on range and range difference statistics. This point will be 
examined in more detail after the development of additional equations 
based on the first order Markov approximations. 

The smooth fractional frequency autocorrelation function (t) 
is given by the inverse Fourier transform of tbe five Markov process 
power spectral densities Sj (w) . The result of each transformation is an 
analytic function whose form is given by equation (4.4.16) • The final 
result is the sum of these, functions 

* (t) » t • ' CA.4.21) 

yy jfi 

For range ami integrated Doppler observations tbe statistical contribu- 
tion duo to random oscillator error is obtained using equation (4.4.21) 
in equation (4.1.19) through (4.1.26). Figure 4.4.9 illustrates tbe 
steps discussed in tbe. development of these statistics from the Allan 
variance model. 

Figures 4.4.6 and 4.4.7 show tbe original transfer function for 
the satellite rubidium oscillator with the smooth least squares and 
assymptotlc approximations. The least squares fit to each subinterval 
of the second model S'* (w) was based on two hundred equally spaced 
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Figure 4.4.6. Satellite Oscillator Transfer 
Function and Sum of Least 
Squares Approximations 
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Figure 4.4.7. Satellite Oscillator Transfer 
Function and Sum of 
Assymptotic Approximations 
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Figure 4. 4. 9. Development of Uanpe and 
Doppler Statistics 
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TABLE 4.4,4, FRACTIONAL FREQUENCY AUTOCORRELATION 
FUNCTION PARAMETERS FOR MARKOV 
PROCESS APPROXIMATIONS 


ASSYMPTOTIC LEAST SQUARES 


OSCILLATOR TYPE 

INTERVAL 

(ALPHAS 

BETA 

(ALPHA) 2 

BETA 

RUBIDIUM (SPEC) 

b 

3.1177*10 24 

1732*10 * 

3.3719*10 24 

1,601*10 9 


h 

0,2020*10 28 

2 032*10 8 

7.0238*10 28 

2,250*10 8 


h 

0,2626*10 28 

1.120*10 4 

7,6240*10 28 

1.262*10 4 


<4 

0.2020*10 26 

0,282*10 4 

7.6245*10 28 

0.947*10 4 



1.0000*10 19 

1.000*10 9 ’ 

1,9343*10 19 

9.631X10 2 ' 

CESIUM (SPEC) 

It 

7.7942*10 27 

1,732*10 7 




u 

1.2922*10 27 

1.677*10 8 




h 

1.2922*10 27 

4.321x10 6 




U 

1,2922*10 27 

1.113*10 6 




h 

4.0000*10 20 

1. 000*1 o 8 * 




■tof„ <= i.ijxiO 4 S'yy|w h | a N o l1.0*10 2 

samples of the function within the subinterval. The parameters 
obtained using each approximation procedure arc given in Table 4.4.4 
for this rubidium oscillator. Since the assymptotic procedure produced 
results comparing favorably with the least squares procedure this method 
was adopted for use. Hence no least squares parameters appear in 
Table 4.4.4 for the cesium oscillator. The assymptotic transfer func- 
tion for the cesium oscillator fractional frequency error and the origi- 
nal power spectra • ensity are shown in Figure 4.4.8. 

4. 4. 1.2 Observation Error Statistics 

Based on Markov Process Approximations 

The first order Markov autocorrelation function, equation 
(4.4.21) and equations (4.1.19) through (4.1.26) give the second order 
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statistics for random range and Integrated Doppler observation errors due 
to each oscillator used In the measurement process. These Integrals may 
be evaluated giving analytical expressions for the variance and covari- 
ance of range and Doppler observations. 


A. 4. 1.2.1 Range Observation Statistics . Let R(t^) and R-(t^) 
be range observations subject to random clock error only. The covari- 
ance between the observations Is given by equation (4.1.19). Using the 
first order Harkov q\ / > ‘ sctlons, the Integration of equation (4.1.19) 
gives the covariance ae 


E[R(t 1 )R(t k )] » E[n(t i )n(t k )] 


, 5 0 r 

, 2 V _JL 

3, 


w-P,(t-t) (4.4.22) 

J 3 


j - 1 M j L 

~3,(t. -t) 


+ e 


y it s' 


'j vw k w i' 


-0 


for t, greater than t, where t„ is the start or reset time of the 
lc is 

clock. The variance of the random range error Is obtained by setting t‘ k 
equal to t^ in equation (4.4,22) 


E[R(t ± )R(t i )] = E[n(t i )n(t i )] 

2 


(4.4.23) 


= c 


i%[ 

3=1 P j *- 


i -3. (t, - t ) 

3 - 1 ) 


i 


The range error T|(t) resulting from the integration of frac- 
tional frequency error y(t) is a statistically non-stationary process. 
An examination of equation (4.4.22) and (4.4.23) reveals terms in these 

expressions which are functions of t. or t. minus t . Thus the 

i Lc s 
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variance, for Instance, increases with time, This is illustrated in 

Figure 4,4.10 for the rubidium clock. The standard error of a range 

measurement based on the use of this clock is given for twenty range 

observations spaced at fifteen minute Intervals starting five minutes, 

one hour, and five hours after the start of the clock. The increase in 

variance is almost linear. An examination of the autocorrelation func™ 

tlon shows that this function, domlnately flat, is similar to a random 

bias having a constant autocorrelation and whose integral is a random 

ramp which increases exactly linearly. Hence a linear growth in 

variance is expected as seen in Figure 4.4.10. The correlation coeffl- 

dents p u between the first range observation and the i'th in 

each of these sequences are given in Figure 4.4.11. As the starting time 

of the sequence from t Increases so does the correlation among the ran- 

s 

dom errors which again is expected since the variance increases with 
time and the errors are correlated. 

Figure 4.4.12 gives the autocorrelation function for the cesium 
clock based on the Markov process approximation and Figures 4.4.13 and 
4.4.14 give the standard error and correlations of range errors based on 
this clock. A comparison of Figures 4.4.10 and 4.4.13 reveals the 
greater stability of the cesium clock. After ten hours of operation the 
standard error of the cesium clock output is approximately 3.5 nano- 
seconds compared to 63 nanoseconds for the rubidium standard. In addi- 
tion the correlations among the cesium clock errors decreases more 
rapidly than the rubidium clock. 

Considering both random clock error sources the total variance 
and correlation of range observations y tp and y tj ) measured by 
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receiver k are given by the equations 


■'•'WW 1 “ E WW 1 + ‘WVV 1 


where the variances and correlations of the random error n are given by 
equations (4,4.22) and (4.4.23). The subscript "a" refers to the satel- 
lite rubidium clock. 

For simultaneous observations of range bv two receivers the 
covariance of the observations R k (t^) an ^ R^(tj) is given by 


E[R lc (t i ) \ (t j )] " * (4.4.26) 

In the above equations the random errors n have zero-mean which is a con- 
sequence of fractional frequency error being zero mean. 


4. 4. 1.2. 2 Integrated Doppler Observation Statistics , Let 
AR(t n ) be an integrated Doppler or range difference measurement over 
the Interval [t.,t ] and AR(t c ) a similar measurement from the same 
receiver over the interval [t^t^]. The covariance of the observations 
Is 


E[AR(t n ),AR(t £ )] “E[n(t a ) -n(t i ),n(t^) -n(tj c )] 


»E[n(t n ) n(t^)] - E[n(t n ) n(t {< ,)3 - E[n(t i ) n(t^)] 
+ E[n( t ± ) n(t fc )] 


2 A fi 


E 

i = i ej L 


• -3. (t 0 - t ) -3. (t, - t ) 

e i ^ n _ e j k n 


(4.4.27) 


-3 (t £ - t i ) t ± ) 
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Figure 4. A. 12. Assymptotic Fractional Frequency Autocorrelation 
Function for Cesium Standard 
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Figure 4.4.13. Standard Error of Range Observations 
Based on Cesium Oscillator 



The variance of n range difference observation la given by 
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E[AR(C n )AU(t n )] «=> c 
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(4.4,20) 


Equations (4,4,27) and (4,4.28) are Independent of the clock epoch t Q . 
The statistics of the range difference error depend only on the Doppler 
Integration Interval or the time difference between observations. Thus 
the random range difference error Is stationary. Expressions analogous 
to equations (4.4.24) through (4,4,26) express the complete statistics 
of range difference observation errors for individual or simultaneous 
observations due to clock error. 


4.4.2 Statistics of Residuals to Polynomial Clock Models 
The statistical characteristics of fractional frequency error 
and Its Integrated effect on range and Doppler observations have been 
discussed in detail. For range observations the total random error is 
due to three sources, two of which are correlated noise processes. The 
total random range error is expressible as 

n(f) a n s (t) + n k (t) + KM (4,4.29) 

where n and rt are the correlated random range errors due to satellite 

G K 

and receiver random clock errors respectively. The quantity £ repre- 
sents receiver white noise as discussed in Section 4,1.4. The total 
integrated Doppler random error over the Integration Interval [t^t^] Is 

An( V " W • n a <t j ) + W ' \ Ct J> + h 
where is the white noise associated with the Doppler measurement pro- 


cedure. 
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Figure 4,4,14. Correlation Coefficients between 
Range 1 and Range 1 (Cesium Clock) 

Depending on the stability of the clock the random range or 
Doppler error components, ri (t) and n (t), may appear quite systematic 
over fixed time Intervals and may be represented by polynomial models of 
varying degree, For short time spans the models for the clocks con- 
sidered In this analysis were taken to be a bias and drift for range 
observations or a frequency bias for Doppler observations. However 
these models and even higher order polynomial models are not sufficient 
to entirely represent this correlated error. Thus knowledge of the 
statistical properties of the deviations of the error from such a model 
becomes Important since these residuals represent an unmodeled part of 
the observation equation after the Inclusion of the polynomial model. 
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Proceeding, equation (4.4.29) Is expressed as follows 


ri(t) « P ms (t) + P^t) + r 8 ( fc ) + ^(t) + (4.4.31) 

where P (t) Is aa m’th degree polynomial chosen to model the corre- 
ms 

lated random error n (t) and P , (t) is an n’th degree polynomial 

s me 

modeling the random process • The statistics of the residuals 

r(t) may be developed from the covariance of the random clock errors 
developed in Section 4.4.1 using the procedure derived in Appendix A 
which develops the mathematics for polynomial approximation to random 
walk segments. Using equation (A. 1.9) the second order statistics of 
the range residuals r(t) to a polynomial model are obtained as 

E[r(t)r T (t) ] » GE[R(t)R T (t)]G T (4.4.32) 

where 

G = [I - A(A T A)“ 1 A T ] (4.4.33) 

and A is the design matrix for the polynomial model selected. The 
E[R(t)R (t)] is the covariance matrix of the random clock error being 
modeled. This covariance is given by equations (4.4.22) and (4.4.23). 

For integrated Doppler observations the statistics of the resi- 
duals to a given degree polynomial model are similarly obtained from 
equations (4.4.32) and (4.4.33) with the use of the covariance matrix 
for integrated Doppler random error due to each system clock, equations 
(4.4.27) and (4.4.28). The equation may be written as 

E[Ar(t)Ar T (t)] = HE[AR(t)AR T (t)]H T (4.4.34) 

where the matrix H is similar to the matrix G of equation (4.4.33) 


with changes due to the choice of the model adopted for clock induced 
random Doppler errors 

T T 

11 « Cl - A"(aV)" 1 A''] . (4.4.35) 

Equation (4.4.30) has the form 

AtlCt^) «=» ^i B (t^) + + Ar s^V + ^k^J^ ' H S (4.4.36) 

after the selection of the polynomial models. 

If the statistics of these residuals were ignored in dynamic 
point positioning adjustments the resulting coordinate covariance 
matrix would be optimistic. An Increase in the degrees of the poly- 
nomial clock models would offset this optimism to some extent since the 
level of unmodeled error would be decreased. However if a rigorous 
adjustment is to be performed then these residual statistics must be 
included in the least squares adjustment weight matrix to account for 
the unmodeled error r(t) or Ar(t) in a statistical rather than para- 
metric fashion. The adjustment should then produce a valid coordinate 
covariance matrix regardless of the order of the polynomial models used 
provided numerical problems are not encountered and the parameters are 
independent and well observed. 

The question of adequacy of a particular polynomial for a given 
data span needs to be addressed. 

4. 4. 2.1 Comments on the Choice of 
Polynomial Error Models 

To determine what degree polynomial model would be best to 
represent random clock error various factors have to be considered. 
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First among these are the length of the data span being modeled and the 
number of additional parameters which need to be introduced into the 
geodetic positioning adjustment. Use of a higher order polynomial will 
reduce the variance of the clock model residuals but may tend to numeri- 
cally weaken the adjustment normal equations. 

To determine how well a given order polynomial model repre- 
sents correlated clock error over a fixed time Interval a series of 
first and second order polynomial fits were made using simulated ran- 
dom clock error. The algorithm of Mediteh [1975] was used to generate 
sequences of clock error which were then converted to range error. The 
polynomial fits were equal weighted least squares approximations to the 
range errors. A sampling rate of one minute was used. From the 
residuals of fit r(t) autocorrelation functions were numerically 
obtained for each approximation using 

$,(t) “ “ 2 r(t )r(t, + r) (4.4.37) 

3 n i=l 1 1 

where n depends on t and the total number of samples. The variance of 
the residuals from each case were avei-aged to determine an overall 
variance for the residuals of fit for both the linear and quadratic 
polynomials . For the rubidium clock three time, intervals were 
considered with a linear polynomial fit. The root mean square 
errors are given in Table 4.4.5. For the cesium clock the 
results indicate that the longer the interval the better the second 
order polynomial performs, as expected. Howeve' this Increase in good- 
ness of fit is less significant as the length of the interval decreases. 
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TABLE 4. 4. 5. RANGE RESIDUAL STANDARD ERRORS BASED ON 
POLYNOMIAL FITS TO SIMULATED CLOCK ERROR 


CLOCK TYPE 

POLYNOMIAL MODEL 

NO. CASES 

RMS ERROR OF FIT (cm) 

4 hr 

ESS 

16 hr 

24 hr 

CESIUM 

LINEAR 

10 

8.7 

11.8 

18.7 

25.6 


QUADRATIC 

10 



9.3 

16.1 

17.4 


2 hr 

4 hr 

8 hr 

RUBIDIUM 

LINEAR 

30 

18.1 

36.5 

66.3 


It is obvious that a tradeoff exists between the level, of model error 
remaining and the number of model parameter required. For instance two 
linear models over 1G hours leaves an 11.8-centimeter sample standard 
error for the residuals, while a single quadratic fit over the same 
interval leaves 15.1 centimeters of expected error. An increase in one 
parameter produces a 22-percent decrease in the expected error. For 

the rubidium clock the expected level of residual error is higher due to 
the poorer short term stability of this clock (see Figure A. 1.3) and is 
comparable only to the cesium if the fit Interval is about one-eighth 
the length. Figures A. 4. 15 and 4.4,16 give examples of the residuals of 
fit for each clock for a linear fit over 8 hours. In Table 4.4.5 the 
length of the rubidium clock cases was limited to 8 hours since one 
clock model for each satellite pass was anticipated for the positioning 
studies to be conducted. 

Finally, the theoretical standard errors for range residuals to 
a linear fit were determined using equation (4.4.32) for these clocks 
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for the same intervals with the exception of the 4-hour case for the 
cusium clock. The results are Riven in Figures 4.4.17 and 4.4.13. 

These figures support the conclusions drawn above, and in addition, 
graphically demonstrate that the statistics of the residuals to the 
clock modeling polynomial are not stationary. The variance of a 
residual depends on the order of the polynomial, the interval length 

and the location within Che sample. The correlation coefficient matrix, 

contoured in Figure 4.4.19 for an 8-hour linear fit for the cesium 
clock, does not have the constant diagonals except for equally spaced 
samples of a stationary statistical process. However, by Theorem A,1 
of Appendix A, the statistics of the residuals will be constant from 
Interval to interval of the same length provided the sampling is per- 
formed equivalently and the same order polynomial is used. 

4. 4. 2. 2 Correlation Between Sets of Residuals 

An examination of equations (4.4.22), (4.4.23), (4.4.27) and 
(4.4.28) shows that the random errors due to oscillator instability are 
correlated over all time. That this is the case is a consequence of 
the error being a random walk or the difference in elements of a random 
walk where the underlying process is fractional frequency error. 

Since correlation between range or Doppler observations is due 
entirely to clock error, it becomes interesting to examine the cor- 
relation between the residuals of two successive polynomial fits to ran- 
dom clock error. If the cross correlations are relatively small, the 
assumption that successive ^servation sets can be taken as statisti- 
cally independent is justified when polynomial clock models are adopted. 
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Figure 4.4.17. Standard Error of Satellite Rubidium 
Clock Residuals Bapud on a Linear Fit 


a 



NUMBER Of SAMPLES 


Figure 4.4.18. Standard Error of Cesium Clock 
Residuals Based on a Linear Fit 



This as sump cion would then permit sets of correlated observations to be 
Introduced Into a geodetic positioning adjustment as Independent blocks 
in a sequential least squares approach. Computationally this Implies 
that the dimension of the observation covariance matrix to be Inverted 
to form the least squares weight matrix Is reasonable. 

To test that assumption the residual covariance matrix was 
computed using equation (A, 2. 5) of Appendix A for two linear fits 
to successive clock error segments for both the cesium and rubidium 
clocks : 


where 


and 


nmm •— TT 1 Bccc/P 

ECrjrJ] E[r r -*] 
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(4.4.39) 


(4.4.40) 


The correlation coefficient matrix for the residuals was computed and 

<_~fP 

the coefficients from the off-diagonal block, E[r^r 2 J were compared to 
the correlation coefficients from the diagonal blocks, E^rJ] and 
*^ r 2 r 2^ * results indicated In all cases that the correlation coef- 

ficients between residuals from two different fits were at least two 
orders of magnitude smaller than the correlation coefficients for 
residuals from the same polynomial fit. These results support the 
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assumption that successive blocks could be treated as independent 
although each block itself would be Internally correlated. 

4.4.3 Weight Matrix 

The Introduction of polynomial clock models tends to statisti- 
cally decouple the residuals from successive polynomial fits; thus., GTS 
range and Doppler observation sets, which are correlated in time by 
random clock error, may be treated as Independent when polynomials are 
adopted to model these random components. Each set is Itself fully 
correlated and the statistics of the residuals to the adopted poly- 
nomial models must be Included in forming the least squares adjustment 
weight matrix as shown In Figure 4.4.20. The size of each correlated 
data set will depend on the time interval over which the models are 
applied which, along with the degree of polynomial, determines the 
variance of the remaining residuals. Since the receiver cesium clock 
has better stability than the satellite rubidium clock the time inter- 
val over which a single receiver polynomial clock model is adopted may 
span multiple intervals of satellite tracking data each with its own 
clock model. This will of course depend not only on clock stability, 
but also on the geometric strength of the observations taken. Figure 
4.4.21 illustrates this concept in which observations within block K 
are assumed statistically Independent of observations within block L, 
each of which Includes range observations from four satellites taken :Ln 
this case simultaneously from two stations. 

The weight matrix, taken as the inverse of the covariance matrix 
of random observation errors, is assumed to be block diagonal wherein 
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each block contalno the occond order statlotlen of the residuals to the 
selected polynomial clock models and of the white noise due to the 
receiver. The weight matrix corresponding to the observation schedule 
of Figure 4.4,21 in shown in Figure 4.4.22, Each diagonal block is the 
inverse of the sum of three covariance matrices. For relative posi- 
tioning using range, observations from two stations observing simultane- 
ous the form of the diagonal block is given by equation (4.4.41) where 
the covariance matrices K[rr ] are based on equation (4.4.32). The form 
of the matrix is Identical for Doppler observations using the covariance 
KfArAr ] given in equation (4.4.34). 

4 , 5 Results of Dynamic Positioning Studies 

The simulated range and integrated Doppler observations 
developed in Section 4.2.1 according to equations (4.2.8) through 
(4.2.10) were used in the sequential least squares adjustment algorithm 
developed In Section 4.3.1 to obtain minimum variance estimates of 
geodetic station coordinates and baseline components using a dynamic 
positioning approach. Observations from three separate station 
groupings wore considered in this analysis. The geodetic coordinates 
of these stations are found in Table 4.2.1. The GPS orbxtal elements 
adopted xn this study are given in Table 3.2.1 referred to the mean 
equator and equinox of 1950.0. 

Solutions were obtained for the geodetic coordinates of indi- 
vidual tracking stations then for baseline components and chord lengths 
from simultaneous observations from pairs of stations. Solutions were 
developed using either range or integrated Doppler observations sepa- 
rately. No solutions based on both observation types were considered 
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alth&ugh this possibility way be available with two frequency Doppler 
receivers having a two frequency ranging capability. For each posi- 
tioning problem the effects of the random and systematic error sources 
ot *> ■! * - tv 4.1 were addressed and the adjustment weighting procedure 
developed In Section 4,4.3 was utilized as a function of the random 
error sources considered and the error models chosen to represent atomic 
clock error. 

Ii.tegrated Doppler observations were assumed to be Independent 
sixty-second measurements aggregated every five minutes, not correlated 
range differences, as In Chapter 3. This latter type of treatment would 
add additional strength to the least squares normal equations enhancing 
the Doppler results presented below. The type of correlations con- 
sidered In this analysis however are those due to the correlated random 
atomic clock error present in both the receiver and satellite clocks. 

Range observations were considered subject to time errors 
and the normal equations included timing parameters In accordance 
with the tracking scenario under consideration. The inclusion 
of such parameters weakens the normal equations as considered in 
Chapter 3. However in actual applications these parameters are neces- 
sary since tracking receiver clocks will be subject to timing offsets 
and drifts with respect to an adopted time system such as GPS system 
time. 

The solutions presented were based on two basic tracking proce- 
dures each with the adaptation of similar modeling for atomic clock 
errors. The first data acquisition procedure consisted of tracking 
satellites over three hour intervals and performing the least squares 


adjustment for station coordinates every nine hours. In this case a 
linear model was used to approximate the error in each separate clock. 
Thus the satellite rubidium clocks we re modeled by a linear function of 
time over three hour intervals and the tracking station's cesium clock 
error was modeled by the same type of function over the nine hour Inter- 
val. Solutions were performed sequentially approximately every nine 
hours with some variation if a tracked satellite's period of observ- 
ability is less than three hours. Observations were utilized only If 
the topocentric elevation angle of the satellite exceeded ten degrees. 

The second tracking scenario reduces the satellite tracking interval to 
one hour with a sequential adjustment of parameters occurring after 
four hours of observation. A similar clock modeling procedure was 
adopted but over the shorter Intervals. This latter tracking procedure 
allows a more rapid sampling of the satellite pass geometries and a bet- 
tor approximation of the random clock error; however, this procedure 
introduces a larger number of parameters of solution over a fixed period 
of site occupation. 

In all cases considered parameters representing the satellite 
orbital elements were introduced into the adjustment with a priori 
weighting consistent with the amplitude of the ephemeris error intro- 
duced, as described in Section 4.1.2. Orbital elements were introduced 
for each satellite tracking interval and corrections to these elements 
at the midpoint of the Interval were estimated as described in Section 
4.3.2. The inclusion of these parameters in the adjustment is con- 
sistent with the approach of Brown [1976] although the modeling proce- 
dure for ephemeris error is different. 
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As a final introductory comment It must be noted that the trans- 
formations between the mean Celestial System of 19,30.0 and the 
earth-fixed coordinate system are assumed known. This implies that no 
errors in precession, nutation, earth rotation, or polar motion are 
introduced into the results. The consequence of errors in these vari- 
ables is of groat Importance in geodesy but are not addressed in this 
study, Therefore in the following it is assumed that an error free 
transformation into the earth-fixed coordinate system exists. 

4.5.1 Dynamic Point Positioning 

4. 5. 1.1 Range Solutions Based on Throe-Hour 
Tracking Intervals 

A limited set of simulations based on two frequency compensated 
range observations were made using observations from Stations 1001. and 
1002 with each selected satellite of the GPS constellation tracked for 
three hours. A sequential adjustment of the earth-fixed Cartesian sta- 
tion coordinates was performed every nine hours over a five day period. 
The complete parameter set Included a linear clock model for each 
satellite rubidium clock for every three-hour Interval of tracking, 
a linear model for the receiver cesium clock for every nine-hour 
Interval, ephemerls parameter corrections for every three hours of 
tracking, and the earth-fixed Cartesian coordinates of the station. 

In addition a tropospheric scaling parameter, as described in Section 
4.1.3, was included for every three hours of observation when 
tropospheric refraction errors were introduced into the observations. 
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To evaluate the of fee. La of random and systematic error a our eon 
on station positioning, simulations were made in which individual 
error sources were Introduced into the range observations and only a 
limited number of parameters were adjusted. First, the effect of 
random cesium clock error on station positioning was examined. Simu- 
lated random cesium clock errors were developed for the cesium clock 
specifications given in Figure 4.1.3 using the algorithm of Meditch 
f 107 5 ] . This random cesium clock error was added to the. geometric, 
ranges to GPS satellite positions according to equation (4.2.8). The 
satellites were selected using the criterion of minimining the trace of 
the station covariance matrix as described in Section 4.2.2, To these 
ranges an optimistic ten centimeters of Gaussian white noise was intro- 
duced, The adjustment parameters included the Cartesian coordinates of 
the station and a first degree, polynomial in time to represent the 
cesium clock error every nine hours. The adjustment weighting was based 
solely on the white, noise statistics and the station coordinates were in 
error initially by 100 meters in latitude, Figure. 4.5.1 gives the error 
in estimated position of station 1001 as a function of time with a 
sequential adjustment in station position performed every nine hours. 
With a random white noise level of ten centimeters the range observa- 
tions would predict standard errors of 0.17, 0.21, and 0.15 centimeters 
for the u, v, and w components of station position in this example. 
According to Figure. 4.4.18 adopting a linear model to represent random 
cesium clock error over nine-hour intervals would leave an unmodeled 
random residual error with approximately a 12 to 16 centimeters standard 
error. The thirteen sets of this random residual error introduced the 
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Figure 4.5,1, Effect of Random Cesium Clock Error on 
Positioning Using Range Observations 
(.Station 1001) 

station position errors in Figure 4.5,1, The residual random clock 
errors represented by the difference between the random clock error and 
the best fitting linear model are correlated as shown in Figure 4.4.10. 
The errors introduced into station positioning by these random clock 
errors, modeled ;:n a linear function, will not average as those intro- 
duced by Gaussian white noise of a similar magnitude, as evidenced by the, 
station position errors shown in Figure 4.5.1. A comparison of these 
errors with the standard errors expected by 10 centimeter Gavissian white 
observation noise, given above, indicates the level of error expected 
from unmodeled cesium clock noise. Errors of similar magnitude although 
different in their distribution were present in the results from station 
1002, The magnitude of this error plays a more critical role for the 

determination of baseline components and is discussed in Section 4.5.2, 
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The* effect of atomic clock error on at at Ion positioning also 
Includes tlie effect of unmodeled random satellite rubidium clock error. 
For the dynamic positioning approach under examination in this analysis 
it is assumed that estimates of each satellite's clock error are pro- 
vided with the ephemerides. For the current study this implies that 
these estimates will take the form of the best linear fit and will pro- 
vide a a estimate of the systematic rubidium clock error over the inter- 
val of satellite tracking utilized in geodetic positioning. In the 
current examples that represents a three-hour Interval of time. 

Assuming for the moment that over this interval the Irian and drift of 
the satellite clock are known, then the question raised is what effect 
will the unmodeled random residual satellite clock error produce in 
station positioning? To obtain an estimate of this error station posi- 
tioning simulations were made introducing this residual rubidium clock 
error into the same geometric ranges used in the previous examples. 

This random residual error was computed by differencing simulated ran- 
dom rubidium clock noise with the best linear least squares fit to the 
noise over the tracking interval. The residuals from such a fit have an 
average standard error of approximately 30 centimeters as seen from 
Figure 4. 4. 17. The rubidium clock noise simulated was consistent with 
the rubidium oscillator Allan variance given in Figure 4.1.3. Ten 
centimeter Gaussian white noise was also introduced into the observa- 
tions representing an unreallstiely optimistic level of random receiver 
noise. The adjustment parameters included the station's earth-fixed 
Cartesian coordinates. Figure 4.5.2 gives the position errors for 
station 1001 as a function of time. The error represents the magnitude 
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Figure 4.5.2. Effect of Residual Random Rubidium Clock Error on 
Positioning Using Range Observations (Station 1001) 

of the difference between actual and estimated station coordinates after 
each nine hour update . The observations were weighted using onlv 
the white noise statistics. Based on a ten centimeter standard 
error of observation the full set of range observations would predict an 
uncertainty in station position of 0.11, 0.12, and 0.09 centimeters for 
the u, v, and w coordinates. The final errors in the station coordi- 
nates were 1.3, 0.8, and 0.6 centimeters after 117 hours of observation. 
Again, scaling the predicted standard errors by 3.0, the error intro- 
duced into station coordinates by the sequences of correlated residual 
rubidium clock error averages although not as rapidly as errors 
introduced by white observation noise of an equivalent variance. In 
this example the residual rubidium clock error even though of higher 
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variance than the residual cesium clock noise produces a smaller posi- 
tioning error. The rate at which the errors in station positioning 
average will depend on the variance and correlation of the residual 
noise process, the number of noise segments introduced, and the correla- 
tions among station position coordinates and the clock modeling para- 
meters introduced into the adjustment, However, as in the previous 
example, the effect of unmodeled residual rubidium clock error on the 
determination of station coordinates is negligible. 

To further refine the estimate of station position error 
Introduced by atomic clock error sources, an adjustment of station 
coordinates and linear clock error models was made in which random 
cesium clock error, residual rubidium clock error and ten centi- 
meters of Gaussian white receiver noise were Introduced into the geo- 
metric ranges. Adjustment weighting was based on the statistical 
modeling developed in Section 4.4 including the fully correlated 
weighting due to unmodeled atomic clock errors. The results of this 
adjustment are given in Figure 4.5.3(a) through 4.5.3(c) for each 
Cartesian coordinate of station 1001. The atomic clock errors intro- 
duced into the range observations were simulated as previously described 
for the cesium and rubidium clocks under consideration. Remembering the 
results of the previous two examples where the resulting standard errors 
of station positioning based only on Gaussian white noise were extremely 
small, it can be seen that the standard error in station position com- 
ponents due to correlated atomic clock error sources ranges from 8 to 11 
centimeters after one day of observation and from 4 to 5 centimeters after 
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Figure 4.5.3. Effect of Atomic Clock Error Sources on Positioning 
Derived from Ranee Observations (Station 1001) 
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five days, The magnitude of this error will be of Importance in the 
determination of baseline components discussed In Section A. 5. 2. 

To obtain an estimate of the effect of the ephemerls error 
described In Section 4.1.2 orbit error was Introduced into station 
positioning simulations using equation (4.2.10). The nominal level of 
ephemerls error utilized throughout this study is given in Table 4.2.2 
and the assumptions regarding its distribution are discussed in Section 
4.2.1. Adjustment results for the Cartesian coordinates of Station 1001 
are given in Figure 4.5.4 where the absolute value of the coordinate 
errors are given. Parameters in the adjustment included only station 
coordinates, ephemerls error modeling being momentarily ignored. 

Ten centimeters of Gaussian white noise were again applied to the obser- 
vations and formed the basis for the adjustment weighting. The results 
indicate that the level of orbit error addressed In this analysis may 
introduce errors into station position of greater than one meter In each 
component even after five days of continuous observation. Modeling of 
the ephemerls error tends to reduce this error. These results are simi- 
lar to those obtained for Station 1002. The errors introduced into 
positioning by each error source are dominated by the effect of errors 
in the satellite ephemerldes. This error will be the limiting factor in 
the overall accuracy to which geodetic station positions may be obtained 
using GPS ranging. 

To define an upper bound on the effect of unmodeled tropospheric 
refraction error on station positioning a five percent error was assumed 
in the predicted tropospheric refraction correction based on the 
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Figure 4.5.4, Effect of Uncompensated Systematic Orbit Error on 
Positioning Ueln<* Range Observations (Station 1001) 

Hopfleld model discussed In Section 4.1.3. This error was taken with a 
constant sign. Observations below ten degrees elevation angle were 
excluded. The results of a positioning simulation for Station 1001 are 
given In Figure 4.5.5 where the adjustment Included only the Cartesian 
coordinates of the station. Ten centimeter Gaussian white noise was 
Included as before. Refraction scaling parameters discussed in Section 
4.1.3 were not Included In the adjustment. The results demonstrate that 
a constant percentage model error In tropospheric refraction of five 
percent can Introduce errors In station position varying between 8 and 
12 centimeters, If the actual modeling error had taken the form of a 
constant percentage for each observation but with a random sign varia- 
tion for each tracking interval, the error In station positioning would 
be considerably less since the sign of the station position errors from 
each Interval of tracking would have variations resulting In better 
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Figure 4.5.5. Effect of Uncompensated Systematic Error in 
Tropospheric Refraction (5?0 on Positioning 
Using Range Observations (Station 1001) 

averaging. Subtracting the mean of each component of station position 
error from the results in Figure 4.5.5 gives an estimate of 4 to 6 
centimeters of variation at 1 day and 1-2 at 5 days which could be 
expected in such a case. Again the magnitude of this error, even con- 
sidering a worst case as in this example, is small in comparison to the 
effect of ephemerls error. 

Finally, the effect of a realistic level of receiver white noise 
is assessed in Figure 4.5.6 in which the standard error of the Cartesian 
station coordinates are given as a function of time. After twenty-four 
hours of continuous observation the standard error of the solution for 
each coordinate is approximately 2 centimeters and reduces exponentially 
to approximately 1 centimeter after five days of observation. A com- 
parison of Figures 4.5.6 and 4,5.3 reveals that, in the absence of 
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systematic errors such an satellite position and tropospheric refrae* 
tlon, Improvement in the receiver range measurement noise will not 
improve the quality of station positioning since, the effect of random 
atomic clock error will dominate. 

Table A .5.1 summarizes the approximate levels of error intro- 
duced into station positioning from the error sources discussed above 
when the satellite tracking interval is three hours. 

For Stations 1001 and 1002 complete simulations of dynamic 
noinL positioning were made using range observations from three hour 
tracking intervale, The error sources introduced into the observations 
consisted of epheraeris error, satellite rubidium clock error, receiver 
cesium clock error, tropospheric refraction error, and one meter of 
Gaussian receiver white noise in accordance with Table A. 2. 2. Va* . jus 
independent sequences of random atomic clock error were utilized in the 
analysis of station positioning for both stations. Figure A. 5. 7 gives 
the standard errors and actual position errors for Station 1001 for one 
ease. The parameters of the adjustment consisted of the full set 
described above weighted according to the level of error Introduced into 
the observations. This set included station coordinates, ephemerls 
parameters for each three hour Interval, a linear error model for the 
receiver clock over every nine-hour interval, a linear error model for 
each satellite clock for every three-hour tracking interval, and a 
tropospheric refraction scaling parameter every three hours. The least 
squares adjustment algorithm Incorporated the fully correlated adjust- 
ment weighting based on random atomic clock error and the Gaussian white 
receiver noise. Initially the station's position was in error by 100 
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TIME (HR) 

Figure 4.5.6. Standard Error of Station 1001 Coordinates 
Obtained from Range Observations Subject 
to One Meter Gaussian White Noise 


TABLE 4.5.1. EFFECT OF ERROR SOURCES ON POSITIONING DERIVED FROM 

RANGE OBSERVATIONS USING A THREE- HOUR TRACKING INTERVAL 


ERROR SOURCE 


TROPOSPHERIC REFRACTION (t5%) 

EPHEMERIS 

RESIDUAL SATELLITE RUBIDIUM CLOCK ERROR 
RANDOM RECEIVER CESIUM CLOCK ERROR 
RECEIVER WHITE NOISE (1m) 


APPROXIMATE COORDINATE 
ERROR (cm) 


1 DAY 
8-12 
150-200 
6 
8 

1. 5-2.0 


5 DAYS 
4-8* 
80-120 
2 
4 

,7-, 9 


* ASSUMING AN AVERAGING DUE TO SIGN VARIATIONS 


u COMPONENT 



TIME turn 




Figure 4.5.7. Complete Simulation of Station 1001 Positioning Using 
Range Observations over Three-Hour Tracking Intervals 
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meters in latitude',, An examination of the solution with all error 
sources Included Indicates that each component of station position can 
he determined with an accuracy of from 1.5 to 2.3 meters after one day 
of observation and from 0.8 to 1.2 meters after five days, 

The station positioning analysis based on a three-hour tracking 
interval was not Immediately extended to Include the other stations In 
Table 4.2.1, Instead consideration was given to Improving the current 
results . 

4. 5. 1.2 Range Solutions Based on One-Hour 
Tracking Intervals 

Taking Into account the results obtained in Section 4.2,2 for 
the optimal selection of satellites for point positioning, improvement 
in the geometric strength of the solution could be obtained by decreasing 
the tracking interval and sampling the satellite constellation geometry 
more rapidly. Thus a second scenario was investigated consisting of 
tracking each selected satellite for one hour and estimating station 
position every four hours. For a fixed interval of site occupation 
this approach introduces additional modeling parameters but allows 
a better sampling of satellite-station geometry. Using this approach 
ephemeris parameters are included for each hour of observation along 
with a linear satellite clock error model. A linear receiver clock 
error model is introduced every four hours. Since the clock modeling 
intervals are reduced the linear models are a better approximation to 
the random noise processes and the residual error statistics are 
reduced. However the inclusion of additional modeling parameters will 
have the opposite effect of weakening the least squares normal equa- 
tions . 
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To determine if an improvement in the prior station positioning 
results was achievable a series of solutions were made for Stations 1001 
and 1002 adopting this new approach. These positioning simulations were 
based on one hour tracking Intervals with i t :.'tal site occupation 
ranging from two to five days and were dts;f ;r,ed to measure the effect of 
random and systematic errors on positioning based on this tracking 

if 

scenario. Table 4.5.2 gives estimates of the effects of these error 
sources in a form comparable with Table 4.5.1. The magnitude of the 
errors introduced are again taken from Table 4.2.2. 


TABLE 4.5.2. EFFECT OF ERROR SOURCES ON POSITIONING DERIVED FROM RANGE 
OBSERVATIONS USING A ONE-HOUR TRACKING INTERVAL 

APPROXIMATE COORDINATE 

ERROR SOURCE ERROR (cm) 


1 DAY 5 DAYS 

TROPOSPHERIC REFRACTION (tS%) 8-12 4-6* 

EPHEMERIS 50-80 25-40 

RESIDUAL SATELLITE RUBIDIUM CLOCK ERROR 4 1 

RANDOM RECEIVER CESIUM CLOCK ERROR 5 2 

RECEIVER WHITE NOISE (1m) 1, 5-2.0 .7:9 


# ASSUMING AN AVERAGING DUE TO SIGN VARIATIONS 


An examination of Tables 4.5.1 and 4.5.2 reveals that this 
change in the observation and modeling procedure reduces the effect of 
two primary error sources, epheme'rls error and random atomic clock 
error. The effect on positioning due to residual tropospheric refrac- 
tion and receiver white noise remain virtually the same. These latter 
effects will be discussed first. 
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Simulations of station positioning wore made In which Gaussian 
white noise with a standard error of one meter was Introduced into the 
geometrical ranges to GPS satellites selected using the criterion which 
produces the smallest trace of the station coordinate covariance 
developed sequentially. The adjustment results for Station 1001 arc 
given In Figure 4. 5. 8 where the standard error and magnitude of the 
station position error are given for each component. The results are 
similar to those given In Figure 4.5.6 for the three-hour tracking 
Interval demonstrating that the effect of receiver instrumental noise 
on positioning averages equivalently for each observation procedure. 

The results indicate, that the error in each component of position due to 
receiver noise is approximately 1.5 to 2.0 centimeters after one day and 
0.7 to 0,9 centimeters after five days of continuous observation. 

For tropospheric refraction the results based on a five percent 
bias in the predicted refraction corrections show station position com- 
ponent errors ranging from 5 to 20 centimeters* With the sign of the 
modeling error taken as constant the error in the station coordinates 
appears as a bias In the range of values just given with variations 
generally on the order of five centimeters. Thus the overall effect of 
residual tropospheric refraction error remains at a level similar to 
that from the prior tracking approach. However with the introduction of 
refraction bias parameters this error Is substantially reduced. Figure 
4.5.9 gives the results of an adjustment with range observations subject 
to a systematic tropospheric refraction error of five percent and ran- 
dom Instrumental noise with a one meter standard error. In addition to 
the Cartesian station coordinates refraction scaling parameters, as 
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8. Position Error for Station 1001 Due to One Meter Gaussian 
White Noise (Range Observations/One-Hour Tracking Interval) 
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9. Position Error for Station 1001 Due to a Five Percent 
Systematic Tropospheric Refraction Error. Adjustment 
Includes Tropospheric Scaling Parameters (Range 
Observations/One-Hour Tracking Interval) 
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given in equation (4.1.57), were Included, one for each hour of 
tracking , The errors in positioning due to refraction were reduced to 
u level of approximately 4 to 6 centimeters after one day and 2 to 4 
centimeters after five days of observation. Thus refraction errors 
will not play a critical role in the determination of earth-fixed 
coordinates from GPS range observations. 

The adoption of a one hour satellite tracking interval with 
linear modeling of random clock error over shorter time intervals 
decreases the effective error In station positioning as mentioned 
earlier. After a linear approximation of random cesium clock error over 
a four hour interval the unmodeled correlated residual errors remaining 
have standard errors of approximately 9 centimeters compared to the 12 
to 16 centimeter standard error after an eight hour linear approxima- 
tion. Similarly residual rubidium clock noise over a one hour interval 
has a standard error of approximately 12 centimeters compared to approx- 
imately 30 centimeters for a three hour fit interval. Thus the expected 
magnitude of the unmodeled clock error will decrease with this alternate 
tracking approach. However the number of model parameters required in 
the adjustment will increase tending to weaken the normal equations for 
station position. Figure 4.5.10 gives an example of the errors in 
station position when random cesium clock error and instrumental 
receiver noise with a standard error of one meter are present in the 
observations. The adjustment parameters included station position and 
a linear receiver clock model for each four-hour interval. The adjust- 
ment weighting was developed using the statistics of the two random 
error sources. Taking into account the results given in Figure 4.5.8, 


197 


COORDINATE ERROR I METERS ) COORDINATE ERROR I METERS I COORDINATE ERROR fMETERSj 


Figure 4 





5.10. Position Error for Station 1001 Due to Random 

Cesium Clock Error (Range Observations/One-Hour 
Tracking Interval) 
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Figure 4.5.10 Indicates that random cesium clock errors, consistent with 
the stability specifications adopted for the receiver oscillator, intro- 
duce approximately 5 centimeters of error in each component of station 
position after one day of observation and 2 centimeters after five days. 

The random rubidium clock error modeled os a linear function of 
time over one-hour Intervals was also considered. This error source 
introduces errors of approximately 4 and 1 centimeters after one and 
five days of continuous observation respectively. With both atomic 
clock random error sources and random Instrumental noise Included in the 
adjustment, using the complete statistical weighting, the standard 
errors of station position were reduced to approximately 60 percent of 
the error present in the three-hour tracking procedure. 

With the selection of a satellite occurring each hour the 
effects of ephemerls error, whose distribution is discussed in Section 
4.2.1, averages to a greater extent than in the three-hour tracking 
scheme. Figure 4.5.11 gives an example of the errors in positioning 
expected from range observations subject to one meter random instrumen- 
tal error when ephemerls errors are present. Adjustment parameters 
include station position and six orbital elements for each one-hour 
Interval. A priori weighting consistent with the amplitude of 
ephemerls error Introduced was included for the orbital elements. The 
expected error in position due to the level of ephemerls error outlined 
in Table 4.2.2 is given in Table 4.5.2 to be 50 to 80 centimeters after 
one day and 25 to 40 centimeters after five days of observation. 

Finally in Figure 4.5.12 results are given for Station 1001 for 
a complete simulation of station positioning in which all error sources 
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Figure 4.5.11. Position Error of Station 1001 Due to Ephemeris Error 
(Range Qbservatlons/One-Hour Tracking Interval) 
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5.12, Range Positioning Results for Station 1001 for 
Complete Simulation (One-Hour Tracking Interval) 


from Table 4,2,2 were Included, All modeling parameters were Included 
in tbc adjustment and the full weight matrix booed on all unmodclcd ran- 
dom error sources was utilized. These results Indicate that the 
expected error In the components of station position range from 70 to 
140 centimeters after one day of observation and from 35 to 60 centi- 
meters after five days. A comparison of Figures 4,5.7 and 4,5.12 gives 
the improvement obtained using the shorter tracking interval. The 
Improvements in the accuracy of the recovered station coordinates for 
Station 1001 were approximately 125, 90, and 30 centimeters for the u, 
v, and w coordinates respectively after one day of observation and 65, 
55, and 20 centimeters respectively after five days of continuous obser- 
vation, Similar gains In accuracy were achieved for Station 1002. 

Since the adoption of the shorter tracking Interval produced a 
significant Increase In the accuracy of the recovered station position, 
simulations of dynamic point positioning were made for all stations In 
Table 4.2.1. These simulations Incorporated all error sources from 
Table 4.2.2, the full set of modeling parameters with a priori weights 
consistent with the level of error Introduced, and the weighting pro- 
cedure developed in Section 4,4.3 for single station tracking. Table 
4.5.4 presents the uncertainties In the geodetic coordinates for all 
stations under Investigation obtained from dynamic point positioning 
using range observations. Table 4,5.3 is provided as a key for tables 
presenting simulation results. For the adopted levels of systematic 
and random errors utilized these results indicate that the geodetic 
coordinates may be recovered to the 100 to 150 centimeter level or 
better after one day of continuous OPS range observations. After five 
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days the accuracy of the recovered coordinates is between 10 to 70 
centimeters. The dominant error source in these results is the satel- 
lite ephemerls. 

Some variations in the results ure evident. For instance, the 
solutions for the polar Stations 1009 and 1010 have a larger standard 
error for latitude and a smaller uncertainty in longitude and height 
than the results obtained in the mid-latitude station group solutions. 
This difference can he explained by examining the change in 
s r a t i on -sat e Hit e geometry. For high latitude stations the maximum 
elevation angle is considerable less. Up to a latitude of 63 degrees 
satellite crossings of the zenith are possible. However for higher 
latitudes the maximum elevation angle decreases to approximately 54.5 
degrees meaning that a larger percentage of the observations will he at 
lower elevation angles. As noted in' Section 4.2.2 increasing the num- 
ber of lower elevation observations Increases the strength of the height 
solution in the presence of timing errors. With lower elevation angle 
observations the strength of the latitude and longitude components will 
depend on the distribution of observing azimuths. The equifcorlal sta- 
tions show a larger uncertainty in height and a lower uncertainty in 
latitude and longitude, again due to the distribution of observing ele- 
vations and azimuths . The Increased frequency of higher elevation obser- 
vations is reflected in the increased height uncertainty. Figure 4.5.13 
gives the positioning results for Station 1011. 

For these adjustment solutions the a posteriori variance of unit 
weight was computed from equation (4.3.32), The square root of this 
quantity, CTq, for the solutions given In Table 4.5.4 ranged from 0,879 
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Range Positioning Results Cor Station 1011 for 
Complete Simulation (One-Hour Tracking Interval) 


to 0.914. For ft leant squares adjustment In which the mathematical 
model for the observation equation Is exact and the second order sta- 
tistics of the raidom processes are modeled correctly In the weighting, 
the theoretical value of the a posteriori variance is unity. Deviation 
from unity is primarily due to error in the above assumptions. The 
range adjustment results riven in Table 4.5. 4 are not scaled by this 
quantity. 

Some specific reasons for the range adjustment a posteriori 
variance not being unity are the following. First, the. Markov process 
transfer functions given in Figures 4.4.7 and 4.4.8 both assign more 
power to certain frequencies than the specified transfer functions given 
in the same figures. For the satellite rubidium oscillator this addi- 
tional power is at frequencies whose wavelength is greater than 100 
seconds. The actual clock noise, sequences simulated using the Medltch 
[1975] algorithm do not contain the same power at these frequencies. 

Thus the second order range statistics will predict observation uncer- 
tainties in excess of their value based on the exact use of the speci- 
fied transfer function. This tends to decrease the a posteriori 
variance. Secondly, the errors introduced into the epbemerls using the 
equations of Section 4.1.2 are periodic in mean anomaly but modeled by 
a constant amplitude correction at the midpoint of each tracking inter- 
val. This modeling difference affects the a posteriori variance 
since the level of error introduced into the observations was smaller 
than the a priori orbital element uncertainty. And finally with a small 
number of degrees of freedom for each tracking interval white and cor- 
related noise sequences will tend to he. fit better than expected causing 
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a decrease in the a posteriori variance. The interpretation of the 
results should take these factors into consideration. 

Finally Table A. 5. 5 gives the decrease in the standard error of 
geodetic coordinates obtained in two previous examples when the variance 
of the Instrumental white noise is decreased and the tropospheric 
refraction is compensated completely. The first case shows that no slg- 
nlficant increase in the accuracy of the adjusted station coordinates 
can be expected by decreasing this Instrumental random error component. 
Decreasing the standard error of this component from 100 to 60 centi- 
meters produces a decrease of only 1.1 centimeters or less as opposed 
to an expected decrease of AO percent based on range measurements sub- 
ject to white measurement noise only. In this case however, with the 
inclusion of the fully correlated statistics for unmodeled atomic clock 
error, the resulting decrease Is marginal. 

In the second case assuming that tropospheric refraction effects 
can be compensated the refraction scaling parameters are excluded from 
the adjustment. The decrease in the standard error of the geodetic 
coordinates ranges from 1.1 to 3.0 centimeters after one day of obser- 
vation and from 0.8 to 1.6 centimeters after five days. The largest 
decrease is in the height uncertainty; although, the net effect on the 
determination of earth-fixed coordinates is minor. 

A. 5. 1.3 Integrated Doppler Solutions 

Solutions based on integrated Doppler or range difference obser- 
vations were examined subsequently. Range differences over five-minute 
intervals were formed by aggregating independent one minute integrated 

Doppler observations with an instrumental measurement uncertainty of 
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three centimeters. The parameters of the adjustment were equivalent 
to the range solution set with the exception of the clock models. For 
integrated Doppler observations the linear clock error models were each 
replaced with a single parameter representing time drift or frequency 
bias. Hphemeris elements, tropospheric refraction corrections, and the 
geodetic station coordinates were retained. 

An Initial solution was made for .Station 1001 using range dif- 
ference observations over three-hour tracking intervals. Observations 
were simulated using equation (4.2.9) and ephemerls error was intro- 
duced into the adjustment using an equation analogous to equation 
(4.2.10). Error sources were taken from Table 4.2.2. Adjustment 
weighting Included both the Instrumental white noise statistics and the 
random clock error statistics developed in Section 4.4.2. The receiver 
clock was modeled over a nine-hour interval as in the range solutions 
based on the same interval of tracking. Table 4.5.6 gives the results 
for tills adjustment. These results Indicate that after one day of 
observation the geodetic coordinate errors can be expected to range from 
125 to 215 centimeters and reduce to from 60 to 100 centimeters after 
five days. 

For this station a simulation based on a ore-hour tracking 
interval was next tried to determine if better results could be obtained 
as in the range case with a receiver clock model adopted every four 
hours. The results from this solution are given at the beginning of 
Table 4.5.7. A comparison of the three and one-hour tracking Interval 
results shows that significant improvement is obtained with the shorter 
tracking interval. This latter tracking procedure allows a better 
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representation of random clock error with the same model, permits better 
sampling of satellite geometry, and produces a more rapid averaging of 
the effects due to systematic and random error sources. The result is 
an uncertainty in station position coordinates in the range of 100 to 
1-50 centimeters after one day and 45 to 05 centimeters after five days 
of observation. 

Since this shorter tracking interval yielded such Improvement In 
the results, station positioning adjustments were made for most of the 
stations in Table 4.2.1. These results are given in Table 4.5.7. Again 
the height uncertainties for the polar Stations 1009 and 1010 are signi- 
ficantly less than for all other stations since the higher occurrence 
of lower elevation observations allows a better separation of height 
and timing errors. However for these stations the latitude and longi- 
tude solutions arc weaker. The results for the mid-latitude stations 
show less variation than the range solution results. In general the 
results Indicate that range difference observations yield position com- 
ponent accuracies of from 85 to 200 centimeters after one day of obser- 
vation and from 40 to 80 centimeters after five days of continuous 
tracking. Variations in the results with location are to be expected 
with the weakest solution for latitude and longitude occurring toward 
the poles. Figure 4.5.14 gives the positioning results for Station 
1013. 

The a posteriori variance of unit weight was computed using 
equation (4.3.32) for each adjustment of Table 4.5.7. The square root 
of this quantity varied from 0.967 to 0.998 indicating more consistency 
in the adjustment modeling and weighting than in the range observation 


simulations . 
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Table 4.5,8 nummarlzoo the effects of various error sources on 
positioning baaed on Integrated Doppler observations using n one-hour 
tracking Interval. Again the ephemerls error dominates the effects of 
all other error sources. Because of the geometric weakness of the Inti 
grated Doppler observations the effect produced by three centimeter 
Instrumental white noise Is much larger than that due to a one meter 
standard error In range. 

TABLE 4.3.8. EFFECT OF ERROR SOURCES ON POS ri'TONTW ITO”RP 
FROM DOPPLER OBSERVATIONS USING A ONE-IIOUR 
TRACKING INTERVAL 

ERROR SOURCE 

TROPOSPHERIC REFRACTION (5%) 

EPHEMER1S 

RESIDUAL RANDOM RUBIDIUM 
CLOCK ERROR 

RANDOM RECEIVER CESIUM 
CLOCK ERROR 

RECEIVER WHITE NOISE 

"ASSUMING AN AVERAGING DUE TO SIGN VAI 
REFRACTION SCALING PARAMETERS WOULD 
APPROXIMATELY D CENTIMETERS. 

And finally Table 4.5.9 gives the reduction In the geodetic 
[f coordinate uncertainties with modifications to the assumed error levels 

introduced into the adjustment for Station 1007. Assuming that tropo- 

I 

t spheric refraction can be accounted for completely either through 

I measurement or modeling and that the scaling parameters are deleted from 



215 
































































the sot of adjustment parameters, Che resulting decrease In position 
uncertainty in negligible ns shown in She cable. Also a decrease in the 
instrumental noiBe level from chrec centimeters to 0.2 conCimeterB pro- 
duces only a minimal reduction in the coordinate uncertainty. Minor 
decreases are realised in these coses because the ephcmcris error 
totally dominates these error sources. Thus for absolute positioning 
additional refinements in the refraction prediction or Improvements in 
the noise level of the receiver will not provide any real Improvement 
unless the ephemeris error is greatly reduced, Asa final example the 
level of ephemeris error adopted in Table 4. 2. 2 was halfed and as expected 
a significant level of improvement in position uncertainty was achieved. 
The uncertainty in the results improved by approximately 45 percent. 

A comment concerning continuous count integrated Doppler is in 
order. In the above analysis one* minute integrated Doppler counts, 
assumed statistically independent, were aggregated to form five minute 
range differences. The instrumental noise thus increased by /5. For a 
continuous count integrated Doppler system this is not true. The five 
minute Doppler counts in that case would still be subject to approxi- 
mately the same white noise level as one minute observations. The 
results presented here consider one-minute observation noise levels of 
from 0.2 to 3 centimeters or 0.45 to 6.7 centimeters for five minute 
aggregated range differences. For continuous count integrated Doppler 
this latter interval would be approximately 0.2 to 3 centimeters, a 
more optimistic but partially overlapping interval. From the results in 
Tables 4,5.7 and 4.5.9 the accuracy of continuous count Doppler utilized 
as independent range differences can be established. 


217 


4.3.2 Determination of BaB eline Components 


In this section results arc presented for the determination of 
baseline components and chord length from simultaneous range and Inte- 
grated Doppler observations from two stations. The least squares normal 
equations now Include the earth-fixed coordinates of each tracking sta- 
tion. After each sequential solution the resulting station coordinate 
covariance matrix is linearly transformed into coordinate differences 
and chord length using equations (3.1.22) and (3.1.25) or analogous 
equations when the coordinates are expressed is geodetic latitude, 
longitude, and height. Satellites are selected using the criterion dis- 
cussed In Chapter 3 and the simulations described in this section 
include Table 4.2.2 error sources. 

4. 5. 2.1 Range Solutions 

As In Section 4.5.1 Initial results were based on the three-hour 
tracking interval. Simultaneous range observations were simulated 
for five days for Station 1001 and 1002. These stations lie on the same 
meridian separated by approximately 100 kilometers as shown In Figure 
4.2.1. Simultaneous observations were excluded from the adjustment if 
the elevation angle from either station was below ten degrees. 

The effect of Individual error sources on baseline components 
was Investigated for these stations by introducing each into the 
adjustment. For this 100 kilometer north-south baseline the results are 
given in Table 4.5.10. A comparison of these results with Table 4.5.1 
demonstrates that the sensitivity of the baseline components to these 
error sources is quite different than for the determination of geodetic 
coordinates from range. Since the baseline distance is small relative 
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TABLE .'i.'j.lO. 


EFFECT OP ERROR SOURCES ON ISAS KUNE COMPONENTS 
DERIVED PROM RANCH OBSERVATIONS USING A 
THRESH -HOUR TRACKING INTERVAL (100 km »ASKIJm v 


ERROR SOURCES 

APPROXIMATE 
COMPONENT 
ERROR (cm) 


1 DAY 

5 DAYS 

TROPOSPHERIC REFRACTION (5%) 

6-8 

2-3 

EPHEMERIS 

1-3 

0.5-1.5 

RESIDUAL RANDOM SATELLITE 
RUBIDIUM CLOCK ERROR 

0.2 

0.1 

RANDOM RECEIVER CESIUM 
CLOCK ERROR 

12 

6 

RECEIVER WHITE NOISE (1 m) 

2-3 

1-1.5 


to the distance to the satellites the effects of errors in the satellite 
ephemeris and clock project almost Identically into the coordinates of 
each station. The transformation into coordinate differences removes 
the majority of the effect. Thus a 1 though satellite position errors 
can contribute 150 to 200 centimeters of uncertainty in station position 
after one day of observation, this same error has only an effect of from 
1 to 3 centimeters on the coordinate differences. This fact precludes a 
requirement for a precise, ephemeris in this application, Figure A, 5. 13 
demonstrates the error in the Cartesian baseline components due to 
ephemeris error. After five days this error can be expected to range 
from 0.5 to 1.5 centimeters. The effect of the satellite clock error 


is likewise minor. 
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Filin' 4 . 5 . 1 5 . Kffoet of Uncompensated Systematic Orbit, Error on 
Baseline Component a Using Range Oh so rvatiovn; 
(Three-Hour Tracking Interval) 


The of foot of tropospheric refraction error will also bo 
reduced if the signature of the error is almost equivalent at both 
sites. In the eases considered here a common five percent error was 
Introduced. For sites separated by up to a few hundred kilometers the 
difference in tropospheric refraction w ill he primarily a function of 
elevation angle difference and the difference in weather conditions. 
Assam inj', the difference is a function of the former, a constant percetv 
tage error will produce approximately the same error at each site and 
the effect on baseline components will be small. In actual applica- 
tions where a more complicated prediction of tropospheric refraction 
exists the baseline component errors may increase to a value greater 
than that given in Table 4.5.10. 

one 

4- ** VJ 


The dominant error source In this application is the Instability 
of the tracking receiver clock. For the cesium oscillator considered In 
tills study the error Introduced Into baseline components can be expected 
to range from 12 to 6 centimeters after one and five days of observation 
respectively. Figure 4. ‘5.16 presents the Cartesian baseline component 
errors as a function of time. These errors tend to average with time 
but at a rate which depends on the stability of the clock. For the 
dynamic determination of baseline components a significant decrease in 
this error can only be achieved by Increasing the stability of the 
receiver oscillator If the tracking Interval is held fixed. 



BASELINE 1001-1002 


TIME (HR) 

Figure 4.5.16. Effect of Random Cesium Clock Error on Baseline 
Components derived from Range Observations 
(Three-Hour Tracking Interval) 


For a complete simulation of relative positioning five days of 
range observations were adjusted sequentially. The parameters of the 
adjustment included two pairs of station coordinates, ephcmerla para- 
meters, a linear satellite clock model every three hours and a linear 
model for each tracking station clock every nine hours, Two trope- 
spheric refraction sealing parameters were introduced for every three 
hours of tracking, The adjustment weighting was based on all random 
errors added to the observations. The weight matrix used every nine 
hours had the form of equation (4.4,41) , Errors were Introduced into 
the range observations according to Table 4.2,11 with the Instrumental 
white noise uncertainty taken as one meter. The result of the adjust- 
ment was Cartesian base, Hite coordinate uncertainties of 23, 21, and ,18 
centimeters after one, day of observation and 10, 11, and 8 centimeters 
after five days. 

The simulation for Stations 1001 and 1002 was repeated using a 
one-hour tracking interval. The uncertainties in the Cartesian base- 
line components after one day of observation were 28, 1,1, and 12 centi- 
meters. After five days of continuous observation the resulting 
standard errors were 12, 7, and 6 centimeters for the Au, Av, Aw compo- 
nents. A comparison of the. trace of the covariance matrix with that 
from the previous three-hour interval simulation shows that the shorter 
tracking Interval produces marginally better results. This is con- 
sistent with the marginal increase in geometric strength for range 
observations demonstrated in Chapter 3. 

For the one-hour tracking procedure the effects of error sources 
on baseline components are given in Table 4.5.11 for the 100 
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kilometer baseline. The major difference between Tables 4.5.11 and 
4,5.10 is the decrease in the effect of random receiver clock noise. 


TABLE 4,5,11. 


EFFECT OF ERROR SOURCES ON BASELINE COMPONENTS DERIVED 
FROM RANGE OBSERVATIONS USING A ONE-HOUR TRACKING 
INTERVAL (100 km BASELINE) 


ERROR SOURCE 

APPROXIMATE 
COMPONENT 
ERROR (cm) 


1 DAY 

5 DAYS 

TROPOSPHERIC REFRACTION (5%) 

6-8 

2-3 

EPHEMERIS 

1-3 

0.5-1. 5 

RESIDUAL RANDOM SATELLITE 
RUBIDIUM CLOCK ERROR 

0.2 

0.1 

RANDOM RECEIVER CLOCK 
ERROR 

8 

3 

RECEIVER WHITE NOISE (1 m) 

2-3 

1 .0-1.5 


Using the one-hour tracking scenario simulations wore performed 
to assess the accuracy to which baseline components and chord lengths 
might be determined using simultaneous GPS range observations from two 
sites. The complete parameter set and weighting based mi all random 
error sources were Included in the adjustment. The resulting uncer- 
tainties in the baseline parameters are given in Table 4.5.12. The 
uncertainty in the chord length d is also expressed in parts per mil- 
lion (ppm) . For baselines less than 300 kilometers in length these 
results indicate that the uncertainty in the latitude component of the 
baseline ranges from between 10.1 and 12.7 centimeters after one day of 
observation and from 4.4 to 5.8 centimeters after five days. The longi- 
tude component uncertainties are slightly weaker ranging from 10.6 to 
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17.4 centimeters after one day and from 4.3 to 7.5 centimeters after 
five days of simultaneous observation. The height difference between 
the stations has the largest uncertainty due to the correlation of height 
error with receiver timing error and refraction. For baselines under 300 
kilometers the height difference uncertainty ranges from 25.2 to 27.7 
centimeters for one day of observation and from 12.0 to 15.0 centi- 
meters after five days. For these baselines the uncertainty in chord 
length ranges from 9.9 to 16.0 centimeters (0.5 to 1.8 ppm) after 
one day and from 4.3 to 6.9 centimeters (0.2 to 0.8 ppm) after five 
days. The chord length uncertainty Increases with baseline distance 
as seen in the results for baselines 1007-1008 and 1015-1016. How- 
ever the relative error in parts per million decreases. The increase 
in the uncertainty Is due to an increasing projection of the ephemeris 
error onto the baseline components. Figure 4.5.17 gives the baseline 
component errors and uncertainties for baseline 1011-1012. The 
chord length between these stations is approximately 100 kilometers. 
The chord uncertainty as a function of time is given in 
Figure 4.5.17(d). 

Two final examples are presented in Table 4.5.13 which show how 
the uncertainty in the results given in Table 4.5.12 are subject to 
change with variations in the simulation. Decreasing the instrumental 
white noise to 60 centimeters decreases the uncertainties of the base- 
line components from 1.5 to 2.4 centimeters after one day of observa- 
tion and from 0.6 to 1.2 centimeters after five days for the 80 
kilometer baseline 1001-1014. The decrease in the uncertainty of the 
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Figure 4.5.17. Baseline Parameter Error for Complete Simulation 

Using Range Observations from Stations 1011 and 1012 
(One-Hour Tracking Interval) 


















chord length ranges from 0,3 ppm to 0,1 ppm for the Interval of site 
occupation in this case. Finally If tropospheric refraction error can 
be measured or predicted with high accuracy then the delation of the 
tropospheric refraction scaling parameters will produce a significant 
decrease In the standard error of the baseline parameters since the 
removal of these parameters will strengthen the normal equations. The 
results for the 200 kilometer baseline 1003-1004 are given In Table 
4,3.13 demonstrating that the measurement of tropospheric refraction 
with a water vapor radiometer may bo required to obtain the best pos- 
sible results using a dynamic approach. 

4. 5. 2. 2 Integrated Doppler Solutions 

Simultaneous Integrated Doppler observations from a pair of 
stations were analyzed to determine the accuracy to which baseline 
parameters can be determined. Adopting a one-hour tracking interval 
the effect of the systematic and random error sources given in Table 
4.2.2 on the vector components of the baseline were evaluated for Sta- 
tions 1001 and 1002. These results are given in Table 4.5.14. As with 
the use of range observations the stability of the tracking receiver 
clock will contribute slgnlf Icantly to the error In this positioning 
problem while satellite ephemeris and clock errors have no significance 
for such short baselines. As mentioned previously a five percent 
unmodeled error in tropospheric refraction having constant sign can 
Introduce errors of up to 50 centimeters in position. However for short 
baselines a large portion of this error is in common at both sites and 
the resulting error In the coordinate differences ranges from 4 to 8 
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TABLE 4.5. 14. EFFECT OF ERROR SOURCES ON BASELINE COMPONENTS 

DERIVED FROM DOPPLER OBSERVATIONS USING A ONE-HOUR 
TRACKING INTERVAL (100 lcm BASELINE) 


APPROXIMATE COMPONENT 

ERROR SOURCE ERROR (cm) 



1 DAY 

5 DAYS 

TROPOSPHERIC REFRACTION 

4.Q 

2-3 

EPHEMERIS 

.5 

,1 

RESIDUAL SATELLITE RUBIDIUM CLOCK ERROR 

.2 

.05 

RANDOM RECEIVER CLOCK ERROR 

10 

4 

RECEIVER WHITE NOISE (3 cm) 

20-25 

8-10 


centimeters after one day to 2 to 3 centimeters after five days of 
observation. In actual applications the signature of this error may 
not be equivalent at each site and the resulting baseline component 
errors may be different. The receiver white noise plays the most 
important role. Because of the geometric weakness of range difference 
observations a 3 centimeter standard error for receiver noise will 
restrict baseline component uncertainties to be more than 20 to 25 
centimeters after one day of observation and from 8 to 10 centimeters 
after five days. Furthermore it will be shown below that reducing the 
receiver noise level will have only limited success in reducing the 
baseline component uncertainties. 

For the case just considered a complete simulation was made to 
determine the uncertainty In the baseline components and chord using 
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five days of continuous observation. The results arc presented In 
Table 4 .5.13. Xn terms of the geodetic coordinate differences the 
uncertainties are 27.7, 41,6, and 45.5 centimeters after one day and 
13.0, 17.1, and 18,2 centimeters after five days for the latitude, 
longitude, and height differoncco. The uncertainty In the chord was 2.8 
ppm and 1,3 ppm after one and five days respectively, Results for the 
same baseline were then obtained using a three hour-satellite tracking 
interval. These uncertainties are given In Table 4,5.16. Comparing the 
24 hour results with those obtained using 27 hours of observation froiii 
three- hour tracking Intervals demonstrates that each tracking procedure 
gives comparable results. After five, days of observation It appears 
that using a longer tracking Interval haw some advantage for deter- 
mining the chord. 

Since the three-hour tracking procedure did not appear to pro- 
duce a significant overall advantage results for other baselines were 
determined using the one-hour tracking Interval and are also presented 
In Table 4.5.15. These results Indicate for baselines less than 500 
kilometers that the latitude difference uncertainty ranges from approx- 
imately 30 centimeters after one day to 13.5 centimeters after five 
days and Is the best determine component of the baseline as in the case 
of range observation. This is due to the fact that the majority of the 
observations are from north or south going pass geometries as shown In 
Figures 4.2.4(a) and (b) . The uncertainty of the longitude component 
of the baselines ranges from 37.3 w 44.8 centimeters after one day 
and from 14,5 to 19.1 centimeters after five days. Height difference 
uncertainty ranges from 35.1 to 44.8 centimeters after one day and from 
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IS. 8 to 19. 1 centimeters after five days. The uncertainty in the chord 
is significantly smaller for north-south baselines for the same reason 
as the latitude component. This is also true i the range observation 
examples although the difference as a function of orientation is less 
pronounced after five days of observation. Finally, the uncertainty 
increases with station separation although the ratio of the uncertainty 
in the chord to its length decreases with the longer baselines given 
here. The increase, in uncertainty is again due to the increased effect 
of ephemeris error. Figure A. 5. 18 gives the results obtained for base- 
line 1015-1016. As is typical with the results from the other baselines 
considered the decrease in parameter variance appears as an exponential 
decay . 

Finally, in Table 4.5.17 various cases are considered in which 
modifications are made to the error sources. Elimination of tropo- 
spheric refraction produces a decrease in the baseline component uncer- 
tainties ranging from 1.0 to 7.3 centimeters after one day of 
observation to 0.7 to 2.7 centimeters after 5 days. The height uncer- 
tainty is decreased to the greatest extent. The chord uncertainty 
decreases by0.2ppm after one day of observation and by 0 .1 ppm after an 
additional four days of observation. Decreasing the ephemeris error 
by 50 percent produces only minor variations in the results as expected. 
And finally the last two cases of Table 4.5.17 show that reducing the 
random receiver noise to 1 centimeter produces a significant increase in 
accuracy but Improvement beyond that level gives only limited success 
since the effect of random receiver clock error begins to dominate the 
resulting parameter uncertainties. 


234 


METERS 


LATITUDE DIFFERENCE 



o 20 to eo ao ioo 120 mo ito 

TIME IHRl 



TIME (HR! 


Figure 4.5.18. Baseline Parameter Error for Complete Simulation 
Using Integrated Doppler Observations from 
Stations 1015 and 1016 





























5. A PRELIMINARY EVALUATION OF SATELLITE 
INTERFEROMETRY FOR BASELINE DETERMINATION 

5 . I Introduction 

In the preceding chapter the accuracy of baseline determinations 
from range and Doppler observations was analyzed considering the effects 
of various error sources. It was found that the range results were pre- 
dominantly influenced by tropospheric refraction modeling error and ran- 
dom receiver clock error while Doppler results were influenced most by 
the same and by random receiver noise. Tropospheric refraction errors 
may be reduced by the use of a water vapor radiometer [MacDoran, 1979] 
and Doppler receiver noise levels may actually be as low as one centi- 
meter [Stanford Telecommunications, Inc., 1978], hence the baseline 
uncertainties obtained from range and Doppler may be enhanced as demon- 
strated in Tables 4.5.13 and 4.5.15. However the resulting baseline 
uncertainties would still be effected by random correlated clock errors 
and, in the case of Doppler, also by the weaker geometric strength of 
the observations themselves. Accuracies on the order of 1 ppm may be 
achieved using these methods if the period of site occupation is at 
least 2 days for range and 5 days for Doppler observation. 

Since neither of these two observational approaches will sup- 
port a rapid first-order determination of baselines, this chapter is 
included to address the utilization of interferometric phase measurements 
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for this application. Although the interferometric proposals discussed 
in Chapter 2 are currently under development, enough information is 
available to support a general estimate of the performance of an inter- 
ferometric approach. The technique examined in this chapter is based on 
the double differencing of interferometric phase from two satellites 
made simultaneously at two sites, This approach has the advantage of 
eliminating most of the clock errors which required polynomial modeling 
in the range and Doppler approaches. The analysis presented here is of 
a preliminary nature Intended to provide a general estimate of the 
accuracy of baseline determination using interferometry. A more 
detailed analysis of the proposed Interferometric procedures of Chapter 
2 should be performed as the specifics of these techniques are refined. 


5 . 2 Double Differencing of Interferometric Phase 

The approach which is Introduced in this section assumes that 
Interferometric phase observations are based on the reconstructed con- 
tinuous wave GPS carrier frequencies. The following observational model 
is adopted for the phase measurement with station 1 observing satellite 

j: 


°« (c) “ *770 V t> -V t> V ,) -° {l i (l:> + e{t j <t) 

- 3<5R i;j (t) + Yj j (t) ] . 


(5.2.1) 


In this equation X is the wavelength of the GPS carrier frequency, 

is the geometric distance between station i and satellite j, m^ is 
the Integer number of full wavelengths comprising R ij * -i 


ij 
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tropospheric refraction error modeled to within a percent 3, and y 
represents uncorrelated measurement error. The quantities eSt, and 
c6tj represent systematic and correlated random time or phase errors, 
converted into units of length, of the geodetic receiver and satellite 
atomic clocks respectively arising from the accumulation of fractional 
frequency error. The wavelength is also a function of time due to 
the Doppler shift caused by the relative motion of the satellite with 
respect to the receiver. Hence, 

X i 

X (t) P J — (5.2.2) 

ij 1 - p/c 

where X^ is the carrier frequency and p is the component of relative 
velocity along the topocentrlc range vector. 

If satellite j is simultaneously observed at station & then 
the difference in phase measured at the two sites is 

A0 (t) « 0^(0 - 0^(t) . (5.2.3) 

Ignoring for the moment the Doppler shift in the carrier frequency and 
assuming the same level of refraction modeling error at both sites, 
equation (5.2.3) may be written as 


A0j(t) 


•xKj (t) ' V° ‘ Kj (t) 

i (5.2.4) 

- cftc^lO-St^e)] -UfrRyO:) - [Y^CO-Y^OO]}. 


Notice that the error in the satellite clock does not appear in 
equation (5.2,4) due to the differencing. 
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If In addition a second satellite k Is simultaneously observed 


at both stations then the double difference is defined as 


AO ik (t) - AOj(t) - AO k (t) 

= - o tJ (« - o lk (t> + o^CO . 


(3.2.5) 


Again ignoring the Doppler shift and assuming that the frequencies A^ 

and A, are equal and that 0 Is constant for all observations, equation 
k 

(3.2.5) may be written as 


“ T s {V° * V° ’ R ik ( 0 + V t) + V (t>x J 

- s [61^(0 -MyCt) 

where 

";|k “ -*lj + + ”'ik - "*k 

and 

Y Ujk 53 Y ij ~ Y £j “ Y lk + Y ik ’ 


(5.2.6) 


(5.2.7) 


(5.2.8) 


In equation (5.2.6) no atomic clock errors appear; thus, the double dif- 
ferencing approach appears to eliminate the timing errors which required 
modeling previously. The Integer term n j k (t) represents the difference 
between a pair of "2k ambiguities" which exist in each single dif- 
ferencing of phase according to equation (5.2.4). This ambiguity repre- 
sents the integer number of full wavelengths comprising the difference 
in the distances between the stations and the satellite. 

If the Doppler shift in frequency is Included, then the double 
difference equation (5.2.6) would be replaced by the substitution of the 
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appropriate equations (3,2,1) Into equation (3.2,3) and a complete can- 
cellation of clock error could not be expected. In the preliminary 
analysis presentee In this chapter the carrier frequency A^(t) will be 
assumed known and equation (3.2.6) will be adopted as the observation 
equation. Introduction of the Doppler shift will some cause additional 
uncertainty in the results depending on the a priori errors in station 
position and satellite position and velocity. 

An additional assumption implied in equation (5.2.6) is that 
simultaneous observations of phase are to be differenced. The recogni- 
tion of simultaneous events depends on accurate time tagging of the 
observations or knowledge of the relative time error between station 
clocks. The first of these Is impossible to achieve and the latter 
requires either portable clock comparisons or the adoption of additional 
parameters In the estimation algorithm. If phase differences are formed 
from observations at two sites having a time of observation difference 
of At seconds, then the error introduced into the double difference is 
given approximately by 
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(5.2.9) 


This equation is obtained from a first order Taylor series expansion of 
equation (5.2.6) assuming that the observations selected from station H 
for differencing are At seconds away from those from station i. An 
examination of equation (5.2.9) reveals that the time synchronization 
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error condo to of a geometric term duo to the position change of each 
satellite relative to station A, a term which is a function of oatel~ 

11 to clock frequency stability and a term due to the variation in tropo- 
(spheric refraction over At, An evaluation of this equation for a 
synchronization error of 200 nanoseconds gives a bound on this error of 
0.03 centimeters. Synchronization to much better than this leve3 could 
be achieved by time tagging observations with the satellite time infor- 
mation encoded in the transmitted GPS signals, 

Adjustment of baseline parameters using double differenced 
interferometric phase observations requires the differentiation of 
equation (5.2,6) with respect to the earth-fixed coordinates of the 
observing stations, the integer n.j c and the constant 8, at a minimum. 
Satellite position also enters into equation (5.2.6) and represents an 
additional set of parameters which strictly should be Included. In the 
results given below corrections to the satellite ephemerldes are not 
Incorporated but the effect of error in satellite position is dis- 
cussed. The partial derivatives used to form the design matrix for the 
least squares adjustment are the following 


8A9 i k (t) 2nf u i (l:> - u l u k (t) - U 1 

3u l _ R ijC) R ik^ 

u ■+ V, w 


3Ae ik (t) 2„ K - U 1 (t) “* - \ (t) 

Su 2 M V t} I W t > 

u *+• V, w 


(5,2.10) 


(5.2.11) 
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The Integer n.. for nn observed satellite pair Is a function of 
j k 

time, changing at each observation time. If the receiver however main- 
tains a count of accumulated phase change over the tracking Interval, 
then the rate of change of n^ is known and only a single Integer 
unknown needs to be incorporated for each Interval of tracking. An 
adjustment based on equations (5.2.10) through (5.2.13) will not pro- 
duce Integer solutions for the n^. Since no constraints ore known 
which will produce an Integer result directly, this initial adjustment 
will provide a set of estimates and variances for the n^. From these 
quantities various test sets of Integers may be formed. The number of 
such sets will depend on the estimates of the n^ and on the magnitude 
of their corresponding variances. For each test set a second least 
squares adjustment would be required utilizing these integers. This 
second adjustment would include a set of absolute constraints fixing 

the n.. . From these adjustments the weighted sum of squares of residuals 
jK 

V T PV may be compared to determine which test set of integers produces a 
minimum, The covariance matrix of the station coordinates from this 
solution may be transformed using equations (3.1.24) or (3.1.25) into 
baseline component and chord length uncertainties. 


r i,3 Effect of Error Sources 


The error sources influencing satellite Interferometry measure.'* 
meats were mentioned tn Chapter 2. For the double differencing approach 
the error sources will be the same except that it appears that most 
error due to frequency instability will be removed, The error sources 
considered In the results presented here are the satellite cphcmerldco , 
tropospheric refraction, and the random error associated with the 
measurement of phase. The magnitudes of the ephemeris and tropospheric 
refraction errors are equivalent to those used in the range and Doppler 
positioning studies as outlined in Table A, 2. 2, The precision of a 
single phase measurement is assumed to be 3 centimeters which was the 
nominal precision adopted for integrated Doppler observations in Chapter 
A, Counselman [1979] estimates the random phase error of the Miniature 
Interferometer Terminals to be less than 1 centimeter. Table 5.3.1 
gives estimates of the effects these error sources have on baseline com- 
ponents and chord length for sites separated by 100 kilometers. These 
results are based on simulations using a total of six hours of observa- 
tion, tracking individual satellites for a fixed one-hour interval. 

TABLE 5.3.1. EFFECT OF ERROR SOURCES ON BASELINE PARAMETERS DERIVED 
FROM SIX HOURS OF DOUBLE DIFFERENCED INTERFEROMETRIC 
PHASE USING A ONE-HOUR TRACKING INTERVAL (100 km 
BASELINES) 

COMPONENT ERROR CHORD ERROR 


ERROR SOURCE 

(cm) 

(cm) 

TROPOSPHERIC REFRACTION (5%) 

2-4 

U i 

EPHEMERIS 

US 

2-3 

RECEIVER WHITE NOISE (3 cm) 

U4 

1 
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5.4 Baseline Determination Results 


Double differenced Interferometric phase observations were simu- 
lated every five minutes for the baselines previously considered in 
Chapter 4. These observations were developed using equation (5.2.6). 
Three initial adjustments were performed using observations from sta- 
tions 1001 and 1002 simulated using a one-half, one, and two hour 
satellite tracking interval. Satellite position error and a five per- 
cent error In tropospheric refraction modeling were Introduced into the 
adjustment. Parameters of the adjustment included -he latitude, longi- 
tude, and height of each station, the Integer n^> and the constant ft for 
each tracking interval. The uncertainty of the latitude, longitude, and 
height components of the baseline and of the chord length obtained from 
these initial adjustments are given in Table 5.4.1 after six hours of 
observation. As the fixed interval for observing a pair of satellites 
is increased from one-half to two hours, there is a marked decrease in 
the parameter uncertainties except for the height component. However, 
w i th additional observations this trend is apparent for height also . 
After ten hours of observation the height component uncertainties are 
7.8, 4.6, and 3.2 centimeters for the three intervals utilized. This 
trend is due to the decrease in the total number of parameters required 
in the. adjustment as the tracking interval is lengthened resulting in a 
general strengthening of the normal equations. 


TABLE ,'i .4.1 . 


VARIATION IN BASELINE PARAMETER UNCERTAINTY WITH SATELLITE 
TRACKING INTERVAL (BASELINE 1001-1002, SIX HOUR RESULTS, 
INITIAL ADJUSTMENT) 


TRACKING INTERVAL 

°A 0 


°Ah 


0.5 hr 

9.3 cm 

12.2 

11.4 

5.6 

1 

t. 3 

8.1 

7.5 

3.1 

2 

2.5 

4.6 

9.0 

1.8 


For the solution based on a one-hour observing interval Table 
5. A. 2 gives the actual and estimated values for the Integers and the 
uncertainty of their solution. It is typical in the shorter tracking 
interval cases for the uncertainty of the estimated n^ t0 exceed 0.5. 
When this occurs, the number of test sets of integers required in sub- 
sequent adjustments may be large. For instance in Table 5.4.2 the 
solution for n^ for the fourth hour of observation was -15.2. With 
the standard error of this solution 0.71 any of the following integer 
values, -13, -14, -15, -16, -17, could be expected as the correct solu- 
tion for this Interval. If all solutions lying within a 95 percent con- 
fidence Interval are considered, the number of possible unique sets of 
integers to be used in subsequent adjustments may be extremely large. 


TABLE 5.4.2. RESULTS OF ADJUSTMENT FOR INTEGERS njk BASED ON INTER- 
FEROMETRIC PHASE MEASUREMENTS AT STATIONS 1001 AND 1002 
USING ONE-HOUR TRACKING INTERVAL 

INTEGER ESTIMATE UNCERTAINTY 

TRACKING INTERVAL n jk h jk a njk 


1 

7 

7.2 

.16 

2 

-6 

- 5.8 

.15 

3 

1 

0.5 

.16 

4 

-16 

- 15.2 

.71 

5 

-13 

-13.0 

.24 

6 

-11 

- 10.7 

.54 
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Since initial adjustments using observations from two-hour 
tracking intervals produced the smallest uncertainty in both baseline 
components and more Importantly in the integers n^, solutions were made 
for all baselines considered in Chapter A using observation schedules 
based on two-hour tracking Intervals. Ephemeris, t , opo spheric refrac- 
tion, and Instrumental errors were added to the observations. The 
satellite pairs were selected for tracking to optimize the trace of the 
baseline parameter covariance matrix given by either equation (3.1.25) 
when the chord was estimated or by an equation analogous to equation 
(3.2.3) when latitude, longitude, and height components were estimated. 
After an initial adjustment the same observational data were utilized in 
a subsequent adjustment in which the correct integer values n^ were 
included and fixed by absolute constraints. In actual practice many 
such solutions may be required. The results of the second adjustment 
are given in Table 5. A. 3. The results after six hours of observation 
indicate that the uncertainty of the baseline components generally 
ranges from be tween 1.0 and A.O centimeters for baselines of 100 kilo- 
meters. These uncertainties increase with baseline length. For shorter 
baselines the height component has the largest uncertainty. The accu- 
racy of the chord length exceeds 0.1 ppm in all cases considered with 
the relative accuracy improving with increasing station separation. 

Although these results do not reflect the uncertainty due to 
ephemeris error, they include the uncertainty due to a five percent 
error in tropospheric refraction and a measurement uncertainty of 3 cen- 
timeters. The ephemeris error will increase the uncertainties of the 
estimated parameters as demonstrated by the error magnitudes given in 
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Table 5.3.1. However even with such Increases the double differencing 
approach appears to be adequate for providing rapid first-order 
determination of baselines. 


6. SUMMARY AND RECOMMENDATIONS 


6 . 1 Precision Comparison 

Some general conclusions were draw In Chapter 3 from an exami- 
nation of the results. For the observation types considered it was 
evident that ranging measurements provided the best geometric strength of 
solution. The two other derived observation types, correlated range 
difference and interferometry , were geometrically weaker although the 
results obtained from these latter procedures can be greatly improved 
upon by increasing the observational precision. Correlated range dif- 
ference observations had best geometric strength when observed satellites 
were tracked over longer time intervals. With this type of tracking 
procedure both the baseline component and chord length uncertainties 
were minimized. For range and interferometric observations shorter 
satellite tracking intervals produced the least uncertainty in the base- 
line parameters. Lengthening the tracking interval for these observa- 
tion types Increased the resulting parameter uncertainties. However the 
rate of increase was smaller than the variation produced in the Doppler 
results by decreasing the satellite tracking interval. And finally the 
interferometry approach became geometrically weaker as the baseline 
length Increased to become a more significant percentage of the distance 
to the satellite; although, the relative error in parts per million, 
decreased for the baselines considered. 
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The analysis presented In Chapter 3 considered the relative geo- 
metric strength of three observation types, two derived from basic 
ranging. The results were based on the assumptions that satellite posi- 
tions in space were known and that the basic ranging measurements were 
subject to uncorrelated stationary random noise. 

6 . 2 Dynamic Point Positionin g 

The range observation results presented In Chapter 4 Indicate 
that such observations from GPS satellites can provide geodetic coordi- 
nates to an accuracy of approximately 85 to 125 centimeters using 
twenty-four hours of continuous observation. Those results were based 
on the use of a one-hour tracking interval, selecting satellites which 
provide the best geometric strength for the solution. If a longer site 
occupation period is utilized, then the uncertainty in the geodetic 
coordinates can be reduced further to approximately 35 to 65 centimeters 
after five days of observation. Since the majority of satellite passes 
are north-south, the estimated latitude has a smaller standard error 
than longitude and height except for stations located toward the poles. 
For these latter stations height uncertainty tends to be smaller since 
a higher frequency of lower elevation observations provide a better 
separation of height and timing errors. If a longer tracking interval 
is utilized, larger uncertainties in estimated position are to be 
expected since the effects systematic satellite position error will 
not average as rapidly. The dominant error source limiting the accuracy 
of geodetic coordinates is this error in satellite position. Thus 
improvement in the receiver noise level and in measurement or modeling 


252 


of tropospheric refraction will yield only minor Increases In accuracy. 
Therefore the geodetic utilization of GPS range observations In a 
dynamic point positioning approach will require satellite ephemerldes 
to be estimated as accurately ns possible. 

Integrated Doppler observations based on independent counts can 
be expected to yield geodetic coordinate uncertainties ranging from 9\ 
to 150 centimeters after twenty-four hours of observation. The uncer- 
tainties will diminish to 45 to 65 centimeters after an additional four 
days of observation. These resul ■ based on a one-hour tracking 
Interval with an expected receiver noise level of 3 centimeters, An 
Increase In the tracking interval to three hours produces a substantial 
increase In the geodetic coordinate uncertainties. Thus, as with 
ranging, the. best procedure Is to track satellites over short intervals 
to obtain stronger geometric strength of solution. Increasing the pre- 
cision of the Doppler receiver or the accuracy of tropospheric refrac- 
tion prediction will produce only a minor change in the results. Again 
the uncertainty Introduced into station position by ephemeris error 
dominates the effects of all other error sources. Reduction of the 
ephemeris error by fifty percent produces a decrease in position uncer- 
tainty of approximately 45 percent. Therefore precise ephemeris compu- 
tation will be required for accurate geodetic positioning using GPS 
Doppler observations . 

The major conclusion which can be stated regarding dynamic point 
positioning using range and Doppler observations from a Global Posi- 
tioning System of navigation satellites Is that the accuracy of esti- 
mated geodetic coordinates will be comparable with the results 
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obtainable with Transit Doppler observations. No major increase in 
accuracy can be anticipated. Thus replacement of the Transit System 
with a Global Positioning System will not be detrimental to the geode- 
tic community since with proper electronic receivers similar levels of 
performance can be expected. The GPS system does offer a distinct 
advantage. This system provides continuous observation thereby 
decreasing the interval of time required to obtain comparable results 
with the Transit System enabling satellite surveying to become a more 
efficient operation. Table 6.2,1 summarizes the effect of systematic 
and random error sources on dynamic point positioning. 

TABLE 6.2.1. EFFECT OF SYSTEMATIC AND RANDOM ERROR SOURCES ON DYNAMIC 
POINT POSITIONING USING ONE-HOUR SATELLITE TRACKING 
INTERVALS 


APPROXIMATE COORDINATE ERROR (cm) 


ERROR SOURCE 

RANGE 

DOPPLER 


1 DAY 

5 DAYS 

1 DAY 

5 DAYS 

TROPOSPHERIC REFRACTION 

10 

5 

10 

5 

EPHEMERIS 

50*80 

25-40 

60-150 

30-70 

RESIDUAL SATELLITE RUBIDIUM 
CLOCK ERROR 

4 

1 

5 

2 

RECEIVER CESIUM CLOCK ERROR 

5 

2 

7 

3 

RECEIVER WHITE NOISE 
(RANGE 1m, DOPPLER 3 cm) 

2 

1 

IB 

B 


6.3 Baseline Determination 


Simultaneous range observations from two stations were utilized 
to determine baseline components and chord length. Solutions based on 
a one-hour tracking interval, selecting satellites which provide the best 
geometry, indicated after one day of observation that the lati- 
tude and longitude components of the baseline have uncertainties 
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of from 10 to 17 centimeters for baselines under 300 kilometers. The 
latitude component was determined with greater accuracy because of the 
frequency of north and south going satellite passes. For baselines 
under 300 kilometers the uncertainty of the height component ranged from 25 
to 28 centimeters. After five days of observation uncertainties in the 
latitude and longitude components were reduced to approximately k to 7 
centimeters and the height component uncertainty to 12 to 15 centi- 
meters. The uncertainty of these components Increased with baseline 
distance reflecting an increasing projection of orbit uncertainty into 
the estimates. For shorter baselines the uncertainty in chord length 
ranged from 10 to 16 centimeters after one day of observation and from 
4 to 7 centimeters after five days. The uncertainty was less for 
north-south baselines and Increased with station separation. However 
for the baselines considered here the relative uncertainty or ratio of 
the uncertainty in the chord to its length decreased with Increased 
baseline distance. The accuracy of 100 kilometer baselines was approxi- 
mately 1 to 1.5 parts per million after one day of observation. An 
increase in the length of the satellite tracking interval slightly 
degraded these results. 

The dominant error sources which will effect the accuracy of 
baseline determination using range observations are the stability of the 
receiver clock and error in refraction prediction. Increasing the 
modeling accuracy of tropospheric refraction will significantly 
increase the accuracy of the baseline parameters. Reducing the receiver 
noise level from 1 meter to 60 centimeters will produce a marginal 
increase in accuracy. 
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The stability of the receiver clock can be Improved by using an 
atomic oscillator with better stability properties. However the oscil- 
lator model chosen for use In this study was typical of cesium oscilla- 
tors having good stability; thus, It Is anticipated that the baseline 
parameter uncertainties attributed to random receiver clock error In 
this study are typical of those expected for an operational survey sys- 
tem. 

With Doppler observations from one-hour satellite tracking 
Intervals the uncertainties in the baseline components ranged from 27 to 
50 centimeters after one day and from 13 to 19 centimeters after five 
days of continuous observation. The chord length uncertainty ranged 
from 28 to 44 centimeters after one day and from 13 to 19 centimeters 
after five days. The latitude component of the baseline was determined 
with the least uncertainty and the chord lengths of north-south 
baselines were determined significantly better. These results are for 
baselines under 200 kilometers and are based on a 3 centimeter receiver 
white noise standard error. Increasing the tracking Interval to three 
hours produced some increase In the accuracy of the chord but the 
results appeared mixed for the baseline component uncertainties . 

The accuracy of the baseline parameters obtained by the geo- 
metrically weaker Doppler observations are improved significantly by 
decreasing the receiver noise level to 1 centimeter. Below that level 
the clock error statistics dominate and further increased precision will 
yield only marginal Improvement, Enhanced modeling or measurement of 
tropospheric refraction would improve the Doppler results but not as 
significantly as for ranging. 
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Assuming an instrumental noise level oil 1 centimeter, uncer- 
tainties In the baseline components would be reduced to approximately 9 
to 14 centimeters after five days of observation. The uncertainty in 
the chord would also be approximately 9 to 14 centimeters after five 
days for baselines under 200 kilometers. The errors limiting the 
accuracy of baseline determination using GPS Doppler observations are 
receiver noise and the stability of the receiver oscillator, 

Simultaneous Interferometric phase observations from two sites, 
twice differenced to eliminate timing errors , were examined as an alter- 
native procedure for the determination of baseline components. The use 
of continuous wave phase measurements requires the Introduction of inte- 
ger unknowns Into the adjustment related to the ambiguity in recognizing 
the exact cycle on which phase measurements were made at the two sites. 
As a consequence Initial and secondary adjustments of the baseline para- 
meters are required. 

Results obtained using a two-hour tracking Interval with a phase 
measurement uncertainty of 3 centimeters revealed that baseline compo- 
nents may be recovered with an uncertainty of from 1.0 to 4.0 centi- 
meters after six hours If sites are separated by up to a few hundred 
kilometers. The uncertainty In the recovered height difference between 
observing sites was larger than the uncertainties in the latitude and 
longicude differences for baselines under a few hundred kilometers. The 
accuracy of the chord length exceeded 0.1 ppm In all cases considered and 
improved with station separation. These results included uncertainty 
due to a five percent error in tropospheric refraction. Probable 
ephemerls error will Increase the uncertainty of the baseline components 
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no indicated in Table 6.3.1; however, even with ouch increaoeo thin 
approach appears to be adequate for rapid first-order determination of 
baselines under 200 kilometers. Table 6,3.1 summarizes the effect of 
various systematic and random error sources on baseline component deter- 
mination, 

6 * 4 Recommendations 

The results presented in Chapter 4 demonstrated that GPS range and 
integrated Doppler observations will provide sufficient accuracy for 
estimation of geodetic coordinates. These observations taken simul- 
taneously at two sites can be utilized to determine baseline parameters 
to better than 15 centimeters after five days of observation. A 
limiting factor for both observational approaches is the inability of 
the receiver oscillator. For certain geodynomlc applications such as 
earthquake prediction accuracies of 10 centimeters or better may be 
required within a short time Interval. GPS range and Doppler observa- 
tions might be capable of providing such accuracies in the future but 
the time Interval required to obtain such results will preclude this 
application. 

Satellite interferometry techniques can be developed which cir- 
cumvent the requirements for high stability frequency standards. This 
lead to the examination of the double differencing of interferometric 
phase. Thus one limiting factor for the range and Doppler approaches Is 
theoretically not a critical limitation for Interferometry, 

There are several Interferometric approaches which have been 
proposed using GPS satellites as radio sources. These proposals 
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have been described In this study and the error sources associated with 
satellite Interferometry have been mentioned. It Is recommended that a 
detailed error analysis of these Interferometry proposals be made to 
determine their effectiveness for determining baseline components. The 
specific details of each need to be examined so that a fair comparison 
Is realized. Further consideration should be given to the long-term 
cost effectiveness of these proposals Including the range and Doppler 
Instrumentation utilized in dynamic point positioning. 
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APPENDIX A 


LEAST SQUARES POLYNOMIAL APPROXIMATION 
OF RANDOM WALK SEGMENTS 


A.l General Polynomial Approximation 

Let (u k ) be a discrete stationary zero-mean stochastic process 
and define {z n )to be its running sum with 

SL 

Z l “ ** U k ‘ (A. 1.1) 

k = I 

The quantity z^ is one element in the discrete random walk sequence 
{z^}. By stationary it is meant that the random process {u^} is one 
whose statistical properties are invariant in time. Further, assume 
that over selected time intervals the random walk {z^} appears to be 

th 

dominated by systematic components enabling {z^ } to be modeled by an m 
degree polynomial P m (t). The difference between z^ and r^(t^) will be 
called the residual r^i 

T Z ^ z l ~ ^ = (A.l. 2) 

where the polynomial model is defined by 

m 

P (t) = Z a . (t - t ) J . (A. 1.3) 

ra 4 = „ 3 0 
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The constant: t In equation (A. 1.3) is arbitrary. The coefficients of 
the approximating polynomial can be determined by a least squares fit of 
P (t) to the random walk elements z sampled within a selected time 
Interval. The least squares solution for this approximation Is 
[Uotlla, 1967] 


_ T -3 T— 
a 13 (A A) A z 


(A. 1.4) 


where 


-T 


[ v 


V 


»a 


m 


— T 

2 a ^ 2 ’ ' ’ ’ ’ . 


The design matrix A is given by 



r**J 

(t i- 



w" 

> 

I! 

1 

/— 1 v 

rr 

to 

I 

V 



1 

( V 





(A. 1.5) 


The covariance for the polynomial coefficients a depends on the choice 
of t Q . In terms of the underlying process {u^}, it i'* given by the 
following equations 


E[&) = (A T A) _1 A T E[z] = 0 

— — T T — 1 T „T T —1 

Etaa X ] « (A A) A E[zz ]A(A A) 

— T 

where the covariance E[zz ] is given by 


(A. 1.6) 


E[zz^j • [R S]E[uu^ 



(A. 1.7) 
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and 


The matrices R and S are given below. 

Notice that the fitting procedure, equation (A. 1,4), is non- 
weighted least squares. The problem considered here is one of approxi- 
mation, not linear estimation, since the a^ are based on samples from 
{z^} not subject to an observation or sampling error. Also the proce- 
dure is Independent of how the z^ are selected within the time interval. 

Using equations (A. 1.2), (A. 1.3), and (A. 1.4), the residual vec- 
tor r can be written as 

- - t -1 T- 

r a z - A(A A) A z 

» [I- A(A T A)“ 1 A T ]i (A. 1.8) 

= Gz . 

The residuals represent the discrepancy between the samples of the .ran- 
dom walk and the approximating polynomial and may be Interpreted as 

"noise" with respect to P (t) . The statistics of these residuals are 

m 

obtained from the statistics of the random walk by the linear transfor- 
mation 

E[rr T ] = GE[zz T ]G T . (A. 1.9) 

This equation is derived using equation (A. 1.8). For the 
residuals 

E[r] = E[Gz] = GE[z] = 0 (A. 1.10) 

since, using equation (A.i'.X), each z^ has zero mean. 
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Thus , 


E[rr T ] o E[Gzz T G T ] » GE[za T ]G T . 


'<ux.ll) 


Therefore given the statistics of {u^} and an m^ 1 degree polynomial 
model P^Ct) to approximate {z^} over a given Interval of time, the 
statistics of the residuals to that model may be developed. 


<£ 

Theorem A. l: The covariance E[rr ] is (i) Independent of the 

epoch of {u^} provided this underlying process is stationary, and is 
(11) Invariant provided the {z r } are sampled in an identical fashion in 
each of two Intervals with comparable polynomial models being adopted. 
The proof is as follows : 

Let 

-T _ t , 

Zj ~ t z^ , , . • . , Z^J 


and 


-T 

Z II 


^ Z N+1 ’ Z N+2 Z 2N^ 


be two identically sampled sequences of the random walk {z^}. Since 


Z N+1 Z N + U N+1 


'2N 


“N 


+ 


“n+i 


+ 


+ u 


2N 


(A. 1.12) 


where 


z 


N 


N 

Z 

k = 1 


*V ’ 


equation (A. 1.12) can be written using 
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as 


» SUj 


Z 11 " Ru I + Su II 


where R Is an N x N matrix of all ones and S Is an H x N 
lower-triangular matrix of ones. Equation (A. 1.13) can be 


z IX a [R S] 


“I 

U IIJ 


By a linear transformation, taking u ° (uj, u n) ’ 


EfZjjijj] =[RS]E[uu T ] 



_ _n> T - -T T 

=> REfUjUjjR + SEf.Ujj.UjjR 

+ RE[UjUjj]S T + SE[UjjUjj]S T . 


From equation (A. 1.9) 

E[r ii r ix- 1 " CE ^ Z IX Z IX^ G 

= GRE[UjUj]R T G T + GSE[UjjUj]R T G T 

rp m m m m tti 

+ GREfUjUjjjS G + GSEfijjjUjjjS G . 


However 


(A. 1.13) 

written as 

(A. 1.14) 

(A, 1.15) 
(A. 1.16) 
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GR => [X-A(A T A)” 1 A T ]R « [I - C]R 


N 

E G 

k « 1 


Ik 


N 

E G. 
k» 1 


Nk 


(A. 1.17) 


Since the coefficients C^, H « 1,2,...,N, are based on the least 
squares approximation of P^t) to the sampled {z^}, it is true that 
(see lemma below) 


N 

E C 
k a l 


fit 


(A. 1.18) 


for every %. Therefore 


GR « 0 


(A. 1.19) 


and equation (A. 1.16) reduces to 


“ GSEtu^^S^ 


(A. 1.20) 


Using the stationarity assumption on {u^} and the result that 


EtSjij] = SE^u^S* 


(A. 1.21) 


equation (A. 1.20) becomes 


EL^ii^i] = GEtz^JG 1 = Elr^] 


(A. 1.22) 


since the matrix G in either case is identical. Thus the quantity 

SE(u,^.jUjj.]S is the only partial sum of the z ^ statistics which is 
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mapped by G Into the statistics of r^. The additional terms In equa» 
tlon (A, 1.16) Involving the R matrix are mapped Into zero by G. 


Lemma A.l : For least squares polynomial approximation 

Z C, = 1 and Z G.. « 0 (A. 1.23) 

i lj i ^ 

where 

G » X - C 

T -1 T 

C » A(A A) V . 


The proof Is as follows: 

The coefficients of the approximating polynomial 


P (0 

m 


m 

Z 

£ = 0 


a » (t ' V 


% 


are determined through a least squares procedure. 
Thus 


and 


..T..-1.T- 
a *= (A A) A z 

- _ T -1 T- 

r=z-Aa=z - A(A A) A z 

« z - Cz ® Gz . 


Consider the matrix product AG 


t T m T -1 T 

A G » A (I-C) = A [I-A(A A) A ] 


T T 

A - A = 0 . 


(A. 1.24) 
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Since by definition the first row of A T consists only of ones because 
the approximating function is a polynomial, it follows from equation 
(A, 1.24) that 


£ G. , « 0(1 « row, j « column) 
1 1J 

for all j . Since G is symmetric 

J G « “ 0 

for all 1. From equation (A. 1.23) 


Z 

i 



= 1 ~ 



Therefore, 


Z G 
i 


lj 


» 1 . 


The above theorem also holds for random walks in which the 
underlying process is continuous. For instance if equation (A. 1.1) is 
replaced by the continuous random walk 


z N = z(t N ) = / u(t)dt , 

fc 0 


(A. 1.25) 


then using equation (A. 1.25), equation (A. 1.13) can be expressed as 
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h 


M to* 

t N+l 

f u(t)dt 


/ u(t)dt 

C 0 


**N 

t 


+ 

# 



C 2N 

/ u(t)dt 


/ u(t)dt 

*0 


> 


w Q + T . 


fch 

The residuals based on the m degree polynomial fit are 

r xx “ (I ~ c)z n 

« (I-C)(Q + T) 

“ (l-G)Q + (I-C)T . 

Since Q is a vector of equal constants and since equation 
as before, it is obvious that 

(I - C)Q « 0 


and therefore 


r n - (I-C)T 
= GT . 


The covariance for the second set of residuals is 


Etrjjr^] = GEfTT T ]G T 

since 


E[T] ** 0 . 


(A. 1.26) 


(A. 1.27) 


(A.1J8) holds 


(A. 1.28) 


(A. 1.29) 


(A. 1.30) 


(A. 1.31) 


269 


Since {u^} is assumed stationary with autocorrelation function 
R^Ct-t'*), the following integral equation Is valid 


t K+J t N+X C N+J " h C N+I “ H J 

/ / R (t- t")dtdtt" « / / R (t* - t*'')dt*dt* > 

t M t M U 0 0 U 

N N 

(A. 1.32) 

tj tj. 

- / / R (t* - t*'')dt*dt*'“ 

0 0 u 

where 

t* - t - t M 

- t N . 


Applying equation (A. 1.32) to each element of the covariance matrix 
E[TT T ], it Is seen that 

E[TT T ] o EtSjZj] (A. 1.33) 


and thus equation (A. 1.30) becomes 


E ^II*II^ * GEtZjZjjG 1 


« Etr^j] , 


(A. 1.34) 


Therefore the theorem is valid in the continuous case. 


A. 2 Correlation Between Residuals from Approximations 
t~ Successive Random Walk Segments 

Consider two elements of the random walk sequence {z^} 
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(A, 2.1) 


and 


N+i 

z N+i a \ 
k a l 


N 

Z N “ Z U k 
N lc«l K 


(A, 2. 2) 


From the assumption that {u^} is a zero-mean process 

E[z n+1 ] a E[z n ] » 0 . (A. 2. 3) 

The correlation between these elements is given by 

E| - z N+i “ E ^ z N Z N ] + E[u N+ 1 Z N ]+ + E ^ U N't*i V 

2 (A. 2. 4) 

=* + E[u^^ Zjj] + •*. + E[Ujj^ Zjj] . 

'N 

The correlations between z^ } and the elements u N+ j depend on the correla- 
tions among the elements of the underlying process {u^}. 

Now consider the following question. If samples or a segment of 
(z n ) are to be modeled by a polynomial of degree m, what correlations 
exist between the residuals from successively fitted segments? Consider 
for example the random quantities z(t^) where 

n 

z(t ) = Z u(t,) n “ 1, ,2& . 

n k=l k 


Suppose the distribution of the z^ is as follows: 
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(Segment 1) 


8(tj) 

* 


a(t A ) 

z(t A+D 

i 

» 

z(t 0 J 

$■*&<* 


(Segment 2) 


Let t" and t** represent, without loss of generality, the midpoint of 
each segment. Xf a polynomial is chosen to model the z R over each seg- 
ment as 


P (t) a Z a,(t~t'') 3 (Segment 1) 

j « o J 

and 

m 

Q (t) « l b.(t-0 J (Segment 2) 

j°o ] 


can it then be assumed that r^ and r^ are uncorrelated? To answer this 
question perform a least squares fit of P^C t) and Q m (t) simultaneously 
to Che {z^} segments and then linearly transform the statistics of {z^} 
to obtain the residual statistics and compare the correlation coefficients 
between the two groups of residuals with those within each group, itie 
equation for this transformation Is 


E [?!?*] 

ECr^j] 

n Q 

E[ Vi ] 

‘■■A 1 

E[r r^T 
M II I- 

- -T , 
L £ r J 


l«* [ Z *j. !jr Z ^ 
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where 


G ° [I~A(A T Ar 1 A T ] 


CA i to • h) 


and 


A « 


A 1 0 

0 A, 


(A, 2. 7) 


A comparison of the correlation coefficients 


p a 

K r.r 


E[r,r . ] 

, ** J 


t j >^ 1 r i Mr J r j j 


C ..2.?) 


of the off~dlagonal blocks, Etr^r^] or Etr^rj] , with those of the 
diagonal blocks, E[r T r^] and E[r TT riL], can be a basis for deciding if 

W m i 4* J* JL 

the sets of residuals may be assumed to be Independent of not. 
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